Spectroscopic polarimetry

ABSTRACT

In the channeled spectroscopic polarimetry, a measurement error of a parameter showing a spectropolarization characteristic of a sample is effectively removed, the error being generated by various variations in retardation of a retarder depending upon the state of the sample. With attention being focused on that the retardation of the retarder may be kept constant by stabilization of an incident direction of light that transmits through the retarder, the retarder was arranged on the light source side with respect to the sample so as to effectively remove an influence relative to a measurement error, such as variations in direction of a light ray due to the sample.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to a technique for stabilizing measurement of spectropolarization characteristics of an object to be measured by use of a channeled spectrum.

2. Description of the Related Art

Light has properties of a “transverse wave”. Based upon the premise of three mutually orthogonal axes (x, y, z), when a propagation direction of light is assumed to be the z-axis direction, a vibration direction of the light is a direction along the x-y plane. The vibration direction of the light within the x-y plane has a bias. This bias of light is referred to as “polarization”. A biased state of light is referred to as a “state of polarization (SOP)” in this specification. Typically, the SOP varies depending upon wavelengths (colors) of light.

When light in some state of polarization is incident on an object to be measured to acquire outgoing light such as transparent light or reflected light and the object to be measured has optical anisotropy, a change in SOP is observed between the incident light and the outgoing light. Acquiring information on anisotropy of the object to be measured from the change in SOP is referred to as “polarimetry”. It is to be noted that causes of such anisotropy may include anisotropy of a molecular structure, presence of stress (pressure), and presence of a local field and a magnetic field.

A measurement in which a change in SOP between the incident light and the outgoing light is obtained with respect to each wavelength and information on anisotropy of an object to be measured is then acquired is especially referred to as “spectroscopic polarimetry”. This spectroscopic polarimetry has an advantage of acquiring a great amount of information as compared to the case of measurement by use of a single wavelength (single color). In the spectroscopic polarimetry, a device for measuring the change in SOP between the incident light and the outgoing light, namely a spectroscopic polarimeter, serves as a key device.

As fields of application of the spectroscopic polarimetry, there are known the field of spectroscopic ellipsometry, the medical field, and the like. In the field of spectroscopic ellipsometry, for example, since thickness as well as a complex refractive index of a thin film can be measured in a nondestructive and non-contact manner, spectroscopic polarimetry has been applied to optical electronic devices, examination/study of semiconductors, and the like. In the medical field, attempts have been made for early detection of glaucoma and a cancer cell since several kinds of cells have polarization characteristics.

As conventional typical spectroscopic polarimetry, rotating-retarder polarimetry and polarization-modulation polarimetry are known.

In these methods, a mechanical or electric polarization control element is used to modulate light to be measured, and from a change in spectral with the modulation, a Stokes parameter or the like is obtained.

However, the following problems and some other problems with the above methods have been pointed out: [1] a mechanical or electric drive unit is required; and [2] it is necessary to repeatedly change a plurality of spectrums while changing conditions of the spectroscopic polarimetry.

In order to solve these problems, channeled spectropolarimetry was previously contrived (cf. “Measurement of spectral distribution of polarized light based on frequency region interference method”, written by Takayuki Katoh, Kazuhiko Oka, Tetsu Tanaka, and Yoshihiro Ohtsuka, preliminary manuscript collection for 34th Academic Lecture Meeting of Hokkaido Branch of Japan Society of Applied Physics, (Hokkaido Branch of Japan Society of Applied Physics, Sapporo, 1998) p. 41).

Further, spectroscopic ellipsometry utilizing the channeled spectropolarimetry has also been reported (cf. “Spectroscopic ellipsometry using channeled spectrum”, written by Kazuhiko Oka and Takayuki Katoh, collected papers of lectures in 26th Study Session on Light Wave Sensing Technology (Light wave Sensing Technology Study Session held by Japan Society of Applied Physics, Dec. 19-20, 2000) pp. 107-114).

A configuration view of an experiment system for explaining the channeled spectroscopic polarimetry is shown in FIG. 26. As apparent from this figure, white light emitted from a xenon lamp 1 is transmitted through a polarizer 2 and a Babinet-Soleil compensator 3, to obtain a light wave having an SOP depending upon a frequency ν. Spectral distributions S₀ (ν), S₁ (ν), S₂ (ν) and S₃ (ν) of the Stokes parameters of the light wave are obtained by a measurement system 4 surrounded by a dashed line in the figure.

Light under measurement is first transmitted through two retarders R1 and R2 having different thicknesses (d1, d2) and an analyzer A, and then incident on a spectrometer 5. Here, a slow axis of the retarder R2 is inclined at an angle of 45° with respect to a slow axis of the retarder R1, while a transmission axis of the analyzer A is arranged in parallel to the slow axis of the retarder R1.

In each of the two retarders R1 and R2, a phase difference created between the orthogonal polarization components depends upon a frequency. Hence, as shown in FIG. 27, a channeled spectrum including three carrier components is obtained from the spectrometer 5 which functions as an optical spectrum analyzer. An amplitude and a phase of each of the carrier components are modulated by the spectrum distribution of the Stokes Parameters of the light under measurement. It is therefore possible to obtain each of the Stokes Parameters by execution of a signal processing with a computer 6 by use of Fourier transformation.

One example of results of an experiment is shown in FIG. 28. This is a result obtained in the case of inclining the Babinet-Soleil compensator 3 at an angle of 30° with respect to the slow axis of the retarder R1. Three solid lines respectively show spectral distributions S₁(ν)/S₀(ν), S₂(ν)/S₀(ν), S₃(ν)/S₀(ν) of the standardized Stokes parameters. It is thereby understood that the SOP depends upon the frequency.

As thus described, according to the channeled spectroscopic polarimetry, it is possible to obtain each spectrally-resolved Stokes Parameter by a frequency analysis (or wavenumber analysis) of characteristics of a spectral intensity. It is reasonably necessary to obtain respective retardations of the two retarders R1 and R2 prior to the frequency analysis. Here, the retardation means a phase difference created between a fast axis component and a slow axis component.

According to the foregoing channeled spectroscopic polarimetry, advantages can be obtained including that: [1] a mechanically movable element such as a rotating retarder is unnecessary; [2] an active element such as an electro-optic modulator is unnecessary; [3] four Stokes Parameters are obtained at once with one spectrum so that a so-called snap shot measurement can be performed; and [4] the configuration is simple, and thus suitable for size reduction.

SUMMARY OF THE INVENTION

However, concerning the foregoing channeled spectroscopic polarimetry, a problem of a relatively large measurement error has been pointed out for the following reasons.

In measurement of a spectropolarization characteristic by the channeled spectropolarimetry, it is necessary to previously calibrate retardation of a retarder. However, when the incident direction of light to be incident on the retarder changes between the time of calibration and the time of sample measurement, a distance of light to pass through the retarder changes, and thereby the retardation changes. This change in retardation between the time of calibration and the time of measurement has caused a measurement error. Further, it has been pointed out that, especially in cases including the case of using a higher-order retarder, retardation widely varies due to variations in light ray direction, or wave surface fluctuation, of light to pass through the retarder, or the like.

Moreover, there are mainly two kinds of methods as follows for investigating properties of an unknown sample by use of the channeled spectropolarimetry: [A] light is reflected on the sample, and by use of an SOP of light acquired from the reflected light, the properties of the sample is investigated; and [B] light is transmitted through the sample, and by use of an SOP of light acquired from the transmitted light, the properties of the light are investigated. The foregoing retardation variation can also be seen in each of those cases. In the following, each of those cases is described.

[A] When light is reflected on the sample so that a spectropolarization characteristic of the sample is measured, it is necessary to keep the incident direction of the wave surface of light to be incident on the retarder constant between the time of pre-calibration and the time of measurement. However, as also apparent from FIG. 24, an incident angle of the light to be incident on the sample varies in many ways due to variations in state of the surface of each sample or arrangement position of the sample, or the like, resulting in variations in incident direction of the wave surface of the light to be incident on the sample, whereby it becomes difficult to keep the retardation of the retarder constant from the time of pre-calibration. It should be noted that in FIG. 24, reference symbol B denotes a sample, reference symbols R1 and R2 respectively denote a first retarder and a second retarder, reference symbol A denotes an analyzer, and an indicator indicates the traveling direction of the light.

[B] Also when light is transmitted through the sample so that a spectropolarization characteristic of the sample is measured, it is necessary to keep the incident direction of the wave surface of light to be incident on the retarder constant between the time of precalibration and the time of measurement. However, the incident direction of the wave surface of the light to be incident on the retarder varies in many waves caused by the variations in light ray direction (cf. FIG. 25A) due to the inclination characteristic of the sample (inclination of the sample surface), scattering of the light ray (cf. FIG. 25B) attributed to a physical characteristic such as a rough surface of the sample, or the like, whereby it becomes difficult to keep the retardation of the retarder constant from the time of pre-calibration. It should be noted that in FIG. 25, reference symbol C denotes a birefringent medium as the sample, reference symbols R1 and R2 respectively denote the first retarder and the second retarder, reference symbol A denotes the analyzer, and an indicator indicates the traveling direction of the light.

The present invention was made by noting the problems as described above, and has an object to solve the problem of variations in incident angle of a retarder seen in the conventional channeled spectroscopic polarimetry, to provide a channeled spectroscopic polarimetry and spectroscopic polarimeter which are capable of measurement with even higher accuracy.

Further objects and action effects of the present invention are readily understood by the skilled in the art by referring to the following description of the specification.

(1) A spectroscopic polarimetry of the present invention includes the steps of: preparing an object to be measured; preparing a polarimetric spectroscope; and obtaining the spectral intensity of the object to be measured by use of the polarimetric spectroscope.

Here, the polarimetric spectroscope includes: a projection optical system; an analyzer for allowing light to transmit therethrough, the light having been emitted from the projection optical system and reflected on or transmitted through the object to be measured; and a means of obtaining the spectral intensity of the light having transmitted through the analyzer. The projection optical system comprises a light source, a polarizer and a plurality of retarders, where the light source, the polarizer and the plurality of retarders are arranged such that light emitted from the light source is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order.

Here, the “object to be measured” is a generic term of an object placed in a light path between the projection optical system and the analyzer. Namely, other than a sample having an unknown spectropolarization characteristic that is intended to be an object of the spectroscopic polarimetry, a polarization element having a known spectropolarization characteristic such as a phase compensator is also included in the “object to be measured” when placed in the light path between the projection optical system and the analyzer.

“A plurality of retarders” include: a retarder which is arranged behind the polarizer with respect to the traveling direction of light and whose principal axis is oriented differently from the transmission axis of the polarizer; and another retarder whose principal axis is oriented differently from the principal axis of the above-mentioned retarder. The “analyzer is an optical element showing a different transmittance from that of the polarization component located in an orthogonal direction to the optical element. The “analyzer” is not restricted to have a plate shape or a film shape. For example, a polarization beam splitter is usable as the “analyzer”.

The “means of obtaining the spectral intensity” may be using a spectrometer or using the light source where a wavelength is scanned. The spectrometer in the case of using the light source where the wavelength is scanned may be one capable of detecting an amount of light received, and a timing for detecting the amount of light received is corresponded to the wavelength of light.

According to the spectroscopic polarimetry of the present invention, since the object to be measured has no influence on the orientation of light to transmit through the retarder, it is possible to perform the spectroscopic polarimetry with high stability.

(2) The spectroscopic polarimetry of the present invention may include a step of obtaining at least one of spectropolarization parameters of the object to be measured by use of the obtained spectral intensity.

In the specification, the “spectropolarization parameter” is used in the meaning of a parameter that expresses the spectropolarization characteristic of the object to be measured.

(3) The plurality of retarders that the projection optical system comprises may be a first retarder and a second retarder. In this case, each of the elements of the projection optical system is arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder.

(4) In the following, three techniques for obtaining the spectropolarization parameter in the case of using two retarders are described. A first technique is one comprising: obtaining, from the obtained spectral intensity, a spectral intensity component (first spectral intensity component) which nonperiodically vibrates with wavenumber and a spectral intensity component (third spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₂(σ)) of the second retarder and not depending upon the retardation (φ₁(σ)) of the first retarder, with wavenumber; and obtaining at least one of spectropolarization parameters by use of each of the spectral intensity components.

According to this method, it is possible to obtain a rate of change in amplitude rate between linearly polarized light components which are orthogonal to each other along the orientation of the principal axis of the first retarder, an intensity attenuation rate of each component, and the like. Further, the optical disposition in this case may be in either a reflection mode or a transmission mode. Namely, light allowed to transmit through the analyzer may be light emitted from the projection optical system and reflected on the object to be measured, light that transmitted through the object to be measured, or light scattered by the object to be measured. Examples of the rate of change in amplitude rate which can be obtained by this method include an arctangent ψ(σ) of a rate of change in amplitude rate which is one of the ellipsometric parameters, and a rate of change in amplitude rate caused by scattering of light due to particles.

(5) A second technique is one comprising: obtaining, from the obtained spectral intensity, at least one of a spectral intensity component (second spectral intensity component) which vibrates at a frequency depending upon a difference between a retardation (φ₁(σ)) of the first retarder and a retardation (φ₂(σ)) of the second retarder with wavenumber, a spectral intensity component (fourth spectral intensity component) which vibrates at a frequency depending upon the sum of the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, and a spectral intensity component (fifth spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₁(σ)) of the first retarder and not depending upon the retardation (φ₂(σ)) of the second retarder, with wavenumber; and obtaining at least one of the spectropolarization parameters of the object to be measured by use of the obtained spectral intensity component.

According to this method, it is possible to obtain a variation in phase difference between linearly polarized light components which are orthogonal to each other along the orientation of the principal axis of the first retarder, and the like. Further, the optical disposition in this case may be in either a reflection mode or a transmission mode. Namely, light allowed to transmit through the analyzer may be light emitted from the projection optical system and reflected on the object to be measured, light that transmitted through the object to be measured, or light scattered by the object to be measured. Examples of the rate of change in phase difference which can be obtained by this method include a phase difference variation Δ(σ) which is one of the ellipsometric parameters, and a variation in phase difference caused by scattering of light due to particles.

(6) A third technique is one comprising: obtaining, from the obtained spectral intensity, at least one of the spectral intensity component (first spectral intensity component) which nonperiodically vibrates with wavenumber and the spectral intensity component (third spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₂(σ)) of the second retarder and not depending upon the retardation (φ₁(σ)) of the first retarder, with wavenumber, and at least one of the spectral intensity component (second spectral intensity component) which vibrates at a frequency depending upon the difference between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, the spectral intensity component (fourth spectral intensity component) which vibrates at a frequency depending upon the sum of the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, and the spectral intensity component (fifth spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₁(σ)) of the first retarder and not depending upon the retardation (φ₂(σ)) of the second retarder, with wavenumber; and obtaining at least one of the spectropolarization parameters of the object to be measured by use of each of the obtained spectral intensity components.

Examples of the spectropolarization parameter that can be obtained by this method include the rate of change in amplitude rate between linearly polarized light components which are orthogonal to each other along the orientation of the principal axis of the first retarder, such as an arctangent ψ(σ) of the rate of change in amplitude rate which is one of the ellipsometric parameters and a rate of change in amplitude rate caused by scattering of light due to particles. Further, the optical disposition in this case can be in either a reflection mode or a transmission mode. Namely, light allowed to transmit through the analyzer may be light emitted from the projection optical system and reflected on the object to be measured, light that transmitted through the object to be measured, or light scattered by the object to be measured. scattering of light due to particles. Other examples of the spectropolarization parameter that can be obtained by this method include an azimuth R of a birefringent medium and a retardation δ(σ).

(7) It is possible to calibrate the retardation of the retarder by use of light itself which is applied to the object to be measured and used for calibration. Meanwhile, it is also possible to calibrate the retardation of the retarder by use of the foregoing polarimetric spectroscope by means of light not applied to the object to be measured, or to separately calibrate the retardation of the retarder without use of the polarimetric spectroscope. In one case of calibrating the retardation of the retarder by use of light for use in measurement, the foregoing measurement method using the two retarders comprises: obtaining the retardation (φ₂(σ)) of the second retarder from the spectral intensity; and obtaining at least one of the spectropolarization parameters of the object to be measured by use of the obtained spectral intensity and the retardation (φ₂(σ)) of the second retarder.

Examples of the spectropolarization parameter that can be obtained by this method include the rate of change in amplitude rate between linearly polarized light components which are orthogonal to each other along the orientation of the principal axis of the first retarder, such as an arctangent ψ(σ) of a rate of a change in amplitude rate as one of the ellipsometric parameters, and a rate of change in amplitude rate caused by scattering of light due to particles.

(8) In another case of calibrating the retardation of the retarder by use of light for use in measurement, the foregoing measurement method using the two retarders includes a step of acquiring data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder, wherein the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder are obtained from the spectral intensity and the data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder, and at least one of the spectropolarization parameters of the object to be measured is then obtained by use of the obtained spectral intensity, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.

The “data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder” is for example a rate between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder for each wavelength.

“Obtaining the retardation” includes the case of obtaining a parameter equivalent to the retardation. In particular, obtaining a complex function including information as to the retardation corresponds to obtaining a parameter equivalent to the retardation.

According to this spectroscopic polarimetry, it is possible to effectively reduce a measurement error of the spectropolarization parameter caused by variations in retardation of the retarder due to a temperature change or other factors.

(9) The retardation of the retarder may be calibrated by use of a reference value for calibration of the retardation as well as light for use in measurement. In one case of calibrating the retardation in such a manner, the measurement method using the two retarders includes the steps of: acquiring data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder; and acquiring a reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and a reference value (φ₂ ^((i))(σ)) for calibration of retardation of the second retarder, wherein the retardation (φ₂ ^((i))(σ)) of the second retarder and the retardation variation (Δφ₂(σ)) of the second retarder from the reference value (φ₂ ^((i))(σ)) for calibration of the retardation are obtained from the obtained spectral intensity; the retardation variation (Δφ₁(σ)) of the first retarder is obtained by use of the obtained retardation variation (Δφ₂(σ)) of the second retarder and data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder; the retardation (φ₁(σ)) of the first retarder is obtained from the reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and the obtained retardation variation (Δφ₁(σ)) of the first retarder, and at least one of the spectropolarization parameters of the object to be measured is obtained by use of the obtained spectral intensity, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.

The “data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder” is for example a rate between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder for each wavelength. To this rate for each wavelength, it is possible to apply the rate between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder for each wavelength so long as a medium of the first retarder is the same as that of the second retarder.

The retardation of the second retarder obtained by use of light for use in measurement is accompanied with phase ambiguity by an integer multiple of 2π. Although this does not have an influence on a calculation error of the spectropolarization parameter by itself, phase unwrapping, which is performed in obtaining the retardation of the first retarder from the retardation of the second retarder, may cause a calculation error of retardation of the first retarder, thereby leading to generation of the calculation error of the spectropolarization parameter. Phase unwrapping is a process of determining a value of retardation of the second retarder such that the value of retardation of the second retarder continuously changes over the range of 2π with respect to a wavelength change. In the case of not using the retardation variation of the second retarder, the retardation of the first retarder is obtained by applying the “data showing the relation between the retardation of the first retarder and the retardation of the second retarder” to the retardation of the second retarder after phase unwrapping. When wavenumber intervals are not sufficiently large during change in value of retardation of the second retarder by 2π as compared to sampling wavenumber intervals, or when noise is included in the measured value of retardation of the second retarder, the retardation of the second retarder after phase unwrapping could be calculated by a wrong unit, 2π. If the retardation of the first retarder is obtained from the retardation of the second retarder including the error by the unit of 2π as thus described, the retardation of the first retarder typically includes an error not by the unit of 2π, and the included error becomes a large error when the spectropolarization parameter is calculated. According to the method for obtaining the retardation variation of the first retarder from the retardation variation of the second retarder, and then obtaining the retarder of the first retardation from the retardation variation of the first retarder and the reference value for calibration of the first retarder, since the retardation variation of the second retarder changes modestly with wavenumber, performing phase unwrapping on the retardation variation of the second retarder is unnecessary, or necessary only in a small frequency, which enables elimination of, or extreme reduction in, the possibility for generation of an error in the retardation variation of the first retarder due to phase unwrapping.

(10) In the spectroscopic polarimetry using the two retarders, the polarizer and the second retarder may be arranged such that an angle between the orientation of the transmission axis of the polarizer and the orientation of a fast axis of the second retarder is 45°.

In the case of arranging the polarizer and the second retarder such that the angle between the orientation of the transmission axis of the polarizer and the orientation of the fast axis of the second retarder is 45°, there is an advantage of simplifying calculation for obtaining the spectropolarization parameter. Meanwhile, in the case of not limiting the angle between the orientation of the transmission axis of the polarizer and the orientation of the fast axis of the second retarder to 45°, there is an advantage of facilitating manufacturing of the optical system since a limit on an assembly error of the optical system is alleviated.

(11) The spectroscopic polarimetry for obtaining at least one of the spectropolarization parameters further includes a step of obtaining a spectral intensity for calibration by use of the polarimetric spectroscope in a state where an object to be measured having an unknown spectropolarization characteristic does not exist in a light path between the projection optical system and the analyzer, wherein at least one of the spectropolarization parameters of the object to be measured is obtained by use of the obtained spectral intensity regarding the object to be measured and the spectral intensity for calibration or data based upon the spectral intensity for calibration.

In obtaining the spectral intensity for calibration, an object to change the spectropolarization state of light may be prevented from existing in the light path between the projection optical system and the analyzer, or an object having a known spectropolarization characteristic may exist in the light path.

(12) In obtaining the spectral intensity for calibration, an analyzer for calibration may be prepared in a position in which light emitted from the projection optical system is received in a state where the object to be measured having an unknown spectropolarization characteristic does not exist in the light path between the projection optical system and the analyzer, and a spectral intensity of light having passed through the analyzer for calibration may be obtained.

When the light path between the projection optical system and the analyzer is bent due to reflection on or inflection in the object to be measured, it is possible to obtain the spectral intensity for calibration by arrangement of the analyzer for calibration in a position in which light emitted from the projection optical system is received in a state where an object to bend the light path as in the case where the object to be measured exists in the light path does not exist. In this case, an object having a known spectropolarization characteristic may also exist in the light path. Further, the analyzer for calibration may be prepared separately from the analyzer for measurement, or the position of the analyzer for measurement may be temporarily changed so that the analyzer for measurement may be used as the analyzer for calibration.

(13) The spectroscopic polarimetry using the spectral intensity for calibration may include a step of obtaining the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder by use of the spectral intensity for calibration, wherein at least one of the spectropolarization parameters of the object to be measured may be obtained by use of the obtained spectral intensity regarding the object to be measured, the retardation (φ₁(σ)) of the first retarder, and the retardation (φ₂(σ)) of the second retarder, which are obtained by means of the spectral intensity for calibration.

(14) In the foregoing measurement method in which data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder is acquired by use of the two retarders and the retardation of the retarder is then calibrated by use of light for use in measurement, the spectral intensity for calibration may be obtained by use of the polarimetric spectroscope in a state where the object to be measured having an unknown spectropolarization characteristic does not exist in the light path between the projection optical system and the analyzer, and the data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder may be obtained by use of the obtained spectral intensity for calibration.

(15) In the foregoing measurement method in which data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder is acquired by use of the two retarders and the retardation of the retarder is then calibrated by use of a reference value for calibration of the retardation and light for use in measurement, the spectral intensity for calibration may be obtained by use of the polarimetric spectroscope in a state where the object to be measured having an unknown spectropolarization characteristic does not exist in the light path between the projection optical system and the analyzer, and the data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder may be obtained by use of the obtained spectral intensity for calibration.

(16) In the spectroscopic polarimetry of the present invention, a spectroscopic quasi-Stokes parameter of the object to be measured may be obtained by use of the obtained spectral intensity.

(17) In the spectroscopic polarimetry for obtaining the spectroscopic quasi-Stokes parameter, the plurality of retarders that the projection optical system comprises may be a first retarder and a second retarder. In this case, each of the elements of the projection optical system is arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of the transmission axis of the polarizer disagrees with the orientation of the principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of the principal axis of the first retarder. Further, the spectroscopic polarimetry may include a step of acquiring data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder. From the obtained spectral intensity, at least one of the spectral intensity component (first spectral intensity component) which nonperiodically vibrates with wavenumber and the spectral intensity component (third spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₂(σ)) of the second retarder and not depending upon the retardation (φ₁(σ)) of the first retarder, with wavenumber, and at least one of the spectral intensity component (second spectral intensity component) which vibrates at a frequency depending upon the difference between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, the spectral intensity component (fourth spectral intensity component) which vibrates at a frequency depending upon the sum of the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, and the spectral intensity component (fifth spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₁(σ)) of the first retarder and not depending upon the retardation (φ₂(σ)) of the second retarder, with wavenumber, may be obtained. The retardation (φ₁(σ)) of the first retarder, the retardation (φ₂(σ)) of the second retarder and the spectroscopic quasi-Stokes parameter may then be obtained by use of the data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder, and each of the obtained spectral intensity components.

The meaning of “obtaining the spectroscopic quasi-Stokes parameter” includes obtaining all or part of four spectroscopic quasi-Stokes parameters M₀, M₁, M₂, M₃ (respective definitonal equations for these parameters are provided in DETAILED DESCRIPTION OF THE INVENTION). Whether or not all the spectroscopic quasi-Stokes parameters are practically obtained is referred to a choice of a person implementing the spectroscopic polarimetry. However, according to the spectroscopic polarimetry of the present invention, it is possible in theory to obtain all the spectroscopic quasi-Stokes parameters.

For obtaining the spectroscopic quasi-Stokes parameter M₀(σ), the first spectral intensity component and a reference amplitude function m₀(σ) are needed. For obtaining the spectroscopic quasi-Stokes parameter M₁(σ), the third spectral intensity component, the retardation of the second retarder and a reference amplitude function m₂(σ) are needed.

For obtaining the spectroscopic quasi-Stokes parameter M₂(σ) and M₃(σ), at least any one of a group of the second spectral intensity component, the retardations of the first and second retarders and a reference amplitude function m⁻(σ), a group of the forth spectral intensity component, the retardations of the first and second retarders, and a reference amplitude function m₊(σ), and a group of the fifth spectral intensity component, the retardation of the first retarder and a reference amplitude function m₁(σ) is needed.

Further, the reference amplitude functions needed for obtaining the spectroscopic quasi-Stokes parameter need to be made usable when the spectroscopic quasi-Stokes parameter is obtained.

According to this spectroscopic polarimetry for obtaining the spectroscopic quasi-Stokes parameter, an active element such as a mechanically moving part or an electric optical modulator for polarization control is not necessary. In the spectroscopic polarimetry, by acquiring a one-time spectral, it is possible in theory to obtain all spectroscopic quasi-Stokes parameter of the object to be measured, and further possible to effectively reduce a measurement error of the spectroscopic quasi-Stokes parameter generated caused by variations in retardation of the retarder due to a temperature change or other factors. By performing further calculation by use of the spectroscopic quasi-Stokes parameters, it is possible to obtain a variety of spectropolarization parameters regarding the object to be measured. Especially in a case where a Mueller matrix of the object to be measured is determined from only two to three parameters at most, an arbitrary spectropolarization parameter can be obtained from the spectroscopic quasi-Stokes parameters.

(18) In the spectroscopic polarimetry for obtaining the spectroscopic quasi-Stokes parameter includes the steps of acquiring data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder by use of the two retarders; and acquiring a reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and a reference value (φ₂ ^((i))(σ)) for calibration of retardation of the second retarder. In the spectroscopic polarimetry, from the obtained spectral intensity, at least one of the spectral intensity component (first spectral intensity component) which nonperiodically vibrates with wavenumber and the spectral intensity component (third spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₂(σ)) of the second retarder and not depending upon the retardation (φ₁(σ)) of the first retarder, with wavenumber, and at least one of the spectral intensity component (second spectral intensity component) which vibrates at a frequency depending upon the difference between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, the spectral intensity component (fourth spectral intensity component) which vibrates at a frequency depending upon the sum of the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, and the spectral intensity component (fifth spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₁(σ)) of the first retarder and not depending upon the retardation (+2(σ)) of the second retarder, with wavenumber, may be obtained. The retardation (φ₂(σ)) of the second retarder and the retardation variation (Δφ₂(σ)) of the second retarder from the reference value (φ₂ ^((i))(σ) for calibration may be obtained by use of the obtained spectral intensity. The retardation variation (Δφ₁(σ)) of the first retarder may be obtained by use of the obtained retardation variation (Δφ₂(σ)) of the second retarder and data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder. The retardation (φ₁(σ)) of the first retarder may be obtained from the reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and the obtained retardation variation (Δφ₁(σ)) of the first retarder. The spectroscopic quasi-Stokes parameter may be obtained by use of each of the obtained spectral intensity components, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.

(19) A polarimetric spectroscope of the present invention comprises: a projection optical system, comprising a light source, a polarizer and a plurality of retarders, where the light source, the polarizer and the plurality of retarders are arranged such that light emitted from the light source is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order; an analyzer for allowing light to transmit therethrough, the light having been emitted from the projection optical system and reflected on or transmitted through the object to be measured; and a means of obtaining the spectral intensity of the light having transmitted through the analyzer.

According to this polarimetric spectroscope, the direction of light that transmits through the retarder is unsusceptible to the object to be measured, enabling spectroscopic polarimetry with high stability.

(20) In this polarimetric spectroscope, the plurality of retarders that the projection optical system comprises may be a first retarder and a second retarder. In this case, each of the elements of the projection optical system is arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder.

(21) In the polarimetric spectroscope using the two retarders, the polarizer and the second retarder may be arranged such that an angle between the orientation of the transmission axis of the polarizer and the orientation of a fast axis of the second retarder is 45°.

(22) The spectroscopic polarimetry of the present invention may comprise: an analyzer for calibration, detachably provided in a position in which light emitted from the projection optical system is received in a state where an object to be measured having an unknown spectropolarization characteristic does not exist in the light path between the projection optical system and the analyzer; and a means of obtaining the spectral intensity of the light having transmitted through the analyzer for calibration.

Here, the whole or part of the “means of obtaining the spectral intensity of the light having transmitted through the analyzer” may double as the “means of obtaining the spectral intensity of the light having transmitted through the analyzer for calibration”.

Also when the light path between the projection optical system and the analyzer is bent due to reflection on or inflection in the object to be measured, the use of this polarimetric spectroscope enables calibration even in a state where an object to bend the light path, in the same manner as in the case of existence of the object to be measured, does not exist in the light path. In this case, an object having a known spectropolarization characteristic may exist in the light path. It is therefore possible to perform calibration by singly using the device, prior to installation of the object to be measured in the device, or installation of the device on the object to be measured.

(23) The polarimetric spectroscope of the present invention may further comprise an optical fiber for projecting light which guides the light emitted from the light source to the polarizer.

According to this spectroscopic polarimetry, the light source can be installed in a position apart from the measurement place, thereby facilitating size reduction in a portion of the polarimetric spectroscope, the portion being used in the vicinity of the measurement place.

(24) In the polarimetric spectroscope of the present invention, the means of obtaining a spectral intensity may comprise a light-reception element or a spectrometer, and may further comprise an optical fiber for receiving light which guides the light having transmitted through the analyzer to the light-reception element or the spectrometer.

According to this polarimetric spectroscope, the light source can be installed in a position apart from the measurement place, thereby facilitating size reduction in a portion of the polarimetric spectroscope, the portion being used in the vicinity of the measurement place.

(25) A spectroscopic polarimeter of the present invention comprises the foregoing polarimetric spectroscope of the present invention and an arithmetic unit for obtaining at least one of spectropolarization parameters of an object to be measured, by use of a spectral intensity of light having transmitted through an analyzer.

(26) In the spectroscopic polarimeter of the present invention, the plurality of retarders that the projection optical system comprises may be a first retarder and a second retarder. In this case, each of the elements of the projection optical system is arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder. Further, the arithmetic unit of this spectroscopic polarimeter is made capable of using data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder. The retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder may be obtained from the spectral intensity of the light having transmitted through the analyzer and the data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder. At least one of the spectropolarization parameters of the object to be measured may then be obtained by use of the spectral intensity of the light having transmitted through the analyzer, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.

(27) In the spectroscopic polarimeter of the present invention, the plurality of retarders that the projection optical system comprises may be a first retarder and a second retarder. The arithmetic unit is made capable of using data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder, a reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and a reference value (φ₂ ^((i))(σ)) for calibration of retardation of the second retarder. The retardation (φ₂(σ)) of the second retarder and the retardation variation (Δφ₂(σ)) of the second retarder from the reference value (φ₂ ^((i))(σ)) for calibration may be obtained from the spectral intensity of the light having transmitted through the analyzer, the retardation variation (Δφ₁(σ)) of the first retarder may be obtained by use of the obtained retardation variation (Δφ₂(σ)) of the second retarder and data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder, the retardation (φ₁(σ)) of the first retarder may be obtained from the reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and the obtained retardation variation (Δφ₁(σ)) of the first retarder. At least one of the spectropolarization parameters of the object to be measured may then be obtained by use of the obtained spectral intensity, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.

(28) An optical device of the present invention comprises: a projection optical system, comprising a polarizer and a plurality of retarders, where the polarizer and the plurality of retarders are arranged such that light incident on the polarizer is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order; and an analyzer for allowing light to transmit therethrough, the light having been emitted from the projection optical system and reflected on or transmitted through the object to be measured.

Such an optical device can be used for the foregoing polarimetric spectroscope.

(29) A light-projection device of the present invention comprises a polarizer and a plurality of retarders, wherein the polarizer and the plurality of retarders are arranged such that light incident on the polarizer is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order.

Such a light-projection device can be used for the foregoing polarimetric spectroscope.

(30) In this light-projection device, the plurality of retarders may be a first retarder and a second retarder. Each of the elements of the light-projection device may be arranged such that light incident on the polarizer transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder.

Next described are a spectroscopic polarimetry, a polarimetric spectroscope and an optical device which make a characteristic of a polarization element or an azimuth angle of an analyzer changeable. The polarization element here refers to a polarization element in a case where the object to be measured is composed of a sample and a polarization element on which light, having transmitted through or been reflected on the sample, is incident. The polarization element is an optical element where incident light and outgoing light have a relation depending upon polarization. For changing the characteristic of the polarization element, it is possible to use, for example, a means of changing the azimuth angle of the polarization element, a means of changing the retardation of the polarization element, and some other means.

(31) In the spectroscopic polarimetry of the present invention, a polarimetric spectroscope further comprising a means of changing the azimuth angle of the analyzer may be used. By use of this polarimetric spectroscope, a spectral intensity regarding the object to be measured may be obtained in a plurality of states where azimuth angles of the analyzer are made different from one another, and at least one of the spectropolarization parameters of the object to be measured may be obtained by use of the obtained spectral intensity.

(32) In the spectroscopic polarimetry of the present invention, an object to be measured including a sample and a polarization element is prepared. The polarimetric spectroscope, further comprising a means of changing the characteristic of the polarization element, is used. By use of this spectroscopic polarimetry, a spectral intensity regarding the object to be measured may be obtained in a plurality of states where characteristics of the polarization element are made different from one another. At least one of spectropolarization parameters of the sample may then be obtained by use of the obtained spectral intensity.

(33) Here, the spectroscopic polarimetry may be used which further comprises a means of changing the azimuth angle of the analyzer in addition to the means of changing the characteristic of the polarization element. By use of this spectroscopic polarimetry, a spectral intensity regarding the object to be measured may be obtained in a plurality of states where characteristics of the polarization element or azimuth angles of the analyzer are made different from one another. At least one of spectropolarization parameters of the sample may then be obtained by use of the obtained spectral intensity.

(34) The polarimetric spectroscope of the present invention may further comprise a means of changing the azimuth angle of the analyzer. This polarimetric spectroscope may be combined with an arithmetic unit for obtaining at least one of spectropolarization parameters of the object to be measured by use of this spectroscopic polarimetry a spectral intensity regarding the object to be measured obtained in a plurality of states where azimuth angles of the analyzer are made different from one another, to give a spectroscopic polarimeter.

(35) In a case where the object to be measured includes a sample and a polarization element, the polarimetric spectroscope of the present invention may further comprise a means of changing the characteristic of the polarization element. This polarimetric spectroscope may be combined with an arithmetic unit for obtaining at least one of spectropolarization parameters of the sample by use of a spectral intensity regarding the object to be measured obtained in a plurality of states where characteristics of the polarization element are made different from one another, to give a spectroscopic polarimeter.

(36) Here, the spectroscopic polarimetry may further comprise a means of changing the azimuth angle of the analyzer in addition to the means of changing the characteristic of the polarization element. This polarimetric spectroscope may be combined with an arithmetic unit for obtaining at least one of spectropolarization parameters of the sample by use of a spectral intensity regarding the object to be measured obtained in a plurality of states where characteristics of the polarization element or azimuth angles of the analyzers are made different from one another, to give a spectroscopic polarimeter.

(37) The optical device of the present invention may further comprise a means of changing the azimuth angle of the analyzer.

(38) In a case where the object to be measured includes a sample and a polarization element, the optical device of the present invention may further comprise a means of changing the characteristic of the polarization element.

(39) Here, the optical device may further comprise a means of changing the azimuth angle of the analyzer in addition to the means of changing the characteristic of the polarization element.

As thus described, by making the characteristic of the polarization element or the azimuth angle of the analyzer changeable, it is possible to obtain a spectral intensity in a plurality of states where characteristics of the polarization element or azimuth angles of the analyzers are made different from one another. Thereby, from the spectral intensity obtained in a state where the number of characteristics of the polarizer or the number of azimuth angles of the analyzer is relatively small, a relatively many kinds of spectropolarization parameters can be obtained, and an error due to an influence such as noise included in the obtained value of the spectropolarization parameter can be reduced. Or, the characteristic of the polarizer or the azimuth angle of the analyzer can be selected such that a specific one or more than one kinds of spectropolarization parameters of the object to be measured or the sample can be obtained with high sensitivity.

According to the present invention, it is possible to perform the spectroscopic polarimetry with high stability since the direction of light that transmits through the retarder is not susceptible to the object to be measured.

BRIEF DESCRIPTION OF THE DPRWINGS

FIG. 1 shows a view for explaining a principle of solving a problem.

FIG. 2 shows a view for explaining a configuration of an optical system device and an azimuth angle of each optical element.

FIG. 3 shows a view for explaining spectroscopic polarimetry.

FIG. 4 shows a view for explaining a relation between a channeled spectrum obtained from a spectrometer and its four components.

FIG. 5 shows a view for explaining a procedure (flows of signal processing) for demodulating a spectrometric quasi-Stokes parameter.

FIG. 6 shows a view for explaining one example of Step 2.

FIG. 7 shows a view for explaining Fourier transformation.

FIG. 8 shows a view for explaining flows of signals for calibration during measurement.

FIG. 9 shows a view for explaining flows of signals in combination of “calibration during measurement” and “measurement of a spectrometric quasi-Stokes parameter”.

FIG. 10 shows a comparative view for explaining methods (No. 1, 2) for calibrating a reference phase function during measurement.

FIG. 11 shows a device configuration view (No. 1) in ellipsometry.

FIG. 12 shows a device configuration view (No. 2) in ellipsometry.

FIG. 13 shows a view for explaining a principle of solving a problem in the case of reflecting light on a sample.

FIG. 14 shows a view for explaining a device configuration in the case of separately installing an optical system for calibration.

FIG. 15 shows a view for explaining a polarization state of light reflected on the sample.

FIG. 16 shows a view for explaining a device configuration (No. 1) in a double refraction measurement.

FIG. 17 shows a view for explaining a device configuration (No. 2) in the double refraction measurement.

FIG. 18 shows a view for explaining a principle of solving the problem in the case of allowing light to transmit through the sample.

FIG. 19 shows a view for explaining a device configuration in the case of arranging a polarization element having a known spectropolarization characteristic in front of or behind the sample.

FIG. 20 shows a configuration view (No. 1) of one example of a spectroscopic polarimeter.

FIG. 21 shows a configuration view (No. 2) of one example of the spectroscopic polarimeter.

FIG. 22 shows a flowchart of a pre-calibration procedure.

FIG. 23 shows a flowchart of a measurement procedure.

FIG. 24 shows a view (No. 1) for explaining variations in wave surface of light that passes through a retarder in the incident direction.

FIG. 25 shows a view (No. 2) for explaining variations in wave surface of the light that passes through the retarder in the incident direction.

FIG. 26 shows a configuration view of an experimental system of channeled spectroscopic polarimetry.

FIG. 27 shows a graph of a Stokes parameter in the experimental system.

FIG. 28 shows a graph of a standardized Stokes parameter in the experimental system.

FIG. 29 shows a device configuration view in the case of measuring a spectropolarization parameter of the sample.

FIG. 30 shows a view of a configuration of a channeled spectroscopic polarization state generator (CSPSG).

DETAILED DESCRIPTION OF THE INVENTION

In the following, one preferred embodiment of the present invention is specifically described with reference to attached drawings (FIGS. 1 to 19).

Chapter 1: Principle of Channeled Spectroscopic Polarimetry

1.1 Constitution of Optical System in Present Invention

FIG. 1 shows a view for explanation comparing a configuration of an optical system in a conventional channeled spectroscopic polarimetry and a configuration of an optical system in a channeled spectroscopic polarimetry of the embodiment of the present invention. The optical system (cf. FIG. 1B) in the conventional channeled spectroscopic polarimetry is comprised of a light source 7, a polarizer P, and a polarimeter. The polarimeter is comprised of two thick retarders R1 and R2, an analyzer A and a spectrometer 8. It is to be noted that reference symbol D denotes a sample through which light is transmitted or on which light is reflected. Here, fast axes of the retarder R1 and the retarder R2 are inclined at an angle of 45° from each other. Meanwhile, a transmission axis of the analyzer A agrees with the fast axis of the retarder R1.

It is to be noted that crossing angles among these three elements may not necessarily be 45°. Measurement is possible even with a different crossing angle, although less efficient to some extent. In short, any crossing angle can be applied so long as principal axes of the adjacent elements are not superposed on each other. A description in this respect is given later. What is important is that each element is fixed and thus not required to be rotated or modulated as in the conventional method.

Light having a broad spectrum is emitted from the light source 7 on the left side of the figure, transmits through the polarizer and is reflected on or transmits through the sample D before being incident on the polarimeter. Thereafter, the light is incident on the polarimeter. A spectral distribution of state of polarization (SOP) of the light having been emitted from the sample D can be expressed by spectrometric Stokes parameters S₀(σ), S₁(σ), S₂(σ), and S₃(σ). Here, σ is a “wavenumber” defined by an inverse number of a wavelength λ. Further, coordinate axes x and y for determining the spectrometric Stokes parameters is taken so as to agree with the fast and late axes of the retarder R1.

The light incident on the spectrometer passes in this order through the retarders R1 and R2 and the analyzer A, and is incident on the polarimeter 8. The Stokes parameters depending upon the wavenumber a are obtained from a spectrum acquired from the polarimeter 8.

However, the optical system shown in FIG. 1B has caused a problem of a variety of changes in incident direction of the wave surface of the light that transmits through the retarder under the influence of the sample to result in generation of an error in spectroscopic polarimetry. The present invention is provided to solve such a problem.

The optical system of the embodiment of the present invention shown in FIG. 1A is comprised of the light source 7, the polarizer P, the retarders R2 and R1, the analyzer A and the spectroscope 8. The light emitted from the light source 7 transmits through the polarizer P, the retarder R2 and the retarder R1 in this order, and is reflected on or transmits through the sample D. The light then transmits through the analyzer A, and is incident on the spectroscope 8. Thereafter, a spectrum of the incident light is acquired in the spectroscope 8, and spectropolarization parameters of the sample and the like are calculated according to a later-described procedure.

As described above, in the present specification, the “spectropolarization parameter” is used in the meaning of a parameter for expressing the spectropolarization characteristic of the object to be measured. This is a generic name of parameters for use in quantitatively expressing a polarization change caused by reflection or transmission of light on or through the object to be measured. Examples of the spectropolarization parameter include ellipsometric parameters ψ(σ), Δ(σ), and a retardation δ(σ) of a double refraction medium. It should be noted that, although the spectropolarization characteristic of the object to be measured is typically expressed completely by 16 elements of a 4×4 Mueller matrix, there are few cases where all of those 16 elements are independent variables. In the spectroscopic polarimetry, all of those elements are often determined from only two to three parameters at most. In practice, those independent parameters may be obtained as the spectropolarization parameters. Moreover, there are applications in which simply obtaining part of the spectropolarization parameters satisfies the case regardless of whether the parameters are independent or non-independent.

Here, it is important that the retarders R2 and R1 have been arranged on the light-source side with respect to the sample D. This makes it possible to keep the incident direction of the wave surface of light to be incident on the retarder constant so as to realize highly-accurate, stable spectroscopic polarimetry. It should be noted that variations in incident direction of the wave surface of light to be incident on the analyzer A has almost no influence on the measurement result. This leads to solving the foregoing problem of the retardation change due to changes in distance and direction of a light ray that passes through the retarder during calibration of the retardation and during sample measurement.

Next, the embodiment of the present invention is specifically described with reference to FIG. 2. This optical system is comprised of the light source 7, the polarizer P, the retarders R2 and R1, the analyzer A and the spectroscope 8. It is to be noted that reference symbol D denotes a sample. Here, the orientations of the fast axes of the retarders R1 and R2 are inclined at an angle of 45° from each other. Meanwhile, the orientation of a transmission axis of the polarizer P agrees with the orientation of the fast axis of the retarder R1. In the figure, the fast axes of the retarders are denoted by “fast” and the slow axes thereof are denoted by “slow”. Further, θ is the azimuth angle of the transmission axis of the analyzer with respect to the fast axis of the retarder R1.

The Mueller matrix of the sample at this time is described as follows. [Mathematical Expression 1] $\begin{matrix} {{M(\sigma)} = \begin{bmatrix} {{\hat{m}}_{00}(\sigma)} & {{\hat{m}}_{01}(\sigma)} & {{\hat{m}}_{02}(\sigma)} & {{\hat{m}}_{03}(\sigma)} \\ {{\hat{m}}_{10}(\sigma)} & {{\hat{m}}_{11}(\sigma)} & {{\hat{m}}_{12}(\sigma)} & {{\hat{m}}_{13}(\sigma)} \\ {{\hat{m}}_{20}(\sigma)} & {{\hat{m}}_{21}(\sigma)} & {{\hat{m}}_{22}(\sigma)} & {{\hat{m}}_{23}(\sigma)} \\ {{\hat{m}}_{30}(\sigma)} & {{\hat{m}}_{31}(\sigma)} & {{\hat{m}}_{32}(\sigma)} & {{\hat{m}}_{33}(\sigma)} \end{bmatrix}} & (1.1) \end{matrix}$

Further, as parameters for effectively expressing a degree of polarization of light, an ellipticity angle, an azimuth angle, and the like, Stokes Parameters are used. The Stokes Parameters are composed of four parameters having definitions as follows:

S₀: total intensity

S₁: difference between intensities of linearly polarized light components with angles of 0° and 90°.

S₂: difference between intensities of linearly polarized light components with angles ±45°.

S₃: difference between intensities of left and right circularly polarized light components.

In a third-dimensional space where the three mutually orthogonal axes are taken as S₁, S₂ and S₃, assuming a sphere with a radius S₀ taking an original point of the axes as a center, an SOP of arbitrary light is expressed as one point in this third-dimensional space and a degree of polarization is expressed by the following expression: Degree of polarization=(distance from original point to point (S ₁ , S ₂ , S ₃))/S ₀=(S ₁ ² +S ₂ ² +S ₃ ²)^(1/2) /S ₀

Here, the Mueller matrix is described with reference to FIG. 3. The Mueller matrix is a matrix for expressing an interaction of light of reflection, transmission, etc. on or through a sample as an object to be measured. As an example, the following case is considered. Light whose SOP is expressed by S(σ) as a state of polarization 1 is incident on the sample, and comes under the influence of a polarization element and an object to be measured such as the sample. Light is then emitted whose SOP is expressed by S′(σ) as a state of polarization 2 (cf. FIG. 3A). Here, the Mueller matrix of the object to be measured is expressed by a relational expression as a 4×4 matrix shown in FIG. 3B (cf. FIG. 3B).

In the following described is a procedure for obtaining a spectropolarization parameter of the sample from the foregoing Mueller matrix or the like.

Before the procedure for obtaining a spectropolarization parameter of the sample is described, characteristics of the retarders R1 and R2 are formulated as preparation for the description. A retarder is an element having the property of changing a phase difference between mutually orthogonal linearly polarized light components before and after transmission of light through the element. An amount of such a change in phase difference is referred to as retardation. Further, coordinate axes taken along the two linearly polarized light directions are referred to as principal axes. Among them, the axis along the linearly polarized light whose phase relatively moves fast is referred to as a fast axis, and the other axis is referred to as a slow axis.

The retardation of a retarder R_(j) (j=1,2) made of a double refraction medium changes with wavenumber a as expressed in the following expression: φ_(j)(σ)=2πd _(j) B(σ)σ=2πL _(j)σΦ_(j)(σ)  (1.2) [Mathematical Expression 2] $\begin{matrix} {L_{j} = {{{\frac{1}{2\pi}\frac{\mathbb{d}\phi_{j}}{\mathbb{d}\sigma}}❘_{\sigma_{0}}} = {d_{j}\left( {{{B\left( \sigma_{0} \right)} + \frac{\mathbb{d}B}{\mathbb{d}\sigma}}❘_{\sigma_{0}}\sigma_{0}} \right)}}} & \left( {1.3a} \right) \\ {{{\Phi_{j}(\sigma)} = {{\left\{ {{\phi_{j}\left( \sigma_{0} \right)} - {2\pi\quad L_{j}\sigma_{0}}} \right\} + {\frac{1}{2}\frac{\mathbb{d}^{2}\phi_{j}}{\mathbb{d}\sigma^{2}}}}❘_{\sigma_{0}}{\left( {\sigma - \sigma_{0}} \right)^{2} + \ldots}}}\quad} & \left( {1.3b} \right) \end{matrix}$ where d_(j) is a thickness of R_(j), and B(σ) is its double refraction. Further, σ₀ indicates a center wavenumber of light under measurement. Hereinafter, the retardation φ_(j)(σ)) of the retarder is referred to as a reference phase function.

Assuming now that dispersion (change rate with wavenumber) of B(σ) is not so large, as seen from the expression (1.2), φ(σ) increases almost linearly with respect to wavenumber σ. Such a property serves as a basis of measurement of the spectropolarization parameter of the sample in a later-described procedure.

1.2 Channeled Spectrum Acquired in Spectrometer

In the “channeled spectroscopic polarimeter” (polarimetric spectroscope) shown in FIG. 2, a spectrum (spectral intensity) acquired in the spectrometer 8 is expressed by the following expression. [Mathematical Expression 3] $\begin{matrix} \begin{matrix} {{P(\sigma)} = {{\frac{1}{2}{m_{0}(\sigma)}{M_{0}(\sigma)}} + {\frac{1}{4}{m\_}(\sigma){{M_{23}(\sigma)}}\cos}}} \\ {\left\{ {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)} + {\arg\left( {M_{23}(\sigma)} \right)}} \right\} +} \\ {{\frac{1}{2}{m_{2}(\sigma)}{M_{1}(\sigma)}\cos\quad{\phi_{2}(\sigma)}} -} \\ {\frac{1}{4}{m_{+}(\sigma)}{{M_{23}(\sigma)}}\cos\left\{ {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)} - {\arg\left( {M_{23}(\sigma)} \right)}} \right\}} \end{matrix} & (1.4) \end{matrix}$ where M ₂₃(σ)=M ₂(σ)+iM ₃(σ)  (1.5) [Mathematical Expression 4] $\begin{matrix} {{M_{0}(\sigma)} = {\frac{1}{2}{{P_{0}(\sigma)}\left\lbrack {{{\hat{m}}_{00}(\sigma)} + {{{\hat{m}}_{10}(\sigma)}\cos\quad 2\theta} + {{{\hat{m}}_{20}(\sigma)}\sin\quad 2\theta}} \right\rbrack}}} & \left( {1.6a} \right) \\ {{M_{1}(\sigma)} = {\frac{1}{2}{{P_{0}(\sigma)}\left\lbrack {{{\hat{m}}_{01}(\sigma)} + {{{\hat{m}}_{11}(\sigma)}\cos\quad 2\theta} + {{{\hat{m}}_{21}(\sigma)}\sin\quad 2\theta}} \right\rbrack}}} & \left( {1.6b} \right) \\ {{M_{2}(\sigma)} = {\frac{1}{2}{{P_{0}(\sigma)}\left\lbrack {{{\hat{m}}_{02}(\sigma)} + {{{\hat{m}}_{12}(\sigma)}\cos\quad 2\theta} + {{{\hat{m}}_{22}(\sigma)}\sin\quad 2\theta}} \right\rbrack}}} & \left( {1.6c} \right) \\ {{M_{3}(\sigma)} = {\frac{1}{2}{{P_{0}(\sigma)}\left\lbrack {{{\hat{m}}_{03}(\sigma)} + {{{\hat{m}}_{13}(\sigma)}\cos\quad 2\theta} + {{{\hat{m}}_{23}(\sigma)}\sin\quad 2\theta}} \right\rbrack}}} & \left( {1.6d} \right) \end{matrix}$ Here, M₀(σ) to M₃(σ) are referred to as spectroscopic quasi-Stokes parameters of the sample. As thus described, the spectroscopic quasi-Stokes parameter is expressed by the sum of values each obtained by multiplying each element of each column of the Mueller matrix of the sample by a coefficient determined by an azimuth of the analyzer. The spectropolarization parameter of the sample can be obtained by simultaneously solving the equations of 1.6a to 1.6d. m₀(σ), m⁻(σ), m₂(σ), and m₊(σ) each shows a ratio of amplitude attenuation due to failure of the spectrometer to follow a fine vibration component. P₀(σ) shows a “spectrum of the light source”. However, in the optical system, attenuation exists caused by the retarder, the polarizer, a lens, a fiber or the like. Therefore, in the present specification, the “spectrum of the light source” P₀(σ) includes the attenuated portion. Further, φ₁ and φ₂ are retardations of the retarder R1 and R2.

The elements of the Mueller matrix of the sample included in M₀(σ) to M₃(σ) are each relative to “each column” of the Mueller matrix M(σ). [Mathematical Expression 5] $\begin{matrix} {{M(\sigma)} = \begin{bmatrix} \begin{matrix} {{\hat{m}}_{00}(\sigma)} \\ {{\hat{m}}_{10}(\sigma)} \\ {{\hat{m}}_{20}(\sigma)} \end{matrix} & \begin{matrix} {{\hat{m}}_{01}(\sigma)} \\ {{\hat{m}}_{11}(\sigma)} \\ {{\hat{m}}_{21}(\sigma)} \end{matrix} & \begin{matrix} {{\hat{m}}_{02}(\sigma)} \\ {{\hat{m}}_{12}(\sigma)} \\ {{\hat{m}}_{22}(\sigma)} \end{matrix} & \begin{matrix} {{\hat{m}}_{03}(\sigma)} \\ {{\hat{m}}_{13}(\sigma)} \\ {{\hat{m}}_{23}(\sigma)} \end{matrix} \\ {{\hat{m}}_{30}(\sigma)} & {{\hat{m}}_{31}(\sigma)} & {{\hat{m}}_{32}(\sigma)} & {{\hat{m}}_{33}(\sigma)} \end{bmatrix}} & (1.7) \end{matrix}$

Information inside the frames (sum of values each obtained by multiplying each element of each column by a coefficient determined by an azimuth angle (θ) of the analyzer A can be demodulated.

While 16 elements exist in the 4×4 Mueller matrix, it is extremely rare that all those elements are independent. In many cases, in polarimetry, only two to three independent parameters at most are included in the Mueller matrix of the sample. Even with a spectrum intensity of the light source included, the total number of parameters required to be measured is often as small as four at most. It is therefore possible by simultaneously solving obtained four equations to obtain up to four parameters, which are independent from each other and show the polarization characteristic of the sample.

For the sake of understanding the property of this expression, Expression (1.2) is substituted therein as follows. [Mathematical Expression 6] $\begin{matrix} \begin{matrix} {{P(\sigma)} = {{\frac{1}{2}{m_{0}(\sigma)}{M_{0}(\sigma)}} +}} \\ {{\frac{1}{4}{m\_}(\sigma){{M_{23}(\sigma)}}{\cos\left\lbrack {{2\pi\quad L_{-}\sigma} + {\Phi_{-}(\sigma)} + {\arg\left\{ {M_{23}(\sigma)} \right\}}} \right\rbrack}} +} \\ {{{\frac{1}{2}{m_{2}(\sigma)}{M_{1}(\sigma)}{\cos\left\lbrack {{2\pi\quad L_{2}\sigma} + {\Phi_{2}(\sigma)}} \right\rbrack}} -}\quad} \\ {\frac{1}{4}{m_{+}(\sigma)}{{M_{23}(\sigma)}}\cos} \\ {\left\lbrack {{2\pi\quad L_{+}\sigma} + {\Phi_{+}(\sigma)} - {\arg\left\{ {M_{23}(\sigma)} \right\}}} \right\rbrack} \end{matrix} & (1.8) \end{matrix}$ where it is found that the following expressions are satisfied. L ⁻ =L ₂ −L ₁  (1.9a) L ₊ =L ₂ +L ₁  (1.9b) Φ⁻(σ)=Φ₂(σ)−Φ₁ (σ)  (1.9c) Φ₊(σ)=Φ₂(σ)+Φ₁ (σ)  (1.9d)

As seen from Expression (1.8), the spectrum P(σ) obtained from the spectrometer contains four components. One of them is a component gently varies with wavenumber σ, and the other three components are quasi-sinusoidal components that vibrate with wavenumber σ. These are schematically shown in FIG. 4.

Here, the central periods of each of the three vibration components are respectively 1/L, 1/L₂ and 1/L₊. The spectrum containing components that periodically finely vibrate with wavenumber (wavelength) as in the figure is referred to as a channeled spectrum.

What needs to be concerned here is that these four components have information of any one of M₀(σ) M₁(σ) and M₂₃(σ). When each component can be separated, it is possible to determine all the spectroscopic quasi-Stokes parameters M₀(σ), M₁(σ), M₂(σ) and M₃ from one spectrum P(σ).

1.3 When Crossing Angle Between Elements is not 45°

Next described is a spectrum acquired in the spectrometer 5 when a crossing angle between the elements is not 45°.

Here also described as a supplemental explanation is a spectrum obtained when a crossing angle between each element in the optical system is not 45°.

It is assumed now that, in the optical system in FIG. 2, the angle formed between the fast axes of the retarders R1 and R2 is θ_(RR) and the angle formed between the fast axis of the retarder R2 and the transmission axis of the polarizer P is θ_(PR). So far, the calculation has been made only in a limited case of θ_(RR)=45° and θ_(PR)=45°. Below, a case where those angles are more common ones is shown.

An expression for the obtained channeled spectrum P(σ) is given as follows. [Mathematical Expression 7] $\begin{matrix} \begin{matrix} {{P(\sigma)} = {{\frac{1}{2}{{m_{0}(\sigma)}\left\lbrack {{M_{0}(\sigma)} + \underset{\_}{\cos\quad 2\quad\theta_{PR}\cos\quad 2\quad\theta_{RR}{M_{1}(\sigma)}}} \right\rbrack}} -}} \\ {\frac{1}{2}\left( {\sin\quad 2\quad\theta_{PR}\sin^{2}\theta_{RR}} \right){m_{-}(\sigma)}{{M_{23}(\sigma)}}\cos} \\ {\left\lbrack {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)} + {\arg\left\{ {M_{23}(\sigma)} \right\}}} \right\rbrack -} \\ {{\frac{1}{2}\left( {\sin\quad 2\quad\theta_{PR}\sin\quad 2\theta_{RR}} \right){m_{2}(\sigma)}{M_{1}(\sigma)}{\cos\left\lbrack {\phi_{2}(\sigma)} \right\rbrack}} +} \\ {\frac{1}{2}\left( {\sin\quad 2\quad\theta_{PR}\cos^{2}\theta_{RR}} \right){m_{+}(\sigma)}{{M_{23}(\sigma)}}\cos} \\ {\left\lbrack {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)} - {\arg\left\{ {M_{23}(\sigma)} \right\}}} \right\rbrack +} \\ {\underset{\_}{\frac{1}{2}\left( {\cos\quad 2\quad\theta_{PR}\sin\quad 2\theta_{RR}} \right){m_{1}(\sigma)}{{M_{23}(\sigma)}}\cos}} \\ \left. \underset{\_}{\left\lbrack {{\phi_{1}(\sigma)} - {\arg\left\{ {M_{23}(\sigma)} \right\}}} \right.} \right\rbrack \end{matrix} & (1.10) \end{matrix}$

When this expression is compared with the spectrum in the previous expression (1.4), namely when the angles θ_(RR) and θ_(PR) are respectively limited to 45° and 45°, the following differences are found in addition to a mere difference in constant multiple of a coefficient. It is to be noted that the different part is indicated by an underline in Expression (1.10).

The component that gently varies with wavenumber a depends not only upon M₀(σ) but additionally upon M₁(σ).

A component that quansi-sinusoidally vibrates according to the phase φ₁(σ), namely a component that vibrates at a central period 1/L₁, is added. It should be noted that this component has information of M₂₃(σ), as in the cases of the two respective components that vibrate according to (φ₂(σ)−φ₁(σ)) and (φ₂(σ)+φ₁(σ)). It means that this term can be treated in the same manner as the other two terms including M₂₃(σ).

Here, conditions for nonappearance of the above two components are considered.

The former term appears in a limited case “where both θ_(RR)≠±45° and θ_(PR)≠±45° are satisfied”. In the meantime, the latter term appears in a case “where θ_(PR)≠±45° “(regardless of whether θ_(RR) agrees with ±45° or not)”. From these, the following fact can be mentioned.

When the fast axis of the retarder R2 and the transmission axis of the polarizer P cross each other at an angle of 45° (i.e. θ_(PR)=±45°), the channeled spectrum is given by Expression (1.4) except for the difference in constant multiple of the coefficient of each term. Here, whether the angle θ_(RR) formed between the principle axes of the retarders R1 and R2 agrees with ±45° or not is irrelevant.

In other words, the channeled spectrum can take the form of Expression (1.4) under a condition that the fast axis of the retarder R2 and the transmission axis of the polarizer P cross each other at an angle of ±45°. Meanwhile, whether the angle formed between the fast axes of the retarders R1 and R2 agrees with ±45° or not is irrelevant.

1.4 When the Number of Retarders is Three or More

The case where the number of retarders is two was described above regarding the spectrum acquired in the spectrometer 5. In the same manner as this case, a spectrum having information of a spectroscopic quasi-Stokes parameter particular to each component is obtained in a case where the number of retarders is three or more. As in the case of the number of retarders is two, each element is separated so as to demodulate all spectroscopic quasi-Stokes parameters from one spectrum P(σ). By simultaneously solving the obtained expressions, it is possible to obtain a spectropolarization parameter of the sample.

1.5 Procedure for Demodulating Spectrometric Stokes Parameter

A specific procedure for demodulating a spectrometric quasi-Stokes parameter is described below with reference to FIG. 5. A brief description of the flow of the procedure is as follows.

Step 1: Each term is separated from Spectrum P(ca).

Step 2: An amplitude and a phase of each component are obtained. (Or equivalent quantities, e.g. a real part and an imaginary part in complex representation are obtained).

Step 3: Reference functions (below) included in the amplitude and phase of each vibration component are removed. [Mathematical Expression 8] ${Reference}\quad{amplitude}\quad{function}\begin{Bmatrix} {m_{0}(\sigma)} \\ {m_{-}(\sigma)} \\ {m_{2}(\sigma)} \\ {m_{+}(\sigma)} \end{Bmatrix}$ ${Reference}\quad{phase}\quad{function}\begin{Bmatrix} {\phi_{1}(\sigma)} \\ {\phi_{2}(\sigma)} \end{Bmatrix}$ Spectrometric Stokes parameters M₀(σ), M₁(σ), M₂(σ), and M₃(σ) are then obtained (These reference functions are peculiar to the polarimeter, depending not upon the sample but only upon parameters of the polarimeter.

Each of the steps is described as follows.

[Step 1]

As described in the previous section, the spectrum P(σ) contains four components. An operation for taking out each component by a signal processing is performed. What is applied to this operation is that each component vibrates at a different period (frequency). With the use of (any one of) a variety of frequency filtering techniques being broadly used in fields of communication engineering, signal analysis and the like, it is possible to separate each component. [Mathematical Expression 9] $\begin{matrix} {{{Component}\quad\lbrack 1\rbrack}\left( {{low}\quad{frequency}\quad{component}} \right)} & \quad \\ {\frac{1}{2}{m_{0}(\sigma)}{M_{0}(\sigma)}} & \left( {1.11a} \right) \\ {{{Component}\quad\lbrack 2\rbrack}\left( {{central}\quad{period}\quad{1/L_{-}}} \right)} & \quad \\ {\frac{1}{4}{m_{-}(\sigma)}{{M_{23}(\sigma)}}{\cos\left\lbrack {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)} + {\arg\left\{ {M_{23}(\sigma)} \right\}}} \right\rbrack}} & \left( {1.11b} \right) \\ {{{Component}\quad\lbrack 3\rbrack}\left( {{central}\quad{period}\quad{1/L_{2}}} \right)} & \quad \\ {\frac{1}{2}{m_{2}(\sigma)}{M_{1}(\sigma)}{\cos\left\lbrack {\phi_{2}(\sigma)} \right\rbrack}} & \left( {1.11c} \right) \\ {{{Component}\quad\lbrack 4\rbrack}\left( {{central}\quad{period}\quad{1/L_{+}}} \right)} & \quad \\ {{- \frac{1}{4}}{m_{+}(\sigma)}{{M_{23}(\sigma)}}{\cos\left\lbrack {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)} - {\arg\left\{ {M_{23}(\sigma)} \right\}}} \right\rbrack}} & \left( {1.11d} \right) \end{matrix}$

Component [1] above is a first spectral intensity component which nonperiodically vibrates with wavenumber. Component [2] is a second spectral intensity component which vibrates at a frequency depending upon a difference between a reference phase function (retardation) φ₁(σ) of the first retarder R1 and a reference phase function (retardation) φ₂(σ) of the second retarder R2 with wavenumber. Component [3] is a third spectral intensity component which vibrates at a frequency depending upon the reference phase function φ₂(σ) of the second retarder and not depending upon the reference phase function φ₁(σ) of the first retarder, with wavenumber. Component [4] is a fourth spectral intensity component which vibrates at a frequency depending upon the sum of the reference phase function φ₁(σ) of the first retarder and the reference phase function φ₂(σ) of the second retarder with wavenumber. It is to be noted that, when the crossing angle between the elements is not 45°, Component [5] is generated. Component [5] vibrates at a frequency depending upon the reference phase function φ₁(σ) of the first retarder and not depending upon the reference phase function φ₂(σ) of the second retarder, with wavenumber.

With reference to Expressions (1.11a) to (1.11d) and Expression (1.5), it is found that M₀ is obtained from Component [1], M₁ is from Component [3], M₂ and M₃ are from Component [2] or Component [4]. It should be noted that, when the crossing angle between the elements is not 45°, M₂ and M₃ are obtained at least one of Component [2], Component [4] and Component [5].

[Step 2]

As for each component separated in Step 1, a “paired amplitude and phase” and a “complex representation” are obtained, as shown in FIG. 6. This can be readily realized by using a variety of demodulation methods in an operation which are common in fields of communication engineering, signal analysis and the like, as in

Step 1. Examples of those methods include:

Amplitude demodulation: rectifying demodulation, envelope demodulation, etc.

Phase demodulation: frequency discrimination, zero-crossing method, etc.

Complex representation demodulation: Fourier transform method (later described), synchronous demodulation, etc.

Here, definitions and basic properties of the “amplitude”, “phase” and “complex representation” of a vibration component are summarized below. As seen from Expressions (1.11 a) to (1.11 d), each of the separated components except for Component [1] takes the form of: a(σ)cos δ(σ). a(σ) and δ(σ) here are respectively referred to as the “amplitude” and “phase” of the vibration component. It is to be noted that, if assuming that the phase Δ₀(σ)=0 (i.e. cos δ₀(σ)=0) also in Component [1], the amplitude of Component [1] can also be defined.

Further, F(σ) having the following relation with the amplitude and the phase is called a complex representation. [Mathematical Expression 10] $\begin{matrix} {{F(\sigma)} = {\frac{1}{2}{a(\sigma)}{\exp\left\lbrack {i\quad{\delta(\sigma)}} \right\rbrack}}} & \left( {1.12a} \right) \\ {\quad{= {\left\lbrack {\frac{1}{2}{a(\sigma)}\cos\quad{\delta(\sigma)}} \right\rbrack + {i\left\lbrack {\frac{1}{2}{a(\sigma)}\sin\quad{\delta(\sigma)}} \right\rbrack}}}} & \left( {1.12b} \right) \end{matrix}$ The real part of F(σ) is formed by dividing the amplitude of the vibration component into halves, and the imaginary part thereof is displaced from the real part at the angle of 90°. It should be noted that in Component [1], the amplitude is not divided into halves since δ(σ)=0, i.e. no imaginary part exists.

What needs to be concerned here is that when either the “paired amplitude and phase” or the “complex representation” can be demodulated, the other one can be immediately calculated by use of the following relational expression. [Mathematical Expression 8] $\begin{matrix} {\left. {``{{{amplitude}\quad a\quad(\sigma)},{{phase}\quad\delta\quad(\sigma)}}"}\rightarrow{``{{complex}\quad{representation}\quad F\quad(\sigma)}"} \right.{{F(\sigma)} = {\frac{1}{2}{a(\sigma)}{\mathbb{e}}^{i\quad{\delta{(\sigma)}}}}}} & (1.13) \\ {\left. {``{{complex}\quad{representation}\quad F\quad(\sigma)}"}\rightarrow{``{{{amplitude}\quad a\quad(\sigma)},{{phase}\quad\delta\quad(\sigma)}}"} \right.{{a(\sigma)} = {2{{F(\sigma)}}}}} & \left( {1.14a} \right) \\ {{\delta(\sigma)} = {\arg\left\lbrack {F(\sigma)} \right\rbrack}} & \left( {1.14b} \right) \end{matrix}$

Namely, demodulation of one of the paired amplitude and phase and the complex representation enables immediate calculation of the other according to need.

When the “amplitude” and “phase” of each component are demodulated, the following results are obtained. [Mathematical Expression 12] $\begin{matrix} {{\bullet\quad{{Component}\quad\lbrack 1\rbrack}\quad\left( {{low}\quad{frequency}\quad{component}} \right)\text{:}}\quad\begin{matrix} \lbrack{Amplitude}\rbrack & {{a_{0}(\sigma)} = {\frac{1}{2}{m_{0}(\sigma)}{M_{0}(\sigma)}}} \\ \lbrack{Phase}\rbrack & {{\delta_{0}(\sigma)} = 0} \end{matrix}} & \left( {1.15a} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 2\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{-}}} \right)\text{:}}\quad\begin{matrix} \lbrack{Amplitude}\rbrack & {{a_{-}(\sigma)} = {\frac{1}{4}{m_{-}(\sigma)}{{M_{23}(\sigma)}}}} \\ \lbrack{Phase}\rbrack & {{\delta_{-}(\sigma)} = {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)} + {\arg\left\{ {M_{23}(\sigma)} \right\}}}} \end{matrix}} & \left( {1.15b} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 3\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{2}}} \right)\text{:}}\quad\begin{matrix} \lbrack{Amplitude}\rbrack & {{a_{2}(\sigma)} = {\frac{1}{2}{m_{2}(\sigma)}{M_{1}(\sigma)}}} \\ \lbrack{Phase}\rbrack & {{\delta_{2}(\sigma)} = {\phi_{2}(\sigma)}} \end{matrix}} & \left( {1.15c} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 4\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{+}}} \right)\text{:}}\quad\begin{matrix} {\lbrack{Amplitude}\rbrack} & {{a_{+}(\sigma)} = {\frac{1}{4}{m_{+}(\sigma)}{{M_{23}(\sigma)}}}} \\ \lbrack{Phase}\rbrack & {{\delta_{+}(\sigma)} = {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)} - {\arg\left\{ {M_{23}(\sigma)} \right\}} + \pi}} \end{matrix}} & \left( {1.15d} \right) \end{matrix}$

On the other hand, when the “complex representation” of each component is demodulated, the following results are obtained. [Mathematical Expression 13] $\begin{matrix} {{\bullet\quad{{Component}\quad\lbrack 1\rbrack}\quad\left( {{low}\quad{frequency}\quad{component}} \right){\text{:}\left\lbrack {{Complex}\quad{representation}} \right\rbrack}\quad{F_{0}(\sigma)}} = {\frac{1}{2}{m_{0}(\sigma)}{M_{0}(\sigma)}}} & \left( {1.16a} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 2\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{-}}} \right){\text{:}\left\lbrack {{Complex}\quad{representation}} \right\rbrack}\quad{F_{-}(\sigma)}} = {\frac{1}{8}{m_{-}(\sigma)}{M_{23}(\sigma)}\exp\quad{{\mathbb{i}}\left\lbrack {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)}} \right\rbrack}}} & \left( {1.16b} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 3\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{2}}} \right){\text{:}\left\lbrack {{Complex}\quad{representation}} \right\rbrack}\quad{F_{2}(\sigma)}} = {\frac{1}{4}{m_{2}(\sigma)}{M_{1}(\sigma)}\exp\quad{\mathbb{i}}\quad{\phi_{2}(\sigma)}}} & \left( {1.16c} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 4\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{+}}} \right){\text{:}\left\lbrack {{Complex}\quad{representation}} \right\rbrack}\quad{F_{+}(\sigma)}} = {{- \frac{1}{8}}{m_{+}(\sigma)}{M_{23}^{*}(\sigma)}\exp\quad{{\mathbb{i}}\left\lbrack {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)}} \right\rbrack}}} & \left( {1.16c} \right) \end{matrix}$

Here, * indicates a complex conjugation. It is to be noted that, for the sake of what is described below, the expressions of the complex representations are rewritten as follows. [Mathematical Expression 14] $\begin{matrix} {{\bullet\quad{{Component}\quad\lbrack 1\rbrack}\quad\left( {{low}\quad{frequency}\quad{component}} \right){\text{:}\quad\left\lbrack {{Complex}\quad{representation}} \right\rbrack}\quad{F_{0}(\sigma)}} = {{K_{0}(\sigma)}{M_{0}(\sigma)}}} & \left( {1.17a} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 2\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{-}}} \right){\text{:}\quad\left\lbrack {{Complex}\quad{representation}} \right\rbrack}\quad{F_{-}(\sigma)}} = {{K_{-}(\sigma)}{M_{23}(\sigma)}}} & \left( {1.17b} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 3\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{2}}} \right){\text{:}\quad\left\lbrack {{Complex}\quad{representation}} \right\rbrack}\quad{F_{2}(\sigma)}} = {{K_{2}(\sigma)}{M_{1}(\sigma)}}} & \left( {1.17c} \right) \\ {{{\bullet\quad{{Component}\quad\lbrack 4\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{+}}} \right){\text{:}\quad\left\lbrack {{Complex}\quad{representation}} \right\rbrack}\quad{F_{+}(\sigma)}} = {{K_{+}(\sigma)}{M_{23}^{*}(\sigma)}}}{where}} & \left( {1.17d} \right) \\ {{K_{0}(\sigma)} = {\frac{1}{2}{m_{0}(\sigma)}}} & \left( {1.18a} \right) \\ {{K_{-}(\sigma)} = {\frac{1}{8}{m_{-}(\sigma)}\exp\quad{{\mathbb{i}}\left\lbrack {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)}} \right\rbrack}}} & \left( {1.18b} \right) \\ {{K_{2}(\sigma)} = {\frac{1}{4}{m_{2}(\sigma)}\exp\quad{\mathbb{i}}\quad{\phi_{2}(\sigma)}}} & \left( {1.18c} \right) \\ {{K_{+}(\sigma)} = {{- \frac{1}{8}}{m_{+}(\sigma)}\exp\quad{{\mathbb{i}}\left\lbrack {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)}} \right\rbrack}}} & \left( {1.18d} \right) \end{matrix}$ [Step 3]

Finally, from the “amplitude” and the “phase” or the “complex representation” obtained in Step 2 above, the spectrometer parameters M₀(σ), M₁(σ), M₂(σ), and M₃(σ) are determined as functions of the wavenumber σ.

The “amplitude” and the “phase” obtained in Step 2 above include, other than the spectrometric Stokes parameters to be obtained, parameters shown below. [Mathematical Expression 15] Parameter  (function)  determined  based  only  upon  characteristic of  polarimeter  itself $\begin{Bmatrix} {m_{0}(\sigma)} \\ {m_{-}(\sigma)} \\ {m_{2}(\sigma)} \\ {m_{+}(\sigma)} \end{Bmatrix}$ and $\begin{Bmatrix} {{\phi_{-}(\sigma)} = {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)}}} \\ {\phi_{2}(\sigma)} \\ {{\phi_{+}(\sigma)} = {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)}}} \end{Bmatrix}$

The former are included in the amplitude while the latter are included in the phase. These parameters provide references in determining spectrometric quasi-Stokes parameters from the amplitude and the phase of each vibration component. Thus, each of these parameters is hereinafter referred to as a “reference amplitude function” or a “reference phase function”. Since these parameters do not depend upon the sample, each of the parameters is subjected to division or subtraction, to be determined as follows.

-   -   M₀(σ) can be determined from [Component [1]].     -   M₂(σ) and M₃(σ) can be determined from (either) [Component [2]]         or [Component [4]].     -   M₁(σ) can be determined from [Component [3]].

Meanwhile, in the case of the “complex representation”, parameters (functions) determined only by the characteristic of the polarimeter itself are K₀(σ), K⁻(σ), K₂(σ), and K₊(σ) which are defined by Expressions (1.18a) to (1.18d). These are, so to speak, “reference complex functions”.

As revealed from Expressions (1.17a) to (1.17d), if the above reference complex functions have been obtained, by division of the complex representation of each vibration component demodulated in Step 2, the parameters can be determined as follows.

-   -   M₀(σ) can be determined from [Component [1]]     -   M₂(σ) and M₃(σ) can be determined from (either) [Component [2]]         or [Component [4]].     -   M₁(σ) can be determined from [Component [3]].

When the angle formed between the retarder R2 and the polarizer P is not 45°, a fifth term that appears can be used in place of [component [2]] and [component [4]]. Namely, the description on lines 2-3 above can be rewritten to:

-   -   M₂(σ) and M₃(σ) can be determined from any one of [component         [2]], [component [4]] and [component [5]].

Next, as one of signal processing methods for demodulating spectrometric quasi-Stokes parameters, a “Fourier transform method” is described with reference to FIG. 7. The use of this method allows efficient concurrent performance of Steps 1 and 2, leading to immediate determination of all complex representations of each vibration component.

In this method, first, the spectrum P(σ) measured with the spectrometer in the channeled spectroscopic polarimeter is subjected to inverse Fourier transformation, to obtain the following correlation function of light incident on the spectrometer. [Mathematical Expression 16] $\begin{matrix} {{{C(h)} = {{A_{0}(h)} + {A_{-}\left( {h - L_{-}} \right)} + {A_{-}^{*}\left( {{- h} - L_{-}} \right)} + {A_{2}\left( {h - L_{2}} \right)} + {A_{2}^{*}\left( {{- h} - L_{2}} \right)} + {A_{+}\left( {h - L_{+}} \right)} + {A_{+}^{*}\left( {{- h} - L_{+}} \right)}}}{where}} & (1.19) \\ {{A_{0}(h)} = {F^{- 1}\left\lbrack {\frac{1}{2}{m_{0}(\sigma)}{M_{0}(\sigma)}} \right\rbrack}} & \left( {1.20a} \right) \\ {{A_{-}(h)} = {F^{- 1}\left\lbrack {\frac{1}{8}{m_{-}(\sigma)}{M_{23}(\sigma)}\exp\quad{\mathbb{i}}\quad{\Phi_{-}(\sigma)}} \right\rbrack}} & \left( {1.20b} \right) \\ {{A_{2}(h)} = {F^{- 1}\left\lbrack {\frac{1}{4}{m_{2}(\sigma)}{M_{1}(\sigma)}\exp\quad{\mathbb{i}}\quad{\Phi_{2}(\sigma)}} \right\rbrack}} & \left( {1.20c} \right) \\ {{A_{+}(h)} = {F^{- 1}\left\lbrack {{- \frac{1}{8}}{m_{+}(\sigma)}{M_{23}^{*}(\sigma)}\exp\quad{\mathbb{i}}\quad{\Phi_{+}(\sigma)}} \right\rbrack}} & \left( {1.20d} \right) \end{matrix}$ As shown in the upper right part of FIG. 7, this correlation function C(h) contains seven components including inverse numbers of a period of each vibration component, 0, ±L⁻, ±L₂, ±L₊ as main components.

Here, appropriate selection of these inverse numbers of the period enables separation of each component, contained in C(h), from each other on the h-axis. When four components as the main components, h=0, L⁻, L₂, and L₊, are taken out and then subjected to the Fourier transformation, the following expressions are satisfied. [Mathematical Expression 17] $\begin{matrix} {{F\left\lbrack {A_{0}(h)} \right\rbrack} = {{\frac{1}{2}{m_{0}(\sigma)}{M_{0}(\sigma)}} = {F_{0}(\sigma)}}} & \left( {1.21a} \right) \\ \begin{matrix} {{F\left\lbrack {A_{-}\left( {h - L_{-}} \right)} \right\rbrack} = {\frac{1}{8}{m_{-}(\sigma)}{M_{23}(\sigma)}\exp\quad{{\mathbb{i}}\left\lbrack {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)}} \right\rbrack}}} \\ {= {F_{-}(\sigma)}} \end{matrix} & \left( {1.21b} \right) \\ \begin{matrix} {{F\left\lbrack {A_{2}\left( {h - L_{2}} \right)} \right\rbrack} = {\frac{1}{4}{m_{2}(\sigma)}{M_{1}(\sigma)}\exp\quad{\mathbb{i}}\quad{\phi_{2}(\sigma)}}} \\ {= {F_{2}(\sigma)}} \end{matrix} & \left( {1.21c} \right) \\ \begin{matrix} {{F\left\lbrack {A_{+}\left( {h - L_{+}} \right)} \right\rbrack} = {{- \frac{1}{8}}{m_{+}(\sigma)}{M_{23}^{*}(\sigma)}\exp\quad{{\mathbb{i}}\left\lbrack {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)}} \right\rbrack}}} \\ {= {F_{+}(\sigma)}} \end{matrix} & \left( {1.21d} \right) \end{matrix}$

As seen from the expressions above, what are obtained in the above operation are just the complex representations of the components [1] to [4] to be obtained in foregoing Step 2. Namely, in the above operation, Steps 1 and 2 are concurrently realized. Hence, when Step 3 is performed using the results of Steps 1 and 2, spectrometric quasi-Stokes parameters are obtained all at once.

1.6 Pre-Calibration: Calibration of Reference Amplitude Function, Reference Phase Function, Reference Complex Function “Prior to Measurement”

As described in the previous section, when a spectrometric quasi-Stokes parameter is determined from a channeled spectrum, it is necessary in Step 3 to in advance determine parameters to be obtained based only upon a characteristic of a polarimeter itself parameters, namely:

“reference amplitude function” m₀(σ), m⁻(σ), m₂(σ), m₊(σ), and “reference phase function” φ₂(σ) and φ₁(σ), or

“reference complex function” K₀(σ), K⁻(σ), K₂(σ), K₊(σ)

The former (“reference amplitude function” and “reference phase function”) and the latter (“reference complex function”) are required in the respective cases of obtaining spectrometric quasi-Stokes parameters from the “amplifier and phase” and the “complex representation” of each vibration component. Since these are functions not depending upon the sample, it is desirable to calibrate the functions at least prior to measurement.

In this section, a procedure for calibrating these reference functions “prior to measurement, i.e. in advance” is described. There are two typical methods as follows.

-   -   [Method 1]: a method for calibrating reference phase functions         and reference amplitude functions from a characteristic of each         element for use in the optical system.     -   [Method 2]: a method for calibrating reference phase functions         and reference amplitude functions by use of a sample having a         known polarization characteristic.

1.6.1 [Method 1]

Method for Calibrating Reference Phase Function and Reference Amplitude Function from Characteristic of Each Element for Use in Optical System

Characteristics of a reference phase function and a reference amplitude function are essentially determined based upon elements for use in a channeled spectroscopic polarimeter. Therefore, optical characteristics of individual elements are repeatedly examined by experiment or calculation to perform calibration of parameters.

1.6.2 [Method 2]

Method for Calibrating Reference Phase Function and Reference Amplitude Function by Use of Sample Having Known Polarization Characteristic.

The reference phase function and the reference amplitude function are in amount determined based not upon the “polarization characteristic of the object to be measured (sample)”, but only upon the characteristic of the “channeled spectroscopic polarimeter”. Accordingly, a “sample having a known polarization characteristic” (light whose measurement result is known)” is inputted into the polarimeter, and using the result of the input, it is possible to calculate backward the reference phase function and the reference amplitude function.

Below, the procedure for calibration is shown. As described at the beginning of this section, the following should be noted.

-   -   When the SOP is obtained from the “amplitude and phase” of each         vibration component, the “reference amplitude function” and the         “reference phase function” are required.     -   When the SOP is obtained from the “complex representation” of         each vibration component, the “reference complex function” is         required.

In the following, the procedures for calibration in the above respective cases are described. Although these procedures are essentially equivalent and different only in calculation method, they are separately put down for the sake of convenience.

A. Calibration Procedure for Separately Obtaining Reference Amplitude Function and Reference Phase Function

In this calibration, first, a “sample having a known polarization characteristic” is prepared, and then inputted into a channeled spectroscopic polarimeter. In this case, spectrometric quasi-Stokes parameters of light are referred to as M₀ ⁽⁰⁾(σ) M₁ ⁽⁰⁾(σ), M₂ ⁽⁰⁾(σ), and M₃ ⁽⁰⁾(σ). When the sampler is subjected to the above-mentioned demodulation means, the amplitude and the phase obtained in Step 2 are expressed as follows according to Expressions (1.15a) to (1.15d). [Mathematical Expression 18] $\begin{matrix} {{\bullet\quad{{Component}\quad\lbrack 1\rbrack}\quad\left( {{low}\quad{frequency}\quad{component}} \right)\text{:}}\quad\begin{matrix} \lbrack{Amplitude}\rbrack & {{a_{0}^{(0)}(\sigma)} = {\frac{1}{2}{m_{0}(\sigma)}{M_{0}^{(0)}(\sigma)}}} \\ \lbrack{Phase}\rbrack & {{\delta_{0}^{(0)}(\sigma)} = 0} \end{matrix}} & \left( {1.22a} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 2\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{-}}} \right)\text{:}}\quad\begin{matrix} \lbrack{Amplitude}\rbrack & {{a_{-}^{(0)}(\sigma)}\quad = \quad{\frac{1}{4}\quad{m_{-}(\sigma)}\quad{{M_{23}^{(0)}(\sigma)}}}} \\ \lbrack{Phase}\rbrack & {{\delta_{-}^{(0)}(\sigma)}\quad = \quad{{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)} + {\arg\left\{ {M_{23}^{(0)}(\sigma)} \right\}}}} \end{matrix}} & \left( {1.22b} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 3\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{2}}} \right)\text{:}}\quad\begin{matrix} \lbrack{Amplitude}\rbrack & {{a_{2}^{(0)}(\sigma)}\quad = \quad{\frac{1}{2}\quad{m_{2}(\sigma)}\quad{M_{1}^{(0)}(\sigma)}}} \\ \lbrack{Phase}\rbrack & {{\delta_{2}^{(0)}(\sigma)}\quad = \quad{\phi_{2}(\sigma)}} \end{matrix}} & \left( {1.22c} \right) \\ {{\bullet\quad{{Component}\quad\lbrack 4\rbrack}\quad\left( {{central}\quad{period}\quad{1/L_{+}}} \right)\text{:}}\quad\begin{matrix} \lbrack{Amplitude}\rbrack & {{a_{+}^{(0)}(\sigma)}\quad = \quad{\frac{1}{4}\quad{m_{+}(\sigma)}\quad{{M_{23}^{(0)}(\sigma)}}}} \\ \lbrack{Phase}\rbrack & {{\delta_{+}^{(0)}(\sigma)}\quad = \quad{{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)} - {\arg\left\{ {M_{23}^{(0)}(\sigma)} \right\}} + \pi}} \end{matrix}{where}} & \left( {1.22d} \right) \\ {{M_{23}^{(0)}(\sigma)} = {{M_{2}^{(0)}(\sigma)} - {{\mathbb{i}}\quad{M_{3}^{(0)}(\sigma)}}}} & (1.23) \end{matrix}$ It is to be noted that this is mere replacement of M₀(σ) to M₃(σ) with M₀ ⁽⁰⁾(σ) to M₃ ⁽⁰⁾(σ).

The amplitude and the phase of each vibration component are determined only by the spectrometric quasi-Stokes parameters, the reference phase functions and the reference amplitude functions. Here, since the spectrometric quasi-Stokes parameters are known in a “case where the sample having a known polarization characteristic”, the remaining reference amplitude functions m₀(σ), m⁻(σ), m₂(σ), m₊(σ), and the remaining reference phase functions φ₁(σ) and φ₂(σ) are determined from the demodulated amplitude and phase. Specifically, these functions are given according to the following expressions: [Mathematical Expression 19] $\begin{matrix} {{m_{0}(\sigma)} = \frac{2{a_{0}^{(0)}(\sigma)}}{M_{0}^{(0)}(\sigma)}} & \left( {1.24a} \right) \\ {{m_{-}(\sigma)} = \frac{4{a_{-}^{(0)}(\sigma)}}{{M_{23}^{(0)}(\sigma)}}} & \left( {1.24b} \right) \\ {{m_{2}(\sigma)} = \frac{2{a_{2}^{(0)}(\sigma)}}{M_{1}^{(0)}(\sigma)}} & \left( {1.24c} \right) \\ {{m_{+}(\sigma)} = \frac{4{a_{+}^{(0)}(\sigma)}}{{M_{23}^{(0)}(\sigma)}}} & \left( {1.24d} \right) \\ \begin{matrix} {{\phi_{-}\quad(\sigma)} = {{\phi_{\quad 2}\quad(\sigma)}\quad - \quad{\phi_{\quad 1}\quad(\sigma)}}} \\ {= {{\delta_{-}^{(0)}\quad(\sigma)} - {\arg\left\{ {M_{\quad 23}^{(0)}\quad(\sigma)} \right\}}}} \end{matrix} & \left( {1.24e} \right) \\ {{\phi_{\quad 2}(\sigma)} = {\delta_{\quad 2}^{(0)}\quad(\sigma)}} & \left( {1.24f} \right) \\ \begin{matrix} {{\phi_{+}(\sigma)} = {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)}}} \\ {= {{\delta_{+}^{90)}(\sigma)} + {\arg\left\{ {M_{23}^{(0)}(\sigma)} \right\}} - \pi}} \end{matrix} & \left( {1.24g} \right) \end{matrix}$ Once these reference functions are obtained (can be calibrated), spectrometric quasi-Stokes parameters of a sample having an unknown spectropolarization characteristic can be determined.

When the case of leaving the analyzer without using the sample is considered as an example, the following expressions are formed where the azimuth angle of the analyzer with respect to the fast axis of the retarder R1 is θ. M ₀ ⁽⁰⁾(σ)=P ₀ ⁽⁰⁾(σ)/2  (1.25a) M ₁ ⁽⁰⁾(σ)=P ₀ ⁽⁰⁾(σ)cos 2θ/2  (1.25b) M ₂ ⁽⁰⁾(σ)=P ₀ ⁽⁰⁾(σ)sin 2θ/2  (1.25c) M ₃ ⁽⁰⁾(σ)=0  (1.25d) Here, P₀ ⁽⁰⁾(σ) is a spectrum of the light source. In this case, the above expressions (1.24a) to (1.24g) are expressed as follows. [Mathematical Expression 20] $\begin{matrix} {{m_{0}(\sigma)} = \frac{2{a_{0}^{(0)}(\sigma)}}{\frac{1}{2}{P_{0}^{(0)}(\sigma)}}} & \left( {1.26a} \right) \\ {{m_{-}(\sigma)} = \frac{4{a_{-}^{(0)}(\sigma)}}{\frac{1}{2}{P_{0}^{(0)}(\sigma)}\sin\quad 2\quad\theta}} & \left( {1.26b} \right) \\ {{m_{2}(\sigma)} = \frac{2{a_{2}^{(0)}(\sigma)}}{\frac{1}{2}{P_{0}^{(0)}(\sigma)}\cos\quad 2\quad\theta}} & \left( {1.26c} \right) \\ {{m_{+}(\sigma)} = \frac{4{a_{+}^{(0)}(\sigma)}}{\frac{1}{2}{P_{0}^{(0)}(\sigma)}\sin\quad 2\quad\theta}} & \left( {1.26d} \right) \\ {{\phi_{-}(\sigma)} = {{{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)}} = {\delta_{-}^{(0)}(\sigma)}}} & \left( {1.26e} \right) \\ {{\phi_{2}(\sigma)} = {\delta_{2}^{(0)}(\sigma)}} & \left( {1.26f} \right) \\ {{\phi_{+}(\sigma)} = {{{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)}} = {{\delta_{+}^{(0)}(\sigma)} - \pi}}} & \left( {1.26g} \right) \end{matrix}$

It is revealed from the above that the reference amplitude function and the reference phase function can be obtained so long as the direction θ and the spectrum P₀ ⁽⁰⁾(σ) of the light source are known in advance. Further, even with P₀ ⁽⁰⁾(σ) unknown, if only the direction θ is known, it serves sufficiently for obtaining part of (essential) polarization parameters.

B. Calibration Procedure for Obtaining Amplitude and Phase Altogether (by Regarding Both as Reference Complex Function) at Once

The above-mentioned method was a method for calculating the “amplitude” and the “phase” of each vibration component separately. However, it may be more convenient (efficient) in some cases to calculate them as the “complex representation” of each vibration component. One example of such calculation may be the case of directly obtaining the “complex representation” (Expressions (1.17a) to (1.17d)), as in the Fourier transform method shown in FIG. 7 above. In such a case, calibration is efficiently performed by calibration of the “complex representation” as it is without separation into the “amplitude” and “phase”.

In the following, mathematical expressions for the above-mentioned case are shown. What needs to be concerned here is that the physical properties of the cases of using “amplitude and phase” and the “complex representation” are completely the same. It is just that in the latter case, a calculation is made using complex numbers, and thus more efficient.

Similarly to the previous section, a case is considered where the sample having a known polarization characteristic is inserted into the channeled spectroscopic polarimeter. A complex representation of each vibration component is obtained according to Expressions (1.17a) to (1.17d) as follows. F ₀ ⁽⁰⁾(σ)=K ₀(σ)M ₀ ⁽⁰⁾(σ)  (1.27a) F ⁻ ⁽⁰⁾(σ)=K ⁻(σ)M ₂₃ ⁽⁰⁾(σ)  (1.27b) F ₂ ⁽⁰⁾(σ)=K ₂(σ)M ₁ ⁽⁰⁾(σ)  (1.27c) F ₊ ⁽⁰⁾(σ)=K ₊(σ)M ₂₃ ⁽⁰⁾⁺(σ)  (1.27d)

Here, the complex functions K₀(σ), K⁻(σ), K₂(σ), and K₊(σ) are in amount (reference complex function) determined based not upon the sample, but only upon the reference amplitude function and the reference phase function, as seen from Expressions (1.18a) to (1.18d). Accordingly, these can be calculated backward as follows. [Mathematical Expression 21] $\begin{matrix} {{K_{0}(\sigma)} = \frac{F_{0}^{(0)}(\sigma)}{M_{0}^{(0)}(\sigma)}} & \left( {1.28a} \right) \\ {{K_{-}(\sigma)} = \frac{F_{-}^{(0)}(\sigma)}{M_{23}^{(0)}(\sigma)}} & \left( {1.28b} \right) \\ {{K_{2}(\sigma)} = \frac{F_{2}^{(0)}(\sigma)}{M_{1}^{(0)}(\sigma)}} & \left( {1.28c} \right) \\ {{K_{+}(\sigma)} = \frac{F_{+}^{(0)}(\sigma)}{M_{23}^{(0)}*(\sigma)}} & \left( {1.28d} \right) \end{matrix}$

Similar to the case of calculating the amplitude and the phase separately, once the above reference complex function is obtained (can be calibrated), spectrometric quasi-Stokes parameters of the sample having an unknown polarization characteristic can be obtained.

It is to be noted that, just for reference, mathematical expressions in the case of leaving the analyzer without using the sample are shown below. [Mathematical Expression 22] $\begin{matrix} {{K_{0}(\sigma)} = \frac{F_{0}^{(0)}(\sigma)}{\frac{1}{2}{P_{0}^{(0)}(\sigma)}}} & \left( {1.29a} \right) \\ {{K_{-}(\sigma)} = \frac{F_{-}^{(0)}(\sigma)}{\frac{1}{2}{P_{0}^{(0)}(\sigma)}\sin\quad 2\theta}} & \left( {1.29b} \right) \\ {{K_{2}(\sigma)} = \frac{F_{2}^{(0)}(\sigma)}{\frac{1}{2}{P_{0}^{(0)}(\sigma)}\cos\quad 2\theta}} & \left( {1.29c} \right) \\ {{K_{+}(\sigma)} = \frac{F_{+}^{(0)}(\sigma)}{\frac{1}{2}{P_{0}^{(0)}(\sigma)}\sin\quad 2\theta}} & \left( {1.29d} \right) \end{matrix}$

Chapter 2: Problems of Channeled Spectroscopic Polarimeter

As described in Step 3 in Section 1.5, for demodulation of spectrometric quasi-Stokes parameters M₀(σ), M₁ (σ), M₂(σ), and M₃(σ) from a measured channeled spectrum P(σ), it is necessary to obtain (calibrate) the following functions in advance. [Mathematical Expression 23] $\left. {\left. {\left. {\begin{matrix} {Reference} \\ {amplitude} \\ {{function},} \end{matrix}\begin{matrix} {m_{0}(\sigma)} \\ {m_{-}(\sigma)} \\ {m_{2}(\sigma)} \\ {m_{+}(\sigma)} \end{matrix}} \right\}\quad\begin{matrix} {Reference} \\ {phase} \\ {{function},{or}} \end{matrix}\begin{matrix} {\phi_{1}(\sigma)} \\ {\phi_{2}(\sigma)} \end{matrix}} \right\}\quad\begin{matrix} {Reference} \\ {complex} \\ {function} \end{matrix}\begin{matrix} {K_{0}(\sigma)} \\ {K_{-}(\sigma)} \\ {K_{2}(\sigma)} \\ {K_{+}(\sigma)} \end{matrix}} \right\}$

However, the reference phase functions φ₁(σ) and φ₂(σ) have the property of varying for a variety of reasons.

When these functions vary, there occurs a problem in that a large error occurs in measured values of the spectropolarization parameters of the sample.

2.1 Cause of Variations in Reference Phase Function 2.1.1 Temperature Change

The reference phase functions φ₁(σ) and φ₂(σ) are amounts (retardation) determined by the retarders R1 and R2 in the spectroscopic polarimeter. This retardation has the property of changing sensitively with respect to a temperature. Hence the phase of the channeled spectrum is displaced due to the temperature change. As a result, due to a temperature rise, a measured value is displaced to cause occurrence of an error therein. Moreover, a similar change occurs with pressure change.

2.1.2 Variations in Wavelength Axis of Spectrometer

When a wavelength to be sampled with the spectrometer is displaced, a problem that is “equivalent” to fluctuations in the reference phase function occurs. When the wavelength to be sampled is displaced during measurement, a similar effect to an effect in lateral displacement of the spectrum is produced. This is an equivalent phase displacement. In particular, in an ordinary spectrometer (type of rotating a diffraction grating with a motor), a wavelength to be sampled is displaced by small degrees (at random) in every measurement, due to backlash of a motor or the like.

2.1.3 Solution Easily Found

For preventing variations in reference phase function of each vibration component, stabilizing a cause of the fluctuations is considered. However, this is very hard to realize. For example, when noting the temperature variation, the accuracy required for wavenumber-distribution of an ellipsometric parameter in spectrometric ellipsometry is in the order of 0.1° or smaller, and for satisfying this, it is necessary to keep the temperature variation within the order of 0.5° C. This necessitates large equipment for temperature stabilization, unfavorably leading to a loss of a variety of advantages (size reduction, non-inclusion of an active element, etc.) of the channeled spectroscopic polarimetry.

Chapter 3: Solution Against Variations in Reference Phase Function

The reference phase functions φ₁(σ) and φ₂(σ) (depending not upon the sample but only upon parameters of the polarimeter) included in the channeled spectrum vary by a variety of factors, which becomes a major contributor to an error. In consideration of this respect, in one embodiment of the present invention, the channeled spectroscopic polarimeter is provided with a function capable of calibrating the reference phase functions φ₁(σ) and φ₂(σ) of each vibration component during measurement (concurrently with measurement) (cf. FIGS. 8 to 10).

3.1 Method for Calibration “During Measurement” (No. 1)

The calibration method described in Section 1.6 was a method for calibration “prior to measurement”. As opposed to this, in the following section, a method for calibration “during measurement” is shown.

3.1.1 Basic Idea

The amplitude and the phase obtained in Step 2 in Chapter 1 during measurement (when light in an unknown SOP is incident on the channeled spectroscopic polarimeter) is shown again below. [Mathematical Expression 24] $\begin{matrix} {{{{\bullet Component}\lbrack 1\rbrack}\left( {{low}\quad{frequency}\quad{component}} \right){\text{:}\quad\lbrack{Amplitude}\rbrack}\quad{a_{0}(\sigma)}} = {{\frac{1}{2}{m_{0}(\sigma)}{{M_{0}(\sigma)}\quad\lbrack{Phase}\rbrack}\quad{\delta_{0}(\sigma)}} = 0}} & \left( {3.1a} \right) \\ {{{{\bullet Component}\lbrack 2\rbrack}\left( {{central}\quad{period}\quad{1/L_{-}}} \right){\text{:}\quad\lbrack{Amplitude}\rbrack}\quad{a_{-}(\sigma)}} = {{\frac{1}{4}{m_{-}(\sigma)}{{{M_{23}(\sigma)}}\quad\lbrack{Phase}\rbrack}\quad{\delta_{-}(\sigma)}} = {{\phi_{2}(\sigma)} - {\phi_{1}(\sigma)} + {\arg\left\{ {M_{23}(\sigma)} \right\}}}}} & \left( {3.1b} \right) \\ {{{{\bullet Component}\lbrack 3\rbrack}\left( {{central}\quad{period}\quad{1/L_{2}}} \right){\text{:}\quad\lbrack{Amplitude}\rbrack}\quad{a_{2}(\sigma)}} = {{\frac{1}{2}{m_{2}(\sigma)}{{M_{1}(\sigma)}\quad\lbrack{Phase}\rbrack}\quad{\delta_{2}(\sigma)}} = {\phi_{2}(\sigma)}}} & \left( {3.1c} \right) \\ {{{{\bullet Component}\lbrack 4\rbrack}\left( {{central}\quad{period}\quad{1/L_{+}}} \right){\text{:}\quad\lbrack{Amplitude}\rbrack}\quad{a_{+}(\sigma)}} = {{\frac{1}{4}{m_{+}(\sigma)}{{{M_{23}(\sigma)}}\quad\lbrack{Phase}\rbrack}\quad{\delta_{+}(\sigma)}} = {{\phi_{2}(\sigma)} + {\phi_{1}(\sigma)} - {\arg\left\{ {M_{23}(\sigma)} \right\}} + \pi}}} & \left( {3.1d} \right) \end{matrix}$

Here, all needed for obtaining the four spectrometric quasi-Stokes parameters are found to be:

-   -   “amplitude” of Component [1]→M₀(σ)     -   “amplitude” and “phase” of one of Component [2] and Component         [4]→M₂(σ) and M₃(G)     -   “amplitude” of Component [3]→M₁(σ)         It is found that the remaining ones as follows are not used for         demodulation of the spectrometric quasi Stokes parameters.     -   “phase” of Component [3]     -   “amplitude” and “phase” of the remaining one of Component [2]         and component [4]

The present inventors and the like found it possible to obtain not only the four spectrometric quasi-Stokes parameters but also the “reference phase functions (φ₁(σ), φ₂(σ), etc.)” all at once through use of the remaining component. In this method, it is meant that calibration can be concurrently performed in the midst of measurement without particular input of light in a known SOP.

3.1.2 Preparation

In order to use the “calibration method during measurement”, the following prior preparation is necessary.

-   -   The reference amplitude functions m₀(σ), m⁻(σ), m₂(σ), and m₊(σ)         are subjected to pre-calibration (cf. FIG. 9)

Since the following method is effective only on the reference phase function, any one of the methods descried in Section 1.6 is to be performed as for the reference amplitude function. It is to be noted that the fluctuations in the reference amplitude function during measurement typically have considerably small magnitude, and are ignorable in many cases. Namely, in contrast to the reference phase function, there is generally almost no need for recalibration of the reference amplitude function during measurement.

-   -   As for the reference phase function, the pre-calibration is not         necessarily required. However, a ratio between φ₁(σ) and φ₂(σ)         needs to be obtained in advance.

Example 1

when the retarders R₁ and R₂ are made of the same medium, the ratio between φ₁(σ) and φ₂(σ) is determined from a ratio between thicknesses of the two retarders.

Example 2

By further pre-calibration of the reference phase function, the ratio between φ₁(σ) and φ₂(σ) is determined (this ratio may be considered not to change during measurement).

Note here that, in cases including a case where the ratio between the retarders R₁ and R₂ changes (e.g. temperatures of the two retarders are different) during measurement, a method described below cannot be used.

3.1.3 Actual Calibration Method

Based upon the above-mentioned idea, a method for actual calibration is described below.

A. Method for Obtaining Reference Phase Function (σ) from Vibration Component [3]

By noting only Vibration component [3], the amplitude and the phase thereof are shown again as follows. [Mathematical Expression 25] $\begin{matrix} \left\{ \begin{matrix} \lbrack{Amplitude}\rbrack & {{a_{2}(\sigma)} = {\frac{1}{2}{m_{2}(\sigma)}{M_{1}(\sigma)}}} \\ \lbrack{Phase}\rbrack & {{\delta_{2}(\sigma)} = {\phi_{2}(\sigma)}} \end{matrix} \right. & (3.2) \end{matrix}$ What needs to be noted here is that the phase δ₂(σ) of this component is one (φ₂(σ)) of the reference phase functions (itself). Namely, when the phase δ₂(σ) of Component [3] is measured, one (φ₂(σ)) of the reference phase functions is immediately determined according to the following expression. φ₂(σ)=δ₂(σ)  (3.3)

This relational expression is constantly satisfied regardless of a polarization characteristic of a measurement sample, meaning that one of the reference phase functions can be immediately obtained from a measured value, even from a channeled spectrum by any kind of sample. This is a calibration method that can be performed utterly concurrently during measurement, and in the case of “using a sample having a known polarization characteristic”, there is no need for performing calibration “prior to measurement or after discontinuation of measurement” as in (Section 1.6). However, it should be noted that at this time, the condition of observing Component [3] at a sufficient SN ratio needs to be satisfied (cf. later-described C).

It is to be noted that, when the “complex representation” is obtained in place of the “paired amplitude and phase” in Step 2 of the “procedure for demodulating spectrometric quasi-Stokes parameters” in Section 1.5, a calculation method, rewritten from the above and described below, may be applied.

From Expression (1.14)b, δ₂(σ) has the following relation with the complex representation F₂(σ) of Component [3]. δ₂(σ)=arg [F ₂(σ(]  (3.4) Therefore, the reference phase function φ₂(σ) can be obtained from the complex representation of Component [3] according to the following expression. φ₂(σ)=arg [F ₂(σ)]  (3.5) It should be noted that what is needed at the time of complex representation is not the reference phase function φ₂(σ) but the reference complex function K₂(σ). Since there is a relation between these two functions as expressed by Expression (1.18c), once φ₂(σ) is determined, K₂(σ) can also be determined (this will be later described in details in F).

B. Method for Obtaining Reference Phase Function (σ) from a Plurality of Vibration Components (Paired [2] and [4], etc.)

The respective phases of Vibration Components [2] and [4] are again shown as follows.

Phase of Component [2]: δ⁻(σ)=φ₂(σ)−φ₁(σ)+arg {M ₂₃(σ)}  (3.6a)

Phase of Component [4]: δ₊(σ)=φ₂(σ)+φ₁(σ)−arg {M ₂₃(σ)}+π  (3.6b) When the one phase is added to the other, φ₁(σ) and arg {M₂₃(σ)} are canceled out, whereby it is found that the following expression is satisfied. [Mathematical Expression 26] $\begin{matrix} {{\phi_{2}(\sigma)} = {{\frac{1}{2}\left\{ {{\delta_{-}(\sigma)} + {\delta_{+}(\sigma)}} \right\}} - \frac{\pi}{2}}} & (3.7) \end{matrix}$

The right side of the above expression means that one (φ₂(σ)) of the reference phase functions can be obtained by taking an average of the phases of Vibration components [2] and [4]. Similarly to Method A, this relational expression can also be constantly satisfied regardless of the SOP of the sample, meaning that one of the reference phase functions can be immediately obtained from a measured value, even from a channeled spectrum by any kind of sample.

Namely, similarly to the case of Method A, this is a “calibration method that can be performed utterly concurrently during measurement”, and in the case of “using a sample having a known polarization characteristic”, there is no need for performing calibration “prior to measurement or after discontinuation of measurement” as in (Section 1.6). However, it should be noted that the condition of observing Components [2] and [4] at a sufficient SN ratio needs to be satisfied (cf. later-described C).

Here, similarly to the case of Method A, a calculation method is described which is used in the case of obtaining the “complex representation” in place of the “paired amplitude and phase” in Step 2 of Section 1.5.

From Expression (1.14b), δ⁻(σ) and δ₊(σ) have the following relations with the complex representations F⁻(σ) and F₊(σ) of Components [2] and [4]. δ⁻(σ)=arg [F ⁻(σ)]  (3.8a) δ₊(σ)=arg [F ₊(σ)]  (3.8b)

Therefore, the reference phase function φ₂(σ) can be obtained from the complex representations of the two components as follows. [Mathematical Expression 27] $\begin{matrix} {{\phi_{2}(\sigma)} = {{\frac{1}{2}\left\{ {{\arg\left\lbrack {F_{-}(\sigma)} \right\rbrack} + {\arg\left\lbrack {F_{+}(\sigma)} \right\rbrack}} \right\}} - \frac{\pi}{2}}} & (3.9) \end{matrix}$ Or, the following expression obtained by rewriting the above expression to a simple formula of the complex function may be applied. [Mathematical Expression 28] $\begin{matrix} {{\phi_{2}(\sigma)} = {\frac{1}{2}{\arg\left\lbrack {{- {F_{-}(\sigma)}}{F_{+}(\sigma)}} \right\rbrack}}} & (3.10) \end{matrix}$

In the optically system (channeled spectroscopic polarimeter) in FIG. 2, an obtained spectrum includes another component having a different period as described above except for the case where the angle formed between the retarder R₂ and the polarizer P is not 45°.

As seen from Expression (1.10), the phase of this component is δ₁(σ)=φ₁(σ) −arg {M₂₃(σ)}”, and similar to the phase terms of above Vibration components [2] and [4]. Hence, even when the above component is combined with Components [2] and [4] (or replaced by one of them), it is possible to calibrate φ₂(σ) in the same manner.

C. Combination of A and B

The two methods (Method A and Method B) described above are methods in which one (φ₂(σ)) of the reference phase functions can be calibrated utterly concurrently during measurement. However, the vibration components used are different between the two methods. What should be concerned here is that the amplitude of Vibration component [3] used in Method A is proportional to M₁(σ), while the amplitudes of Vibration Components [2] and [4] used in Method B are proportional to the following. |M ₂₃(σ)|=√{square root over (M ₂ ²(σ)+M ₃ ²(σ))}  [Mathematical Expression 29]

Since the polarization characteristic of the sample is unknown, there is no guarantee that the spectrometric quasi-Stokes parameters are constantly sufficiently large for phase measurement for each component. For example, when M₁(σ) is small, obtaining φ₂(σ) by Method A using the phase of this component might result in occurrence of a large error. For solving this problem, adaptive combination of Methods A and B is desired. Specifically, a value of φ₂(σ) with more certainty can be obtained by selecting, or weighting up and balancing, results of the two methods.

D. Combination of A and B (No. 2)

One idea for efficiently combining A and B is shown below. This is a method in which direct calculation is possible without particular separation by case. It should be noted that, in this part (Method D), three complex representation functions F⁻(σ) F₂(σ) and F⁻(σ) of Components [2] to [4] are used for calculation. When a calculation is to be made from the “paired amplitude and phase” of each vibration component, this pair may once be changed to the “complex representation” according to Expression (1.13), and then the following calculation procedure may be performed.

As a preparation for explaining this method, first, the following two expressions are derived and the properties thereof are described. By transforming Expression (3.5), the following expression can be obtained. 2φ₂(σ)=arg[F ₂ ²(σ)]  (3.11) Meanwhile, by doubling both sides of Expression (3.10), the following expression can be obtained. 2φ₂(σ)=arg[−F ⁻(σ)F ₊(σ)]  (3.12) It is found from the comparison between the above two expressions that the complex function in the brackets on the right side of each of the expressions has the same argument 2φ₂(σ).

“Appropriate weighting functions α(σ) and β(σ) which have the same argument” were respectively multiplied by the above two complex functions, and then the obtained two terms were added together.

[Mathematical Expression 30] α(σ)[F ₂ ²(σ)]β(σ) [−F⁻(σ)F ₊(σ)]  (3.13) The argument of this expression is constantly equivalent to 2φ₂(σ)+argα(σ). It is found that, through use of the properties thus described, φ₂(σ) can be obtained according to the following expression even if one of M₁ and M₂₃ decreases. [Mathematical Expression 31] $\begin{matrix} \begin{matrix} {{\phi_{2}(\sigma)} = {{\frac{1}{2}\arg\left\{ {{{\alpha(\sigma)}\left\lbrack {F_{2}^{2}(\sigma)} \right\rbrack} + {{\beta(\sigma)}\left\lbrack {{- {F_{-}(\sigma)}}{F_{+}(\sigma)}} \right\rbrack}} \right\}} -}} \\ {\frac{1}{2}{\arg\left\lbrack {\alpha(\sigma)} \right\rbrack}} \end{matrix} & (3.14) \end{matrix}$

There are a variety of specific ways to select α(σ) and β(σ).

The simplest way to select α(σ) and β(σ) is making the two functions the same constant (1). In this case, an expression for obtaining the reference phase function φ₂(σ) is shown below. [Mathematical Expression 32] $\begin{matrix} {{\phi_{2}(\sigma)} = {\frac{1}{2}\arg\left\{ {\left\lbrack {F_{2}^{2}(\sigma)} \right\rbrack + \left\lbrack {{- {F_{-}(\sigma)}}{F_{+}(\sigma)}} \right\rbrack} \right\}}} & (3.15) \end{matrix}$

E. Calculation of φ₁(σ)

Since fluctuations in φ₁(σ) are considered to be similar to those in φ₂(σ), it is possible to obtain φ₁(σ) from a measured value of φ₂(σ) by proportional calculation (e.g. by using a thickness ratio).

F. Calculation of Reference Complex Function

In the demodulation in Step 2 of the “procedure for demodulating spectrometric quasi-Stokes parameters” in Section 1.5, when (not the “paired amplitude and phase” but) the “complex representation” is obtained, what are needed ultimately in the operation of Step 3 for obtaining the spectrometric quasi-Stokes parameters are not the reference phase functions φ₁(σ) and φ₂(σ) but the reference complex functions K₀(σ), K⁻(σ), K₂(σ), and K₊(u). However, these can also be immediately obtained through use of the relations of Expressions (1.18a) to (1.18d) if the reference phase functions φ₁(σ) and φ₂(σ) have been obtained by the procedures up to E above.

3.2 Method for Calibrating Reference Phase Function “During Measurement” (No. 2) 3.2.1 Basic Idea

In the same idea as described in the previous section 3.1, “only a variation” in reference phase function can be obtained.

In the previous method (in the previous section 3.1), the “reference amplitude function” was obtained in the pre-calibration, and it was not particularly necessary to obtain the “reference phase function”. However, as appeared from Section 3.2, those two functions can be calibrated almost concurrently. It is thus possible to obtain an “initial value of the reference phase function in pre-calibration” so as to only track the variation thereof during measurement.

Advantages in this case include the following.

-   -   Slightly additional phase displacement part which might be         generated due to characteristics of the spectrometer or the         signal processing system can be removed.     -   Burdensome phase unwrapping is not necessary.     -   Since a variation in phase difference itself is small, a dynamic         range in calculation can be made small. Further, as a result of         this, a calculation error can be relatively made small in many         cases.

Accordingly, “obtaining only the variation in reference phase function” has its own meaning.

The following is described as supplement for the foregoing explanation. As shown in FIG. 10, the two methods have different factors of an error in the case of obtaining φ₁(c from φ₂(σ). Namely, as shown in FIG. 10A, it is necessary to perform phase unwrapping in the case of obtaining φ₁(σ) from φ₂(σ). This phase unwrapping is a major factor of the error. Especially when a period is at high frequency as compared with sampling, noise is included in the period, or the like, wrong phase unwrapping might be performed. With wrong phase wrapping performed, an error becomes an integer multiple of 2π, leading to calculation of a wrong phase. Further, this error affects a broad wavenumber region. This difference is essentially caused by that a solution of an arg operator (or an arctan operator) for obtaining an argument has ambiguity by the integer multiple of 2π. As opposed to this, as shown in FIG. 10B, it is not necessary in obtaining Δφ₁(σ) from Δφ₂(σ) to perform phase unwrapping since the variation Δφ₂(σ) from the initial value of the reference phase function is small. This allows the measurement error to be relatively small.

3.2.2 Preparation

The use of the “calibration method during measurement” is based upon the premise of pre-calibration of both the “reference amplitude function” and the “reference phase function” prior to measurement. It is to be noted that, as for the phase, highly accurately obtaining the reference phase function is not necessary since a phase value can be corrected later by the variance and the measurement error.

3.2.3 Actual Calibration Method

The basic idea on the calibration method is completely the same as in Section 3.1. There thus exists a calculation method corresponding to all A to E described in Section 3.1.3. Hence, in this section, an idea different from that in the previous section is shown and mathematical expressions are mainly cited in the following description.

First, a couple of symbols are defined. The reference phase functions obtained by the pre-calibration are defined as φ₁ ^((i))(σ) and φ₂ ^((i))(σ). Reference complex functions corresponding to these reference phase functions are expressed as follows according to Expressions (1.18a) to (1.18d). [Mathematical Expression 33] $\begin{matrix} {{K_{0}^{(i)}(\sigma)} = {\frac{1}{2}{m_{0}(\sigma)}}} & \left( {3.16a} \right) \\ {{K_{-}^{(i)}(\sigma)} = {\frac{1}{8}{m_{-}(\sigma)}\exp\quad{{\mathbb{i}}\left\lbrack {{\phi_{2}^{(i)}(\sigma)} - {\phi_{1}^{(i)}(\sigma)}} \right\rbrack}}} & \left( {3.16b} \right) \\ {{K_{2}^{(i)}(\sigma)} = {\frac{1}{4}{m_{2}(\sigma)}\exp\quad{\mathbb{i}}\quad{\phi_{2}^{(i)}(\sigma)}}} & \left( {3.16c} \right) \\ {{K_{+}^{(i)}(\sigma)} = {{- \frac{1}{8}}{m_{+}(\sigma)}\exp\quad{{\mathbb{i}}\left\lbrack {{\phi_{2}^{(i)}(\sigma)} + {\phi_{1}^{(i)}(\sigma)}} \right\rbrack}}} & \left( {3.16d} \right) \end{matrix}$ Assuming that the reference phase functions changed during measurement as follows. φ₁(σ)=φ₁ ^((i))(σ)+Δφ₁(σ)  (3.17a) φ₂(σ)=φ₂ ^((i))(σ)+Δφ₂(σ)  (3.17b) Below described are methods for obtaining the variations Δφ₁(σ) and Δφ₂(σ) of the reference phase functions or changes in the reference complex functions corresponding to those variations.

A. Method for Obtaining Reference Phase Function φ₂ from Vibration Component [3]

As described in Method A in the previous section, the phase of Component [3] is expressed as follows. δ₂(σ)=φ₂(σ)=φ₂ ^((i))(σ)+Δφ₂(σ)  (3.18) Here, the variation in φ₂(σ) can be obtained as Δφ₂(σ)=δ₂(σ)−φ₂ ^((i))(σ)  (3.19) Namely, this means that, once the phase δ₂ of Component [3] is measured, one of the variations Δφ₂(σ) in reference phase functions can be immediately determined.

It is to be noted that in Step 2, when not the “paired amplitude and phase” but the “complex representation” is to be obtained, the variation is obtained according to the following expressions. [Mathematical Expression 34] $\begin{matrix} {{\delta_{2}(\sigma)} = {\arg\left\lbrack {F_{2}(\sigma)} \right\rbrack}} & \left( {3.20a} \right) \\ {{\phi_{2}^{(i)}(\sigma)} = {\arg\left\lbrack {K_{2}^{(i)}(\sigma)} \right\rbrack}} & \left( {3.20b} \right) \\ {From} & \quad \\ {{\Delta\quad{\phi_{2}(\sigma)}} = {{\arg\left\lbrack {F_{2}(\sigma)} \right\rbrack} - {\arg\left\lbrack {K_{2}^{(i)}(\sigma)} \right\rbrack}}} & (3.21) \\ {,{or}} & \quad \\ {{\Delta\quad{\phi_{2}(\sigma)}} = {\arg\left\lbrack \frac{F_{2}(\sigma)}{K_{2}^{(i)}(\sigma)} \right\rbrack}} & (3.22) \end{matrix}$

B. Method for Obtaining Reference Phase Function φ₂(σ) from a Plurality of Vibration Components (Paired [2] and [4], etc.)

In the method for obtaining the variation in φ₂(σ) from the phase of each of Vibration components [2] and [4], the variation is obtained according to the following expression. [Mathematical Expression 35] $\begin{matrix} {{\Delta\quad{\phi_{2}(\sigma)}} = {\left\lbrack {{\frac{1}{2}\left\{ {{\delta_{-}(\sigma)} + {\delta_{+}(\sigma)}} \right\}} - \frac{\pi}{2}} \right\rbrack - {\phi_{2}^{(i)}(\sigma)}}} & (3.23) \end{matrix}$

When not the “paird amplitude and phase” but the “complex representation” is to be obtained, the variation is obtained according to the following expression. [Mathematical Expression 36] $\begin{matrix} {{\Delta\quad{\phi_{2}(\sigma)}} = {\frac{1}{2}\begin{Bmatrix} {{\arg\left\lbrack {F_{-}(\sigma)} \right\rbrack} + {\arg\left\lbrack {F_{+}(\sigma)} \right\rbrack} -} \\ {{\arg\left\lbrack {K_{-}^{(i)}(\sigma)} \right\rbrack} - {\arg\left\lbrack {K_{+}^{(i)}(\sigma)} \right\rbrack}} \end{Bmatrix}}} & (3.24) \end{matrix}$ Or, the following expressions obtained by rewriting the above expression using a simple formula of the complex function may be applied. [Mathematical Expression 37] $\begin{matrix} {{\Delta\quad{\phi_{2}(\sigma)}} = {\frac{1}{2}\left\{ {{\arg\left\lbrack \frac{F_{-}(\sigma)}{K_{-}^{(i)}(\sigma)} \right\rbrack} + {\arg\left\lbrack \frac{F_{+}(\sigma)}{K_{+}^{(i)}(\sigma)} \right\rbrack}} \right\}}} & (3.25) \\ {,{or}} & \quad \\ {{\Delta\quad{\phi_{2}(\sigma)}} = {\frac{1}{2}{\arg\left\lbrack {\frac{F_{-}(\sigma)}{K_{-}^{(i)}(\sigma)}\frac{F_{+}(\sigma)}{K_{+}^{(i)}(\sigma)}} \right\rbrack}}} & (3.26) \end{matrix}$ It is to be noted that, as noted at the end of Section 3.1.3, the same idea as shown above can be applied to the case of using another term.

C. Combination of A and B

As in the case described in the previous section, adaptive combination of Methods A and B is effective also in the case of obtaining only the “variation” in reference phase functions. It should be noted that a description of the combination is completely the same as that in the previous section and it is thus omitted.

D. Combination of A and B (No. 2)

One of desired mathematical expressions in the case of obtaining only the variation is as follows. [Mathematical Expression 38] $\begin{matrix} {{\alpha(\sigma)} = \left\lbrack \frac{1}{K_{2}^{(i)}(\sigma)} \right\rbrack^{2}} & \left( {3.27a} \right) \\ {{\beta(\sigma)} = {- \frac{1}{{K_{-}^{(i)}(\sigma)}{K_{+}^{(i)}(\sigma)}}}} & \left( {3.27b} \right) \end{matrix}$ Since arg[α(σ)]=arg[β(σ)]=2φ₂(σ) in the above expressions, the variation can be obtained as follows. [Mathematical Expression 39] $\begin{matrix} {{\Delta\quad{\phi_{2}(\sigma)}} = {\frac{1}{2}\arg\left\{ {\left\lbrack \frac{F_{2}(\sigma)}{K_{2}^{(i)}(\sigma)} \right\rbrack^{2} + {\frac{F_{-}(\sigma)}{K_{-}^{(i)}(\sigma)}\frac{F_{+}(\sigma)}{K_{+}^{(i)}(\sigma)}}} \right\}}} & (3.28) \end{matrix}$

E. Calculation of Δφ₁(σ)

Fluctuations in Δφ₁(a are considered similar to those in Δ₂(σ). It is thus possible to obtain Δ₁ (σ) from a measured value of Δ₂(σ) by a comparative calculation using, for example, a thickness ratio.

F. Calculation of Reference Complex Function

In the demodulation of each vibration component in Step 2, when not the “paired amplitude and phase” but the “complex representation” is obtained, what are needed ultimately in obtaining the spectrometric quasi-Stokes parameters (operation of Step 3) are not the reference phase functions φ₁(σ) and φ₂(σ) but the reference complex functions K₀(σ), K⁻(σ), K₂(σ), and K₊(σ).

If the reference phase function variations Δφ₁(σ) and Δφ₂(σ) have been obtained by the procedures up to E above, the reference complex functions can be immediately obtained as follows.

[Mathematical Expression 40] K ₀(σ)=K₀ ^((i))(σ)  (3.29a) K⁻(σ)=K ⁻ ^((i))(σ)e ^(i[Δφ) ² ^((σ)−Δφ) ¹ ^((σ)])  (3.29b) K ₂(σ)=K ₂ ^((i))(σ)e ^(iΔφ) ² ^((σ))  (3.29c) K ₊(σ)=K ₊ ^((i))(σ)e ^(i[Δφ) ² ^((σ)+Δφ) ¹ ^((σ)])  (3.29d)

Chapter 4: Specific Embodiment of Present Invention

4.1 Case of Spectroscopic Polarimetry Performed by Reflecting Light on Sample

An embodiment of the optical system in the case of measuring a spectropolarization characteristic of a sample by reflecting light on the sample is described in detail with reference to FIGS. 11 to 13. In this case, as shown in FIG. 11 and FIG. 12, the optical system includes the light source 7, the polarizer P, the retarder R2, the retarder R1, the analyzer A, and the spectroscope 8. It is to be noted that reference symbol B denotes a sample on which light is reflected. Further, light emitted from the light source 7 is transmitted through the polarizer P, the retarder R2 and the retarder R1 in this order. The light is then incident on the sample B in a slanting direction, to be reflected thereon. Thereafter, the light is transmitted through the analyzer A and then received in the spectroscope 8. Attention should be paid that this device is configured by arranging an optical element on the light source side with respect to the sample in the configuration of the conventional optical system as shown in FIG. 1B. Further, the orientation of the transmission axis of the polarizer P agrees with the orientation of the principal axis of the retarder R1, and the orientations of the fast axes of the retarder R1 and the retarder R2 are inclined at an angle of −45° from each other. θ shows the azimuth angle of the transmission axis of the analyzer with respect to the fast axis of the retarder R1. Further, the incident plane of the light agrees with the orientation of the fast axis of the retarder R1. It is to be noted that the device including the polarizer P, the retarder R2, the retarder R1 and the analyzer A is referred to as a channeled spectrometer unit.

At the same time, according to the optical system as shown in FIG. 11 and FIG. 12, a stable measurement can be realized since the incident direction of the wave surface of light that transmits through the retarders R1 and R2 is not susceptible to the sample, as shown in FIG. 13. That is, it is possible to solve the problem of the retardation change between the time of calibrating the retardation of the retarder and the time of measuring the sample due to changes in distance and direction of a light ray that passes through the retarder.

Further, an ellipsometric parameter and the like are determined from the light incident on the spectroscope 8. A procedure for such determination is described below.

Here, an arctangent of a rate of change in amplitude rate of a p-polarized light (light having an SOP where the polarizing direction runs parallel with the incident plane) and an s-polarized light (light having an SOP where the polarizing direction runs vertical to the incident plane) is expressed as ψ(σ), and a phase difference is expressed as Δ(σ). Then, a Mueller matrix of an isotropic medium can be described as follows. [Mathematical Expression 41] $\begin{matrix} {{{{M(\sigma)} =}\quad}\begin{bmatrix} 1 & {{- \cos}\quad 2{\Psi(\sigma)}} & 0 & 0 \\ {{- \cos}\quad 2{\Psi(\sigma)}} & 1 & 0 & 0 \\ 0 & 0 & {\sin\quad 2{\Psi(\sigma)}\cos\quad{\Delta(\sigma)}} & {\sin\quad 2{\Psi(\sigma)}\sin\quad{\Delta(\sigma)}} \\ 0 & 0 & {{- \sin}\quad 2{\Psi(\sigma)}\sin\quad{\Delta(\sigma)}} & {\sin\quad 2{\Psi(\sigma)}\cos\quad{\Delta(\sigma)}} \end{bmatrix}} & (4.1) \end{matrix}$ It is thereby considered that ψ(σ) and Δ(σ) are obtained as ellipsometric parameters.

Here, the following expressions are satisfied according to Expressions (1.5a) to (1.5d). [Mathematical Expression 42] $\begin{matrix} {{M_{0}(\sigma)} = {\frac{1}{2}{{P_{0}(\sigma)}\left\lbrack {1 - {\cos\quad 2{\Psi(\sigma)}\cos\quad 2\quad\theta}} \right\rbrack}}} & \left( {4.2a} \right) \\ {{M_{1}(\sigma)} = {{- \frac{1}{2}}{{P_{0}(\sigma)}\left\lbrack {{{- \cos}\quad 2{\Psi(\sigma)}} + {\cos\quad 2\quad\theta}} \right\rbrack}}} & \left( {4.2b} \right) \\ {{M_{2}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\sin\quad 2{\Psi(\sigma)}\cos\quad{\Delta(\sigma)}\sin\quad 2\theta}} & \left( {4.2c} \right) \\ {{M_{3}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\sin\quad 2{\Psi(\sigma)}\sin\quad{\Delta(\sigma)}\sin\quad 2\theta}} & \left( {4.2d} \right) \end{matrix}$ Here, assuming that θ=45°, the following expressions are satsfied. [Mathematical Expression 43] $\begin{matrix} \begin{matrix} {{M_{0}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}}} \\ {{M_{1}(\sigma)} = {{- \frac{1}{2}}{P_{0}(\sigma)}\cos\quad 2{\Psi(\sigma)}}} \\ {{M_{2}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\sin\quad 2{\Psi(\sigma)}\cos\quad{\Delta(\sigma)}}} \\ {{M_{3}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\sin\quad 2{\Psi(\sigma)}\sin\quad{\Delta(\sigma)}}} \end{matrix} & (4.3) \end{matrix}$

Since M₀(σ) to M₃(σ) which are obtained by demodulation processing include three unknowns; the spectrum P₀(σ) and the ellipsometric parameters ψ(σ) and Δ(σ), as expressed in the above expressions, it is possible to determine the ellipsometric parameters ψ(σ) and Δ(σ).

Moreover, one example of a case where modulating at least one “component” shown in Expressions (1.11a) to (1.11d) gives a useful application is described here.

For example, it is assumed that there is a sample, only a value of Δ(σ) of which is desired to be measured. In such a case, by obtaining only M₂₃(σ), two equations of M₂(σ) and M₃(σ) are obtained. Solving these equations can lead to calculation of the ellipsometric parameter Δ(σ) as follows. [Mathematical Expression 44] $\begin{matrix} {{\Delta\quad(\sigma)} = {\tan^{- 1}\left\{ \frac{M_{3}(\sigma)}{M_{2}(\sigma)} \right\}}} & (4.4) \end{matrix}$

As described before, M₀(σ) is obtained from the first spectral intensity, M₁(σ) is from the third spectral intensity, and M₂(σ) and M₃(σ) are from at least one of the second spectral intensity, the fourth spectral intensity, and the fifth spectral intensity. Therefore, with reference to Expression (4.3), ψ(σ) which is obtained from the expressions of M₀(σ) and M₁(σ) can be obtained from the first spectral intensity and the third spectral intensity. In this case, even if the reference phase function is to be calibrated during measurement, calibration of only φ₂(σ) sufficiently meets the need. ψ(σ) can also be obtained from at least one of the first spectral intensity and the third spectral intensity and at least one of the second spectral intensity, the fourth spectral intensity and the fifth spectral intensity. Since Δ(σ) can be obtained from M₂(σ) and M₃(σ), at least one of the second spectral intensity, the fourth spectral intensity and the fifth spectral intensity can be obtained.

4.2 Case of Spectroscopic Polarimetry Performed by Transmitting Light Through Sample

An embodiment of the optical system in the case of measuring a spectropolarization characteristic of a sample by transmitting light through the sample is described in detail with reference to FIGS. 16 to 18. In this case, as shown in FIG. 16 and FIG. 17, the optical system includes the light source 7, the polarizer P, the retarder R2, the retarder R1, the analyzer A, and the spectroscope 8. It is to be noted that reference symbol C denotes a sample through which light is transmitted. Further, light emitted from the light source 7 is transmitted through the polarizer P, the retarder R2 and the retarder R1 in this order. The light is then incident on the sample C in a vertical direction, to be transmitted therethrough. Thereafter, the light is transmitted through the analyzer A and then received in the spectroscope 8. Note should be taken that this device is configured by arranging an optical element on the light source side with respect to the sample in the configuration of the conventional optical system as shown in FIG. 1B. Here, the orientation of the transmission axis of the polarizer P agrees with the orientation of the fast axis of the retarder R1, while the orientations of the fast axes of the retarder R1 and the retarder R2 are inclined at an angle of 45° from each other. θ shows the azimuth angle of the transmission axis of the analyzer A with respect to the fast axis of the retarder R1.

Moreover, according to the optical system as shown in FIG. 16 and FIG. 17, since the incident direction of the wave surface of light that transmits through the retarders R1 and R2 is not susceptible to the inclination characteristic (cf. FIG. 18A) and the surface state (cf. FIG. 18B) of the sample, a stable measurement can be realized without restrictions of being unable to measure the surface shape of the sample and a disturbing substance (living body, etc.). Namely, it is possible to solve the problem of the retardation change between the time of calibrating the retardation of the retarder and the time of measuring the sample due to changes in distance and direction of a light ray that passes through the retarder.

Below described is a procedure for obtaining an azimuth angle R of a birefringent axis of a sample (birefringent medium) with respect to the fast axis of the retarder R1, as well as a retardation δ(σ), from light incident on the spectroscope 8.

When the azimuth angle of the birefringent axis of the sample (birefringent medium) with respect to the fast axis of the retarder R1 is R, and the retardation of the sample (birefringent medium) is δ(σ), a Mueller matrix for expressing the sample (birefringent medium) is described as follows. [Mathematical Expression 45] $\begin{matrix} {{M(\sigma)} = \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & {{\cos^{2}2R} + \quad{\cos\quad{\delta(\sigma)}\sin^{2}2R}} & {\cos\quad 2R\quad\sin\quad 2{R\left( {1 - {\cos\quad{\delta(\sigma)}}} \right)}} & {{- \sin}\quad{\delta(\sigma)}\sin\quad 2R} \\ 0 & {\cos\quad 2R\quad\sin\quad 2{R\left( {1 - {\cos\quad{\delta(\sigma)}}} \right)}} & {{\sin^{\quad 2}2R} + {\cos\quad{\delta(\sigma)}\cos^{\quad 2}2R}} & {\cos\quad 2R\quad\sin\quad{\delta(\sigma)}} \\ 0 & {\sin\quad{\delta(\sigma)}\sin\quad 2R} & {{- \cos}\quad 2R\quad\sin\quad{\delta(\sigma)}} & {\cos\quad{\delta(\sigma)}} \end{bmatrix}} & (5.1) \end{matrix}$ Here, the following expressions are satisfied according to Expressions (1.5a) to (1.5d). [Mathematical Expression 46] $\begin{matrix} {{M_{0}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}}} & \left( {5.2a} \right) \\ \begin{matrix} {{M_{1}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\left\{ {{\cos\quad 2R\quad\cos\left( {{2R} - {2\quad\theta}} \right)} +} \right.}} \\ \left. {\cos\quad{\delta(\sigma)}\sin\quad 2R\quad{\sin\left( {{2R} - {2\quad\theta}} \right)}} \right\} \end{matrix} & \left( {5.2b} \right) \\ \begin{matrix} {{M_{2}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\left\{ {{\sin\quad 2R\quad{\cos\left( {{2R} - {2\quad\theta}} \right)}} -} \right.}} \\ \left. {\cos\quad{\delta(\sigma)}\cos\quad 2R\quad{\sin\left( {{2R} - {2\quad\theta}} \right)}} \right\} \end{matrix} & \left( {5.2c} \right) \\ {{M_{3}(\sigma)} = {{- \frac{1}{2}}{P_{0}(\sigma)}\quad\sin\quad{\delta(\sigma)}\quad\sin\quad\left( {{2R} - {2\quad\theta}} \right)}} & \left( {5.2d} \right) \end{matrix}$ Here, assuming that 0=45°, the following expressions are satisfied. [Mathematical Expression 47] $\begin{matrix} {{M_{0}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}}} & \left( {5.3a} \right) \\ {{M_{1}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\cos\quad 2R\quad\sin\quad 2{R\left( {1 - {\cos\quad\delta\quad(\sigma)}} \right)}}} & \left( {5.3b} \right) \\ {{M_{2}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\left( {{\sin^{2}2R} + {\cos\quad{\delta(\sigma)}\cos^{2}2R}} \right)}} & \left( {5.3c} \right) \\ {{M_{3}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\cos\quad 2R\quad\sin\quad{\delta(\sigma)}}} & \left( {5.3d} \right) \end{matrix}$

Since M₀(σ) to M₃(σ) which are obtained by demodulation processing include three unknowns: the spectrum P₀(σ) of the light ray, the azimuth angle R of the birefringent axis, and the retardation δ(σ) of the sample (birefringent medium), it is possible to determine the azimuth angle R of the birefringent axis and the retardation δ(σ) of the sample (birefringent medium).

Moreover, there exists a case where what needs to be obtained may be either the azimuth angle R of the birefringent axis or the retardation δ(σ) of the sample (birefringent medium). For example, in the cases of a liquid crystal, a polymer film, and the like, the orientation of the birefringent axis can be determined when the azimuth angle thereof is obtained. Further, since the azimuth angle R can be obtained only from M₁(σ) and M₂(σ) while M₂(σ) is obtained from the second, fourth or fifth spectral intensity, M₃(σ) can be obtained concurrently with M₂(σ). Further, if M₂(σ) and M₃(σ) are known, M₀(σ) (obtained from the first spectral intensity) may be used in place of M₁(σ) (obtained from the third spectral intensity). Eventually, the necessary spectral intensities are either the first or third spectral intensity and any one of the second, fourth and fifth spectral intensities. This also applies to necessary spectral intensities for obtaining δ(σ).

4.3 Pre-Calibration

Next, pre-calibration of the device configuration which was described in Sections 4.1, 4.2 and the like is described with reference to FIG. 14 and FIG. 15.

FIG. 14 shows a device configuration necessary for pre-calibration. This device is comprised of the light source 7, the spectroscope 8, a unit 9 for pre-calibration, a unit 12 for measurement, an optical fiber 10 for use in pre-calibration, and an optical fiber 11 for use in measurement. The unit 12 for measurement includes the polarizer P, the retarder R2, the retarder R1, and an analyzer A2 for measurement. The unit 9 for pre-calibration includes an analyzer A1 for pre-calibration. It is to be noted that the analyzer A1 for pre-calibration is set to have a known polarization angle.

According to this device, at the time of pre-calibration, light is emitted from the light source 7, and transmitted through the polarizer P, the retarder R2 and the retarder R1 in this order. The light is then transmitted through the analyzer A1 for pre-calibration included in the unit 9 for pre-calibration, to be incident on the spectroscope 8 through the optical fiber 10. On the other hand, when an object to be measured such as a sample exists, a measurement is performed by the methods described in Sections 4.1 and 4.2, using the unit 12 for measurement.

What is of importance here is that there is no need for using the same one analyzer “for pre-calibration” and ‘for measurement’. This is because the analyzer has the (disappearing) characteristic of being susceptible to fluctuations due to variations in incident angle of a transmitted light ray, as above described.

Therefore, the unit 9 for pre-calibration (light-reception part) including the analyzer A1 for pre-calibration set to have a known polarization angle becomes movable to place where calibration is readily performed. Thereby, the advantage of being able to perform precalibration in a place which is not a measurement place is obtained. Further, the advantages of time reduction and the like are also obtained simultaneously.

Further, as for the pre-calibration, a problem with the channeled spectroscopic polarimetry described in “Spectroscopic ellipsometry using channeled spectrum”, written by Kazuhiko Oka and Takayuki Katoh, collected papers of lectures in 26th Study Session on Light Wave Sensing Technology (Light wave Sensing Technology Study Session held by Japan Society of Applied Physics, Dec. 19-20, 2000) pp. 107-114 is shown in FIG. 15. In FIG. 15, light is incident from the lower left of the figure. The light is transmitted through the polarizer P, reflected on the sample B, and then transmitted through the retarder R.

In this calibration, it is necessary to apply known linearly polarized light to the channeled spectroscopic polarimeter at the time of calibration. Namely, it is necessary to adjust an SOP of light to be incident on the sample so that the light has the SOP of the known linearly polarized light after reflection on the sample. This has required changing the light to be incident on the sample to the p-polarized light (light having an SOP where the polarizing direction runs parallel with the incident plane) and an s-polarized light (light having an SOP where the polarizing direction runs vertical to the incident plane), limited to the time of calibration. (Light other than the p-polarized light and the s-polarized light undesirably becomes elliptically polarized light after reflection on the sample. Since the SOP of the elliptically polarized light depends upon a refraction index of the sample, surface roughness of the sample, or the like, treating the elliptically polarized light as calibration light might cause generation of a measurement error, which is inconvenient in calibration.) Further, it has been necessary to adjust the polarizer P again so as to have a known rotation angle other than the p-polarized light and the s-polarized light at the time of measurement. This has required a system for adjusting the rotation angle of the polarizer, such as a stage, preventing size reduction in the channeled spectropolarimetric unit.

However, according to the calibration method of the present embodiment, pre-calibration is not accompanied by reflection on the sample as apparent from FIG. 14, thereby giving the advantage of eliminating the need for arranging a polarization angle adjustment system on the light-projection side (retarder side), to permit size reduction in light-projection side unit.

4.4 Case where Object to be Measured Includes Known Polarization Element Other Than Sample

There are some cases where the object to be measured includes a polarization element (e.g. quarter wave plate) having a known polarization characteristic, and FIG. 19 shows a device configuration of an optical system in such cases. In this figure, the optical system includes the light source 7, the polarizer P, the retarder R2, the retarder R1, the analyzer A and the spectroscope 8. The object to be measured includes the sample D and a known polarization element E. Light emitted from the light source 7 is transmitted through the polarizer P, the retarder R2 and the retarder R1 in this order. The light is then reflected on or transmitted through the object to be measured, and transmitted through the analyzer, to be incident on the spectroscope 8.

In the arithmetic processing for obtaining a spectroscopic quasi-Stokes parameter of the sample in this case, the sample and the known sample are regarded and measured as one object to be measured, and the effect of the known polarization element may be eliminated from the obtained equation (spectroscopic quasi-Stokes parameter). In the following, as an example, a case is considered where a quarter wave plate is arranged in a position after the sample such that the slow axis thereof runs parallel with the fast axis of the first retarder. (In FIG. 19, the known polarization element E is regarded as the quarter wave plate.)

The sample and quarter wave plate are regarded as one object to be measured and the Mueller matrix thereof is described as follows. [Mathematical Expression 48] $\begin{matrix} {{M^{\prime}(\sigma)} = \begin{pmatrix} {\begin{matrix} {{\hat{m}}_{00}^{\prime}(\sigma)} \\ {{\hat{m}}_{10}^{\prime}(\sigma)} \\ {{\hat{m}}_{20}^{\prime}(\sigma)} \end{matrix}\begin{matrix} {{\hat{m}}_{01}^{\prime}(\sigma)} \\ {{\hat{m}}_{11}^{\prime}(\sigma)} \\ {{\hat{m}}_{21}^{\prime}(\sigma)} \end{matrix}\begin{matrix} {{\hat{m}}_{02}^{\prime}(\sigma)} \\ {{\hat{m}}_{12}^{\prime}(\sigma)} \\ {{\hat{m}}_{22}^{\prime}(\sigma)} \end{matrix}\begin{matrix} {{\hat{m}}_{03}^{\prime}(\sigma)} \\ {{\hat{m}}_{13}^{\prime}(\sigma)} \\ {{\hat{m}}_{23}^{\prime}(\sigma)} \end{matrix}} \\ \begin{matrix} {{\hat{m}}_{30}^{\prime}(\sigma)} & {{\hat{m}}_{31}^{\prime}(\sigma)} & {{\hat{m}}_{32}^{\prime}(\sigma)} & {{\hat{m}}_{33}^{\prime}(\sigma)} \end{matrix} \end{pmatrix}} & (6.1) \end{matrix}$

Spectroscopic quasi-Stokes parameters are obtained which are described using Mueller matrix elements in the frame of Expression (6.1). Here, it is assumed that a Mueller matrix of the sample essentially required to be known is as follows. [Mathematical Expression 49] $\begin{matrix} {{M(\sigma)} = \begin{pmatrix} {{\hat{m}}_{00}(\sigma)} & {{\hat{m}}_{01}(\sigma)} & {{\hat{m}}_{02}(\sigma)} & {{\hat{m}}_{03}(\sigma)} \\ {{\hat{m}}_{10}(\sigma)} & {{\hat{m}}_{11}(\sigma)} & {{\hat{m}}_{12}(\sigma)} & {{\hat{m}}_{13}(\sigma)} \\ {{\hat{m}}_{20}(\sigma)} & {{\hat{m}}_{21}(\sigma)} & {{\hat{m}}_{22}(\sigma)} & {{\hat{m}}_{23}(\sigma)} \\ {{\hat{m}}_{30}(\sigma)} & {{\hat{m}}_{31}(\sigma)} & {{\hat{m}}_{32}(\sigma)} & {{\hat{m}}_{33}(\sigma)} \end{pmatrix}} & (6.2) \end{matrix}$ Then, the relation between this Mueller matrix and a Mueller matrix M′(σ) of the object to be measured is as follows. [Mathematical Expression 50] $\begin{matrix} \begin{matrix} {{M(\sigma)} = \begin{pmatrix} {{\hat{m}}_{00}(\sigma)} & {{\hat{m}}_{01}(\sigma)} & {{\hat{m}}_{02}(\sigma)} & {{\hat{m}}_{03}(\sigma)} \\ {{\hat{m}}_{10}(\sigma)} & {{\hat{m}}_{11}(\sigma)} & {{\hat{m}}_{12}(\sigma)} & {{\hat{m}}_{13}(\sigma)} \\ {{\hat{m}}_{20}(\sigma)} & {{\hat{m}}_{21}(\sigma)} & {{\hat{m}}_{22}(\sigma)} & {{\hat{m}}_{23}(\sigma)} \\ {{\hat{m}}_{30}(\sigma)} & {{\hat{m}}_{31}(\sigma)} & {{\hat{m}}_{32}(\sigma)} & {{\hat{m}}_{33}(\sigma)} \end{pmatrix}} \\ {= \begin{pmatrix} {{\hat{m}}_{00}^{\prime}(\sigma)} & {{\hat{m}}_{01}^{\prime}(\sigma)} & {{\hat{m}}_{02}^{\prime}(\sigma)} & {{\hat{m}}_{03}^{\prime}(\sigma)} \\ {{\hat{m}}_{10}^{\prime}(\sigma)} & {{\hat{m}}_{11}^{\prime}(\sigma)} & {{\hat{m}}_{12}^{\prime}(\sigma)} & {{\hat{m}}_{13}^{\prime}(\sigma)} \\ {{\hat{- m}}_{30}^{\prime}(\sigma)} & {{\hat{- m}}_{31}^{\prime}(\sigma)} & {{\hat{- m}}_{32}^{\prime}(\sigma)} & {{\hat{- m}}_{33}^{\prime}(\sigma)} \\ {{\hat{m}}_{20}^{\prime}(\sigma)} & {{\hat{m}}_{21}^{\prime}(\sigma)} & {{\hat{m}}_{22}^{\prime}(\sigma)} & {{\hat{m}}_{23}^{\prime}(\sigma)} \end{pmatrix}} \end{matrix} & (6.3) \end{matrix}$ This means that the spectroscopic quasi-tokes parameters obtained using this polarimeter are determined by the elements on the first row, the second row and the fourth row of Expression (6.2) of the Mueller matrix of the sample which is essentially required to be known. Namely, the elements associated with the spectroscopic quasi-Stokes parameters are different from those in the case without the quarter wave plate. However, the number of parameters obtained is four, which is the same as in the case without the quarter wave plate. It should thus be noted that some of the spectropolarization parameters of the sample can be calculated also in the case with the quarter wave plate. As apparent form the above, the use of the quarter wave plate or the like permits for changing the relational expression of the Mueller matrix of the sample and the spectroscopic quasi-Stokes parameters to be measured. Actively using this fact can for example lead to enhancement of measurement sensitivity with respect to the specific spectropolarization parameter. As an example of such use, the case of calculating the retardation of the sample shown in Section 4.2 is considered. In the case without the quarter wave plate, when the azimuth angle R almost agrees with the azimuth angle θ of the analyzer, sin(2R−2θ) in Expressions (5.2a) to (5.2d) become a value close to zero, the spectroscopic quasi-Stokes parameter has almost no sensitivity with respect to the retardation δ. This results in degradation in measurement accuracy of δ. On the other hand, in the case with the quarter wave plate, the spectroscopic quasi-Stokes parameters are obtained from the first row, the second row and the fourth row in Expression (5.1) of the Mueller matrix of the sample. [Mathematical Expression 51] ${M_{0}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}}$ $\begin{matrix} {{M_{1}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\left\{ {{\left( {{\cos^{2}\quad 2R}\quad + {\cos\quad{\delta(\sigma)}\sin^{2}2\quad R}} \right)\cos\quad 2\quad\theta} -} \right.}} \\ \left. {\sin\quad{\delta(\sigma)}\sin\quad 2R\quad\sin\quad 2\quad\theta} \right\} \end{matrix}$ $\begin{matrix} {{M_{2}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\left\{ {\cos\quad 2R\quad\sin\quad 2R\quad\left( {1 - {\cos\quad{\delta(\sigma)}\cos\quad 2\quad\theta} -} \right.} \right.}} \\ \left. {\sin\quad{\delta(\sigma)}\cos\quad 2R\quad\sin\quad 2\quad\theta} \right\} \end{matrix}$ ${M_{3}(\sigma)} = {\frac{1}{2}{P_{0}(\sigma)}\left\{ {{\sin\quad{\delta(\sigma)}\sin\quad 2R\quad\cos\quad 2\quad\theta} + {\cos\quad\delta\quad(\sigma)\sin\quad 2\quad\theta}} \right\}}$ Hence the sensitivity with respect to the retardation δ can be obtained regardless of the azimuth angle R. It is further found that the sensitivity becomes constant especially when θ=45°.

5. Demodulation of Analyzer and Known Polarization Element and Expansion of Number of Measurable Spectropolarization Parameters by Modulation

As described in the previous chapters, the use of measurement principle of the present invention allows concurrent and independent measurement of four spectroscopic quasi-tokes parameters in one spectrum measurement. This links to the characteristic of being capable of concurrently determining a plurality of spectropolarization parameters regarding an object to be measured (i.e. sample included therein).

However, there are some cases where measurement is not sufficiently performed due to insufficient information obtained only from the above-mentioned four spectroscopic quasi-Stokes parameters, depending upon properties of the sample. Examples of such a case include a case where the polarization parameter necessary to concurrently measure exceeds four.

As shown in Expression (1.1), the Mueller matrix for expressing the SOP of the sample has sixteen elements, and depending upon the sample, each of those sixteen element may be a different value. For example, in the case of measuring the spectropolarization characteristic of the sample by reflecting light on the sample, the condition of the sample was being an isotropic medium in Expression (4.1), but when the sample is an anisotropic medium, the sixteen elements are expressed by up to seven independent parameter equations. This was demonstrated by G. E. Jellison, Jr, (“Handbook of ellipsometry”, edited by H. G Thompkins and E. A. Irene, William Andrew Publishing, P. 244). Further, there is even a case where, when the sample is a heterogeneous medium, the sixteen elements are all independent parameters in measurement of transmitted light or reflected light on or through the sample.

In the present chapter described is expansion of the principle of the present invention for the above-mentioned case where a large number of parameters are required to be obtained. It is to be noted that, although applying this expansion technique causes loss of the characteristic of the channeled spectroscopic polarimetry of “needing no mechanical or active polarization control element”, another advantage of “needing to measure a necessary spectrum only an extremely few times” is generated which is not included in the characteristics of the corresponding conventional method.

5.1 Relation Between Spectroscopic Quasi-Stokes Parameter and Mueller Matrix of Sample

Prior to description on a method for expanding the principle, as preparation for performing the method, a relational expression between the spectroscopic quasi-Stokes parameters and Mueller matrix of the sample is derived. A case is considered here where the object to be measured is composed of the sample D and the known polarization element E, as in FIG. 19. Mueller matrixes of the sample D and the known polarization element E are described as follows. [Mathematical Expression 52] $\begin{matrix} {{M_{sample}(\sigma)} = \begin{bmatrix} {{\hat{m}}_{s\quad 00}(\sigma)} & {{\hat{m}}_{s\quad 01}(\sigma)} & {{\hat{m}}_{s\quad 02}(\sigma)} & {{\hat{m}}_{s\quad 03}(\sigma)} \\ {{\hat{m}}_{s\quad 10}(\sigma)} & {{\hat{m}}_{s\quad 11}(\sigma)} & {{\hat{m}}_{s\quad 12}(\sigma)} & {{\hat{m}}_{s\quad 13}(\sigma)} \\ {{\hat{m}}_{s\quad 20}(\sigma)} & {{\hat{m}}_{s\quad 21}(\sigma)} & {{\hat{m}}_{s\quad 22}(\sigma)} & {{\hat{m}}_{s\quad 23}(\sigma)} \\ {{\hat{m}}_{s\quad 30}(\sigma)} & {{\hat{m}}_{s\quad 31}(\sigma)} & {{\hat{m}}_{s\quad 32}(\sigma)} & {{\hat{m}}_{s\quad 33}(\sigma)} \end{bmatrix}} & \left( {7.1a} \right) \\ {{M_{pc}(\sigma)} = \begin{bmatrix} {{\hat{m}}_{p\quad 00}(\sigma)} & {{\hat{m}}_{p\quad 01}(\sigma)} & {{\hat{m}}_{p\quad 02}(\sigma)} & {{\hat{m}}_{p\quad 03}(\sigma)} \\ {{\hat{m}}_{p\quad 10}(\sigma)} & {{\hat{m}}_{p\quad 11}(\sigma)} & {{\hat{m}}_{p\quad 12}(\sigma)} & {{\hat{m}}_{p\quad 13}(\sigma)} \\ {{\hat{m}}_{p\quad 20}(\sigma)} & {{\hat{m}}_{p\quad 21}(\sigma)} & {{\hat{m}}_{p\quad 22}(\sigma)} & {{\hat{m}}_{p\quad 23}(\sigma)} \\ {{\hat{m}}_{p\quad 30}(\sigma)} & {{\hat{m}}_{p\quad 31}(\sigma)} & {{\hat{m}}_{p\quad 32}(\sigma)} & {{\hat{m}}_{p\quad 33}(\sigma)} \end{bmatrix}} & \left( {7.1b} \right) \end{matrix}$ The Mueller matrix of the object to be measured is given by the product of the above two Mueller matrixes, i.e. M(σ)=M _(pc)(u)M _(sample)(σ)  (7.2) Expression (7.2) is written as follows using the element of each matrix. $\begin{matrix} {{{\hat{m}}_{ij}(\sigma)} = {\sum\limits_{k = 0}^{3}{{{\hat{m}}_{pik}(\sigma)}{{\hat{m}}_{skj}(\sigma)}}}} & (7.3) \end{matrix}$ where i and j are integers from 0 to 3.

When Expression (7.3) is substituted into Expressions (1.6a) to (1.6d), the following expression is derived for relating the spectroscopic quasi-Stokes parameter M₁(σ) of the object to be measured (where I=0 . . . 3) to the Mueller matrix of the sample. [Mathematical Expression 53] $\begin{matrix} \begin{matrix} {{M_{l}(\sigma)} = {\frac{1}{2}{{P_{0}(\sigma)}\begin{bmatrix} {{\sum\limits_{k = 0}^{3}{{{\hat{m}}_{p0k}(\sigma)}{{\hat{m}}_{skl}(\sigma)}}} +} \\ {\sum\limits_{k = 0}^{3}{{\hat{m}}_{plk}\left( \sigma\quad \middle| \quad{{{{\hat{m}}_{skl}(\sigma)}\cos\quad 2\quad\theta} +} \right.}} \\ {\sum\limits_{k = 0}^{3}{{{\hat{m}}_{p2k}(\sigma)}{{\hat{m}}_{skl}(\sigma)}\sin\quad 2\quad\theta}} \end{bmatrix}}}} \\ {= {{P_{0}(\sigma)}{\sum\limits_{k = 0}^{3}{{{\hat{a}}_{k}(\sigma)}{{\hat{m}}_{skl}(\sigma)}}}}} \end{matrix} & (7.4) \end{matrix}$ where $\begin{matrix} {{{\hat{a}}_{k}(\sigma)} = {{\frac{1}{2}{{\hat{m}}_{p0k}(\sigma)}} + {\frac{1}{2}{{\hat{m}}_{plk}(\sigma)}\cos\quad 2\quad\theta} + {\frac{1}{2}{{\hat{m}}_{p2k}(\sigma)}\sin\quad 2\quad\theta}}} & (7.5) \end{matrix}$ The next relational expression using vectors and columns is immediately derived from Expression (7.4). $\begin{matrix} \begin{matrix} {\begin{bmatrix} {M_{0}(\sigma)} & {M_{1}(\sigma)} & {M_{2}(\sigma)} & {M_{3}(\sigma)} \end{bmatrix} = \begin{matrix} {{P_{0}(\sigma)}\begin{bmatrix} {{\hat{a}}_{0}(\sigma)} & {{\hat{a}}_{1}(\sigma)} & {{\hat{a}}_{2}(\sigma)} & {{\hat{a}}_{3}(\sigma)} \end{bmatrix}} \\ \begin{bmatrix} {{\hat{m}}_{s\quad 00}(\sigma)} & {{\hat{m}}_{s\quad 01}(\sigma)} & {{\hat{m}}_{s\quad 02}(\sigma)} & {{\hat{m}}_{s\quad 03}(\sigma)} \\ {{\hat{m}}_{s\quad 10}(\sigma)} & {{\hat{m}}_{s\quad 11}(\sigma)} & {{\hat{m}}_{s\quad 12}(\sigma)} & {{\hat{m}}_{s\quad 13}(\sigma)} \\ {{\hat{m}}_{s\quad 20}(\sigma)} & {{\hat{m}}_{s\quad 21}(\sigma)} & {{\hat{m}}_{s\quad 22}(\sigma)} & {{\hat{m}}_{s\quad 23}(\sigma)} \\ {{\hat{m}}_{s\quad 30}(\sigma)} & {{\hat{m}}_{s\quad 31}(\sigma)} & {{\hat{m}}_{s\quad 32}(\sigma)} & {{\hat{m}}_{s\quad 33}(\sigma)} \end{bmatrix} \end{matrix}} \\ {= {{{P_{0}(\sigma)}\begin{bmatrix} {{\hat{a}}_{0}(\sigma)} & {{\hat{a}}_{1}(\sigma)} & {{\hat{a}}_{2}(\sigma)} & {{\hat{a}}_{3}(\sigma)} \end{bmatrix}}{M_{sample}(\sigma)}}} \end{matrix} & (76) \end{matrix}$ The above expression relates the four spectroscopic quasi-Stokes parameters M₀(σ), M₁(σ), M₂(σ), and M₃(σ) directly to the Mueller matrix M_(sample)(σ) of the sample. Here, the element â_(k)(σ)(k=0 . . . 3) of the vector for relating the spectroscopic quasi-Stokes parameters to the Mueller matrix is an amount determined only by the characteristic of the known polarization element E and the azimuth angle θ of the analyzer A (azimuth angle of the transmission axis of the analyzer A with respect to the fast axis of the retarder R1), as apparent from the definition of Expression (7.5), and does not depend upon the “Mueller matrix of the sample”. It is to be noted that the above argument can be applied as it is to the case without the known polarization element E after the sample D in FIG. 19 if Expression (7.1) is replaced by a 4×4 unit matrix.

5.2 (Method for Increasing Number of Measurable Parameters (and Method for Completely Measuring Mueller Matrix)

As described in the previous chapters, it is possible to concurrently and independently measure the spectroscopic quasi-Stokes parameters M₀(σ), M₁(σ), M₂(σ), and M₃(σ) by the method of the present invention. By forming an equation based upon Expression (7.6) using the obtained spectroscopic quasi-Stokes parameters, it is possible to obtain up to four spectroscopic quasi-Stokes parameters regarding the sample, as shown in the example of Chapter 4. [{circumflex over (α)}₀(σ) {circumflex over (α)}₁(σ) {circumflex over (α)}₂(σ) {circumflex over (α)}₃(σ)]  [Mathematical Expression 54]

However, there could be cases where a large number of spectropolarization parameters need to be obtained, or the above equation is difficult to solve. In such cases, the number of equations can be increased by changing a vector for giving a coefficient in several times and repeating measurement at each of the changes. Since this coefficient vector depends only upon the characteristic of the known polarization element E and the azimuth angle θ of the analyzer A, controlling either one of them allows changing of the vector. With the number of equations increased, it is possible to increase the number of independent spectropolarization parameters which are concurrently obtained.

Here, especially the case of performing measurements under four different conditions regarding the known polarization element E and the analyzer A is considered. When spectroscopic quasi-Stokes parameters obtained in the respective cases are differentiated by means of superscript characters (p)=0 . . . 3, the following relational expressions of matrixes are established, as seen from Expression (7.6). [Mathematical Expression 55] $\begin{matrix} {\begin{bmatrix} {M_{0}^{(0)}(\sigma)} & {M_{1}^{(0)}(\sigma)} & {M_{2}^{(0)}(\sigma)} & {M_{3}^{(0)}(\sigma)} \\ {M_{0}^{(1)}(\sigma)} & {M_{1}^{(1)}(\sigma)} & {M_{2}^{(1)}(\sigma)} & {M_{3}^{(1)}(\sigma)} \\ {M_{0}^{(2)}(\sigma)} & {M_{1}^{(2)}(\sigma)} & {M_{2}^{(2)}(\sigma)} & {M_{3}^{(2)}(\sigma)} \\ {M_{0}^{(3)}(\sigma)} & {M_{1}^{(3)}(\sigma)} & {M_{2}^{(3)}(\sigma)} & {M_{3}^{(3)}(\sigma)} \end{bmatrix} = {{{P_{0}(\sigma)}\begin{bmatrix} {{\hat{a}}_{0}^{(0)}(\sigma)} & {{\hat{a}}_{1}^{(0)}(\sigma)} & {{\hat{a}}_{2}^{(0)}(\sigma)} & {{\hat{a}}_{3}^{(0)}(\sigma)} \\ {{\hat{a}}_{0}^{(1)}(\sigma)} & {{\hat{a}}_{1}^{(1)}(\sigma)} & {{\hat{a}}_{2}^{(1)}(\sigma)} & {{\hat{a}}_{3}^{(1)}(\sigma)} \\ {{\hat{a}}_{0}^{(2)}(\sigma)} & {{\hat{a}}_{1}^{(2)}(\sigma)} & {{\hat{a}}_{2}^{(2)}(\sigma)} & {{\hat{a}}_{3}^{(2)}(\sigma)} \\ {{\hat{a}}_{0}^{(3)}(\sigma)} & {{\hat{a}}_{1}^{(3)}(\sigma)} & {{\hat{a}}_{2}^{(3)}(\sigma)} & {{\hat{a}}_{3}^{(3)}(\sigma)} \end{bmatrix}}{{M_{sample}(\sigma)}.}}} & (7.7) \end{matrix}$ In this expression, the first matrix on the right side above is: $\begin{matrix} {{N(\sigma)} = \begin{bmatrix} {{\hat{a}}_{0}^{(0)}(\sigma)} & {{\hat{a}}_{1}^{(0)}(\sigma)} & {{\hat{a}}_{2}^{(0)}(\sigma)} & {{\hat{a}}_{3}^{(0)}(\sigma)} \\ {{\hat{a}}_{0}^{(1)}(\sigma)} & {{\hat{a}}_{1}^{(1)}(\sigma)} & {{\hat{a}}_{2}^{(1)}(\sigma)} & {{\hat{a}}_{3}^{(1)}(\sigma)} \\ {{\hat{a}}_{0}^{(2)}(\sigma)} & {{\hat{a}}_{1}^{(2)}(\sigma)} & {{\hat{a}}_{2}^{(2)}(\sigma)} & {{\hat{a}}_{3}^{(2)}(\sigma)} \\ {{\hat{a}}_{0}^{(3)}(\sigma)} & {{\hat{a}}_{1}^{(3)}(\sigma)} & {{\hat{a}}_{2}^{(3)}(\sigma)} & {{\hat{a}}_{3}^{(3)}(\sigma)} \end{bmatrix}} & (7.8) \end{matrix}$ When this matrix has an inverse matrix N⁻¹(σ), it becomes possible to inverse transform Expression (7.7), so as to obtain the Mueller matrix M_(sample)(σ) as: $\begin{matrix} {{M_{sample}(\sigma)} = {\frac{1}{P_{0}(\sigma)}{{{N^{- 1}(\sigma)}\begin{bmatrix} {M_{0}^{(0)}(\sigma)} & {M_{1}^{(0)}(\sigma)} & {M_{2}^{(0)}(\sigma)} & {M_{3}^{(0)}(\sigma)} \\ {M_{0}^{(1)}(\sigma)} & {M_{1}^{(1)}(\sigma)} & {M_{2}^{(1)}(\sigma)} & {M_{3}^{(1)}(\sigma)} \\ {M_{0}^{(2)}(\sigma)} & {M_{1}^{(2)}(\sigma)} & {M_{2}^{(2)}(\sigma)} & {M_{3}^{(2)}(\sigma)} \\ {M_{0}^{(3)}(\sigma)} & {M_{1}^{(3)}(\sigma)} & {M_{2}^{(3)}(\sigma)} & {M_{3}^{(3)}(\sigma)} \end{bmatrix}}.}}} & (7.9) \end{matrix}$ This expression means that, by changing either one or both of the characteristic of the known polarization element E and the azimuth angle θ of the analyzer A to perform measurements of at least four kinds of spectroscopic quasi-Stokes parameter, all of the sixteen elements of the Mueller matrix of the sample can be concurrently and independently determined. It is however necessary at this time to control the known polarization element E and the analyzer A so that N(σ) has an inverse matrix. It is to be noted that, if the number of conditions for measurement increases more than four, the number of equations can be increased. If the idea of least-squares is employed, an error due to an influence of noise or the like can be reduced. On the contrary, when the number of measurements is less than four, all the sixteen elements cannot be independently determined, but still a large number of spectroscopic quasi-Stokes parameter expressions can be obtained as compared with the case of not changing the characteristic of the known polarization element E and the azimuth angle θ of the analyzer A, resulting in that a larger number of spectropolarization parameters of the sample can be obtained.

Here, attention should be given to that the present measurement method has a large advantage as compared with the conventional method even in the case of moving the known polarization element E and the analyzer A. Although “mechanical or active polarization control” is necessary for certain, the number of necessary measurement steps is significantly different from that of the conventional method. According to the present method, four spectroscopic quasi-Stokes parameters are obtained in one measurement, meaning that in terms of the number of parameters to be obtained, the number of measurements can be reduced to the order of a quarter of the number of measurements in the conventional method. For example, when all of the 16 elements of the Mueller matrix are to be measured, only four times of measurements are required at the minimum in the present method, whereas in the conventional method, the spectrum measurement needs to be repeated sixteen times at the minimum, and usually 20 to 30 times. This means that the present method has large advantages in terms of reduction in measurement time or simplification of the measurement system.

In addition, the characteristic of the known polarization element E and the azimuth angle θ of the analyzer A can be changed by a variety of methods. A first method is changing either or both of azimuth angles θ of the two elements. Such a change may be made by actually rotating the element, or replacing the element by an element having a different azimuth angle θ, or inserting a Faraday cell or the like in a position before the element so as to magneto-optically rotate the azimuth angle θ of the element from the installation orientation. A second method is introducing a compensator capable of moduclation such as an electro-optic modulator, a photoelastic modulator, or a liquid-crystal optical modulator, to change a retardation of one of parameters for determining a Mueller matrix of this element. A third method is combining the above-mentioned methods. (It is to be noted that the above-mentioned methods do not limit the method to be performed.) Further, the known polarization element E is not necessarily composed of a single element. For example, the known polarization element E may be constituted by combination of a plurality of compensators capable of modulation.

5.3 Example

The case of rotating a compensator (lower-order retarder or zero-order retarder) as the known polarization element E is shown. When the retardation of this compensator is δ_(c)(σ), and the azimuth thereof is θ_(c),

[Mathematical Expression 56]

â_(k)(σ)(k=0 . . . 3) are expressed as follows. $\begin{matrix} {{{\hat{a}}_{0}(\sigma)} = \frac{1}{2}} & \left( {7.10a} \right) \\ \begin{matrix} {{{\hat{a}}_{1}(\sigma)} = {{\frac{1}{2}\cos\quad 2\theta_{C}\cos\quad 2\left( {\theta - \theta_{C}} \right)} +}} \\ {\frac{1}{2}\cos\quad{\delta_{C}(\sigma)}\sin\quad 2\quad\theta_{C}\sin\quad 2\left( {\theta - \theta_{C}} \right)} \end{matrix} & \left( {7.10b} \right) \\ \begin{matrix} {{{\hat{a}}_{2}(\sigma)} = {{\frac{1}{2}\sin\quad 2\theta_{C}\cos\quad 2\left( {\theta - \theta_{C}} \right)} +}} \\ {{- \frac{1}{2}}\cos\quad{\delta_{C}(\sigma)}\cos\quad 2\quad\theta_{C}\sin\quad 2\left( {\theta - \theta_{C}} \right)} \end{matrix} & \left( {7.10c} \right) \\ {{{\hat{a}}_{3}(\sigma)} = {{- \frac{1}{2}}\sin\quad{\delta_{C}(\sigma)}\sin\quad 2\left( {\theta - \theta_{C}} \right)}} & \left( {7.10d} \right) \end{matrix}$ As apparent from this expression, changing the azimuth angle θ of the analyzer or the azimuth angle θ_(c) of the compensator enables control of the following vectors which give coefficients. [{circumflex over (α)}₀(σ) {circumflex over (α)}₁(σ) {circumflex over (α)}₂(σ) {circumflex over (α)}₃(σ)] That is, channeled spectrums may be repeatedly measured while either one of the azimuth angles of the two elements is changed.

For example, when θ_(c) is changed to four different angles: −45°, 0°, 30°, and 60°, on the basis of δ_(c)=90° and θ=45°, a matrix N(σ) given by Expression (7.8) is expressed as follows. [Mathematical Expression 57] $\begin{matrix} {N = {\frac{1}{2}\begin{bmatrix} 1 & 0 & 1 & 0 \\ 1 & 0 & 0 & 1 \\ 1 & \frac{\sqrt{3}}{4} & \frac{3}{4} & \frac{1}{2} \\ 1 & {- \frac{\sqrt{3}}{4}} & \frac{3}{4} & {- \frac{1}{2}} \end{bmatrix}}} & (7.11) \end{matrix}$ Here, the inverse matrix of this matrix is given as follows. $\begin{matrix} {N^{- 1} = {2\begin{bmatrix} {- 1} & 0 & 1 & 1 \\ {- \frac{4}{\sqrt{3}}} & {- \frac{2}{\sqrt{3}}} & \frac{5}{\sqrt{3}} & \frac{1}{\sqrt{3}} \\ 2 & 0 & {- 1} & {- 1} \\ 1 & 1 & {- 1} & {- 1} \end{bmatrix}}} & (7.12) \end{matrix}$ By substituting spectroscopic quasi-tokes parameters obtained by four times of measurements and N⁻¹ mentioned above into Expression (7.9), it is possible to determine all the sixteen elements of the Mueller matrix of the sample for each wavenumber.

It should be noted that in actual measurement, a retardation δ_(c)(σ) of a compensation is a function of a wavenumber a and is not constant. However, even in such a case, the measurement remains unsusceptible since Expression (7.9) is calculated for each wavenumber. Further, according to demonstration similar to that in the case of “spectroscopic polarization state measurement by retarder rotating method” shown in a reference document (“Polarized light” written by D. Goldstein, Mercel Dekker Inc., p. 555), it is possible to demonstrate that 132° is the optimum value as the retardation δ_(c)(σ). The closer to this value the retardation is, the less susceptible to measurement noise the measurement can be.

In the meantime, the number of equations can be increased also by rotating the azimuth angle θ of the analyzer. However, it is necessary to note that rotating only the analyzer prevents the matrix N(σ) from having an inverse matrix due to the property of the analyzer. If all the sixteen Mueller matrix elements are to be obtained, it is at least necessary to rotate the compensator.

Chapter 6 Channeled Spectroscopic Polarization State Generator

As the embodiment of the present invention, it was described in Chapter 1 that the optical system is comprised of the light source 7, the polarizer P, the retarders R2 and R1, the analyzer A and the spectroscope 8, and spectropolarization parameters of the sample and the like are calculated by analysis of a spectrum of incident light acquired in the spectroscope 8 in the foregoing procedure. Meanwhile, when the role of the light-projection part (the light source 7, the polarizer P, and the retarders R2 and R1) of the optical system is considered, this part can be defined as a “spectroscopic polarization state generator” for emitting light having a modulated SOP. This is especially named a “channeled spectroscopic polarization state generator (hereinafter referred to as CSPSG)”. In this chapter, the optical implication of this generator is described.

FIG. 30 shows the configuration of the channeled spectroscopic polarization state generator (CSPSG). This optical system is configured to allow light emitted from the light source 7 to transit through the polarizer P and the retarders R1 and R2. The constituent elements in this configuration are the same as those starting with the light source and ending with the component before the sample in FIG. 2. Further, the azimuth angle of the element is the same as in FIG. 2. At this time, light emitted from the CSPSG is light having an SOP modulated along the wavenumber axis. A stokes spectrum S_(PSG)(σ) emitted from the CSPSG is expressed below by calculation using the Mueller matrix. [Mathematical Expression 58] $\begin{matrix} {{S_{PSG}(\sigma)} = {\begin{pmatrix} {S_{0}(\sigma)} \\ {S_{1}(\sigma)} \\ {S_{2}(\sigma)} \\ {S_{3}(\sigma)} \end{pmatrix} = {\frac{1}{2}{P_{o}(\sigma)}\begin{pmatrix} 1 \\ {\cos\left( {\phi_{1}(\sigma)} \right)} \\ {\sin\quad{\phi_{1}(\sigma)}\sin\quad{\phi_{2}(\sigma)}} \\ {\sin\quad{\phi_{1}(\sigma)}\cos\quad{\phi_{2}(\sigma)}} \end{pmatrix}}}} & (8.1) \end{matrix}$ Here, φ_(j)(σ) (j=1, 2) is the retardation of the retarder formulated by Expression (1.2). Assuming that dispersion of a birefringence B(σ) of the retarder is not so large, φ_(j)(σ) increases almost linearly with respect to wave number a, as seen from Expression (1.2). By substitution of Expression (1.2) into Expression (8.1), the following expression is formed. [Mathematical Expression 59] $\begin{matrix} \begin{matrix} {{S_{PSG}(\sigma)} = \begin{pmatrix} {S_{0}(\sigma)} \\ {S_{1}(\sigma)} \\ {S_{2}(\sigma)} \\ {S_{3}(\sigma)} \end{pmatrix}} \\ {= {\frac{1}{2}{P_{o}(\sigma)}\begin{pmatrix} \begin{matrix} \begin{matrix} 1 \\ {\cos\left\{ {{2\pi\quad L_{1}\sigma} + {\Phi_{1}(\sigma)}} \right\}} \end{matrix} \\ {\frac{1}{2}\left\lbrack {{\cos\left\{ {{2\pi\quad L_{-}\sigma} + {\Phi_{-}(\sigma)}} \right\}} - {\cos\left\{ {{2\pi\quad L_{+}\sigma} + {\Phi_{+}(\sigma)}} \right\}}} \right\rbrack} \end{matrix} \\ {\frac{1}{2}\left\lbrack {{\sin\left\{ {{2\pi\quad L_{-}\sigma} + {\Phi_{-}(\sigma)}} \right\}} - {\sin\left\{ {{2\pi\quad L_{+}\sigma} + {\Phi_{+}(\sigma)}} \right\}}} \right\rbrack} \end{pmatrix}}} \end{matrix} & (8.2) \end{matrix}$ It is found from Expression (8.2) that the light emitted from the CSPSG is light having been modulated along a wavenumber axis having three characteristics as follows. (a) S₁(σ) is modulated in a quasi-sinusoidal manner at a period 1/L₁. (b) S₂(σ) and S₃(σ) are both composed of two components modulated in a quasi-sinusoidal manner at a period 1/L and a period 1/L. In (c) S₂(σ) and S₃(σ), the quasi-sinusoidal components at the same period has initial phases 90° different from each other. Therefore, the light emitted from the CSPSG can be considered as light having four spectroscopic Stokes parameters modulated at a period or phase independently different from one another. It can thereby be said that this CSPSG is a complete spectroscopic polarization state generator. The present invention can be defined as having a configuration for obtaining a spectropolarization parameter of an object to be measured, formed by combination, of the light source of the above-mentioned CSPSG as a complete spectroscopic polarization state generator with the light source, the analyzer and the spectrometer.

Example 1

In the following, a preferred example of the present invention is specifically described with reference to FIGS. 20 to 23. FIG. 20 shows a configuration view of one example of a spectroscopic polarimeter. As shown in this figure, this device comprises a light-projection side unit 200 and a light-reception side unit 300. It is to be noted that reference numeral 400 denotes a sample.

The light-projection side unit 200 comprises: a power source 201; a light source 202 that is turned on by power feeding from the power source 201; a pinhole plate 203 arranged on the front face side of the light source 202 in the light emitting direction; a collimator lens 204 for collimating light transmitting through the pinhole of the 203; a shutter 205 which is arranged on the front face side of the collimator lens 204 and opens and closes to transmit or block the transmitted light; a polarizer 206 on which the light having transmitted through the shutter is incident; and a second retarder 207 and a first retarder 208, through which the light having transmitted through the polarizer transmits in this order.

The light after passage of the first retarder 208 is emitted from the light-projection side unit 200 and applied to the sample 400. The light transmitted through or reflected on the sample 400 is incident on the light-reception side unit 300.

On an incident light channel in the light-reception side unit 300, an analyzer 301, and a spectrometer 302 intervene in this order. Here, a respective angle between the first retarder 208 and the analyzer 301 is set to be a known angle.

The spectrometer 302 comprises: a diffraction grating 302 a for spatially dispersing the incident light; a CCD 302 b with the light-reception face on which light spatially dispersed by the diffraction grating 302 a is incident; and an A/D converter 302 c for converting light-reception output from the CCD 302 b into a digital signal. The digital light-reception output signal obtained from the A/D converter 302 c is taken out from the spectrometer 302, and then processed in a computer 303 such as a personal computer (PC).

As widely known, the computer 303 comprises: an arithmetic processing part 303 a comprised of a microprocessor and the like; a memory part 303 b comprised of an ROM, an RAM, an HDD and the like; and a measurement result output part 303 c comprised of a display, a printer, a variety of data output devices, a communication device, and the like.

Next, FIG. 21 shows a more specific configuration view regarding a sensor head part of the spectroscopic polarimeter. A sensor head part 100 comprises: a light-projection part 110 for emitting light; a light-reception part 120 for receiving light having been reflected on or transmitted through a sample; and a housing 130 for protecting the light-projection part 110 and the light-reception part 120. It is to be noted that reference numeral 50 denotes a sample.

The light-projection part 110 includes: a fiber-optic cable 111 for allowing light emitted from a light source (not shown) to transmit therethrough; a cable head 112 for allowing transmitted light from the fiber-optic cable 111 to transmit therethrough; a collimator lens (light-projection lens) 115 for collimating transmitted light from the cable head 112; a polarizer 116 which is arranged on the front face side of the collimator lens 115 and allows incident light to transmit therethrough; a second retarder 117 and a first retarder 118, through which light emitted from the polarizer transmits in this order; and an optical system holding member 113 and a fixing member 114 which are used to install the optical system in the housing 130. It should be noted that a solid line 119 is a light-projection axis of the light that transmits through the light-projection part 110.

The light-reception part 120 includes: an analyzer 122 for allowing light reflected on or transmitted through the sample 50 to transmit therethrough; a light-reception lens 123 for condensing transmitted light from the analyzer 122; a cable head 126 for allowing light having transmitted through the light-reception lens 123 to transmit therethrough; a fiber-optic cable 127 connected to a spectrometer (not shown); and a fixing member 124 and an optical system holding member 125 which are used to install the optical system in the housing 130. It should be noted that a solid line 121 is a light-reception axis of the light that is reflected on or transmits through the sample 50.

Next, FIG. 29 shows a device configuration view in the case of installing a polarization element having a known polarization parameter in a position after the sample and measuring the polarization parameter of the sample. As compared with FIG. 21, FIG. 29 is configured to install a compensator (known polarization element) 140 held in a hollow motor 141 between the sample 50 and the analyzer 122. Here, reference numeral 142 denotes electrical wiring for motor drive. A spectroscopic quasi-Stokes parameter can be measured under a plurality of conditions by rotating the hollow motor 141 to control the azimuth angle of the compensator 140. Further, the hollow motor 141 is fixed to the housing 130 to form an integral configuration as an element of the light-reception part 120. The hollow motor 141 is controlled by the arithmetic processing part 303 a in FIG. 20. It is to be noted that in the case of not rotating the compensator 140 but rotating the analyzer 122, the hollow motor 141 may be replaced by a fixing hardware of the compensator 140, to provide a hollow motor for rotating the analyzer 122. Further, both the compensator 140 and the analyzer 122 may be made individually rotatable.

Next, FIG. 22 shows a flowchart of a pre-calibration procedure. As shown in this figure, the precalibration procedure is started with application of light to the device in Step 2201. However, in this device, the relative angle between the first retarder 208 and the analyzer 301 is a known angle, and an element for changing an SOP of light is not arranged between the first retarder 208 and the analyzer 301.

Next, in Step 2202, a spectral intensity of transmitted light from the analyzer 301 is measured using the spectrometer. Here, the shutter 205 may be utilized for reduction in influence of unnecessary light, such as lost light. Specifically, a spectrum of the unnecessary light can be canceled out by taking a difference in spectrum between when measured with the shutter open and when measured with the shutter closed.

Next, in Step 2203, the spectral intensity of the transmitted light received is forwarded from the spectrometer to the computer 303, to be provided to calculation in the arithmetic processing part 303 a.

Next, in Step 2204, reference phase functions and reference amplitude functions are calculated by the action of the arithmetic processing part 303 a.

Next, in Step 2205, the calculated reference phase functions and reference amplitude functions are stored in the memory part 303 b, whereby the pre-calibration procedure is completed.

FIG. 23 shows a flowchart of a measurement procedure. As shown in the figure, the measurement procedure is started with application of light to the device in Step 2301.

Next, in Step 2302, light is reflected on or transmitted through the sample 400 using the spectrometer 302, and thereafter, a spectral intensity of transmitted light having transmitted through the analyzer 301 is measured. Here, the shutter 205 can be utilized for reduction in influence of unnecessary light, such as lost light. Specifically, the spectrum of the unnecessary light can be canceled out by taking a difference in spectrum between when measured with the shutter open and when measured with the shutter closed.

Next, in Step 2303, the spectral intensity of the transmitted light is forwarded from the spectrometer 302 to the computer 303, to be provided to processing in the arithmetic processing part 303 a. At this time, in the case of implementing the procedure described in Chapter 5, the optical system described by use of FIG. 29 is used and the azimuth angle of the compensator 140 or the analyzer 122 is changed to acquire a spectral intensity more than once.

Next, in Step 2304, in the computer 303, the arithmetic processing part 303 a acquires reference phase functions and reference amplitude functions from the memory part 303 b.

Next, in Step 2305, in the computer 303, the arithmetic processing part 303 a calculates reference phase function variations (Δφ₁ and Δφ₂) by use of the measured spectral intensity, the reference phase functions and the reference amplitude functions.

Next, in Step 2306, in the computer 303, the arithmetic processing part 303 a calculates spectrometric Stokes parameters by use of the measured spectral intensity, the reference phase functions, the reference amplitude functions and the reference phase function variations.

Next, in Step 2307, in the computer 303, the arithmetic processing part 303 a outputs the spectrometric Stokes parameters of the sample 400. Examples of the measurement result output part 303 c may include a memory, a hard disc, and other processing parts (calculating part for film thickness, complex refractive index, etc.).

As described above, in the spectroscopic polarimetry of the present example, spectrometric Stokes parameters of the sample are calculated through the pre-calibration procedure shown in FIG. 22 and the measurement procedure shown in FIG. 23 in the system constitution shown in FIGS. 20, 21 and 29. 

1. A spectroscopic polarimetry, comprising the steps of: preparing an object to be measured; preparing a polarimetric spectroscope which includes a projection optical system, comprising a light source, a polarizer and a plurality of retarders, where the light source, the polarizer and the plurality of retarders are arranged such that light emitted from the light source is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order, an analyzer for allowing light to transmit therethrough, the light having been emitted from the projection optical system and reflected on or transmitted through the object to be measured, and a device for obtaining the spectral intensity of the light having transmitted through the analyzer; and obtaining the spectral intensity of the object to be measured by use of the polarimetric spectroscope.
 2. The spectroscopic polarimetry according to claim 1, further comprising a step of obtaining at least one of spectropolarization parameters of the object to be measured by use of the obtained spectral intensity.
 3. The spectroscopic polarimetry according to claim 2, wherein the plurality of retarders that the projection optical system comprises are a first retarder and a second retarder, and the light source, the polarizer, the first retarder and the second retarder are arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder.
 4. The spectroscopic polarimetry according to claim 3, wherein the step of obtaining at least one of spectropolarization parameters is a step which comprises: obtaining, from the spectral intensity, a spectral intensity component (first spectral intensity component) which nonperiodically vibrates with wavenumber and a spectral intensity component (third spectral intensity component) which vibrates at a frequency depending upon a retardation (φ₂(σ)) of the second retarder and not depending upon a retardation (φ₁(σ)) of the first retarder, with wavenumber; and obtaining at least one of spectropolarization parameters by use of each of the spectral intensity components.
 5. The spectroscopic polarimetry according to claim 3, wherein the step of obtaining at least one of spectropolarization parameters is a step which comprises: obtaining, from the spectral intensity, at least one of a spectral intensity component (second spectral intensity component) which vibrates at a frequency depending upon a variation between the retardation (φ₁ (σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, a spectral intensity component (fourth spectral intensity component) which vibrates at a frequency depending upon the sum of the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, and a spectral intensity component (fifth spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₁(σ)) of the first retarder and not depending upon the retardation (φ₂(σ)) of the second retarder, with wavenumber; and obtaining at least one of the spectropolarization parameters of the object to be measured by use of the obtained spectral intensity component.
 6. The spectroscopic polarimetry according to claim 3, wherein the step of obtaining at least one of spectropolarization parameters is a step which comprises: obtaining, from the spectral intensity, at least one of the spectral intensity component (first spectral intensity component) which nonperiodically vibrates with wavenumber and the spectral intensity component (third spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₂(σ)) of the second retarder and not depending upon the retardation (φ₁(σ)) of the first retarder, with wavenumber, and at least one of the spectral intensity component (second spectral intensity component) which vibrates at a frequency depending upon the difference between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, the spectral intensity component (fourth spectral intensity component) which vibrates at a frequency depending upon the sum of the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, and the spectral intensity component (fifth spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₁(σ)) of the first retarder and not depending upon the retardation (φ₂(σ)) of the second retarder, with wavenumber; and obtaining at least one of the spectropolarization parameters of the object to be measured by use of each of the obtained spectral intensity components.
 7. The spectroscopic polarimetry according to claim 3, wherein the step of obtaining at least one of spectropolarization parameters is a step which comprises: obtaining the retardation (φ₂(σ)) of the second retarder from the spectral intensity; and obtaining at least one of the spectropolarization parameters of the object to be measured by use of the spectral intensity and the retardation (φ₂(σ)) of the second retarder.
 8. The spectroscopic polarimetry according to claim 3, further comprising a step of acquiring data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder, wherein the step of obtaining at least one of spectropolarization parameters is a step which comprises: obtaining the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder from the spectral intensity and the data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder; and obtaining at least one of the spectropolarization parameters of the object to be measured by use of the spectral intensity, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.
 9. The spectroscopic polarimetry according to claim 3, further comprising the steps of: acquiring data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder; and acquiring a reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and a reference value (φ₂ ^((i))(σ)) for calibration of retardation of the second retarder, wherein the step of obtaining at least one of spectropolarization parameters is a step which comprises: obtaining, from the spectral intensity, the retardation (φ₂(σ)) of the second retarder and the retardation variation (Δφ₂(σ)) of the second retarder from the reference value (φ₂ ^((i))(σ)) for calibration; obtaining the retardation variation (Δφ₁(σ)) of the first retarder by use of the obtained retardation variation (Δφ₂(σ)) of the second retarder and data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder; obtaining the retardation (φ₁(σ)) of the first retarder from a reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and the obtained retardation variation (Δφ₁(σ)) of the first retarder; and obtaining at least one of the spectropolarization parameters of the object to be measured by use of the spectral intensity, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.
 10. The spectroscopic polarimetry according to claim 3, wherein the polarizer and the second retarder are arranged such that an angle between the orientation of the transmission axis of the polarizer and the orientation of a fast axis of the second retarder is 45°.
 11. The spectroscopic polarimetry according to claim 2, further comprising a step of obtaining a spectral intensity for calibration by use of the polarimetric spectroscope in a state where an object to be measured having an unknown spectropolarization characteristic does not exist in a light path between the projection optical system and the analyzer, wherein the step of obtaining at least one of spectropolarization parameters is a step of obtaining at least one of the spectropolarization parameters of the object to be measured by use of the spectral intensity regarding the object to be measured and the spectral intensity for calibration or data based upon the spectral intensity for calibration.
 12. The spectroscopic polarimetry according to claim 11, wherein the step of obtaining the spectral intensity for calibration is a step of preparing an analyzer for calibration in a position in which light emitted from the projection optical system is received in a state where the object to be measured having an unknown spectropolarization characteristic does not exist in the light path between the projection optical system and the analyzer.
 13. The spectroscopic polarimetry according to claim 11, further comprising a step of obtaining the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder by use of the spectral intensity for calibration, wherein the step of obtaining at least one of spectropolarization parameters is a step of obtaining at least one of the spectropolarization parameters of the object to be measured by use of the spectral intensity regarding the object to be measured, the retardation (φ₁(σ)) of the first retarder, and the retardation (φ₂(σ)) of the second retarder, which are obtained by use of the spectral intensity for calibration.
 14. The spectroscopic polarimetry according to claim 8, wherein the step of acquiring data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder is a step which comprises: obtaining the spectral intensity for calibration by use of the polarimetric spectroscope in a state where the object to be measured having an unknown spectropolarization characteristic does not exist in the light path between the projection optical system and the analyzer; and obtaining the data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder by use of the obtained spectral intensity for calibration.
 15. The spectroscopic polarimetry according to claim 9, wherein the step of acquiring data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder is a step which comprises: obtaining the spectral intensity for calibration by use of the polarimetric spectroscope in a state where the object to be measured having an unknown spectropolarization characteristic does not exist in the light path between the projection optical system and the analyzer; and obtaining the data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder by use of the obtained spectral intensity for calibration.
 16. The spectroscopic polarimetry according to claim 1, comprising a step of obtaining the spectroscopic quasi-tokes parameter of the object to be measured by use of the obtained spectral intensity.
 17. The spectroscopic polarimetry according to claim 16, wherein the plurality of retarders that the projection optical system comprises are a first retarder and a second retarder, the light source, the polarizer, the first retarder and the second retarder are arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of the transmission axis of the polarizer disagrees with the orientation of the principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of the principal axis of the first retarder, the spectroscopic polarimetry further comprises a step of acquiring data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder, and the step of obtaining the spectroscopic quasi-tokes parameter comprises: obtaining, from the obtained spectral intensity, at least one of the spectral intensity component (first spectral intensity component) which nonperiodically vibrates with wavenumber and the spectral intensity component (third spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₂(σ)) of the second retarder and not depending upon the retardation (φ₁(σ)) of the first retarder, with wavenumber, and at least one of the spectral intensity component (second spectral intensity component) which vibrates at a frequency depending upon the difference between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, the spectral intensity component (fourth spectral intensity component) which vibrates at a frequency depending upon the sum of the retardation (φ₁(σ) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, and the spectral intensity component (fifth spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₁(σ)) of the first retarder and not depending upon the retardation (φ₂(σ)) of the second retarder, with wavenumber; and obtaining the retardation (φ₁(σ)) of the first retarder, the retardation (φ₂(σ)) of the second retarder and the spectroscopic quasi-Stokes parameter by use of the data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder, and each of the obtained spectral intensity components.
 18. The spectroscopic polarimetry according to claim 16, wherein the plurality of retarders that the projection optical system comprises are a first retarder and a second retarder, the light source, the polarizer, the first retarder and the second retarder are arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder, the spectroscopic polarimetry further comprises the steps of: acquiring data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder; and acquiring a reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and a reference value (φ₂ ^((i))(σ)) for calibration of retardation of the second retarder, and the step of obtaining the spectroscopic quasi-Stokes parameter comprises: obtaining, from the obtained spectral intensity, at least one of the spectral intensity component (first spectral intensity component) which nonperiodically vibrates with wavenumber and the spectral intensity component (third spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₂(σ)) of the second retarder and not depending upon the retardation (φ₁(σ)) of the first retarder, with wavenumber, and at least one of the spectral intensity component (second spectral intensity component) which vibrates at a frequency depending upon the difference between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, the spectral intensity component (fourth spectral intensity component) which vibrates at a frequency depending upon the sum of the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder with wavenumber, and the spectral intensity component (fifth spectral intensity component) which vibrates at a frequency depending upon the retardation (φ₁(σ)) of the first retarder and not depending upon the retardation (φ₂(σ)) of the second retarder, with wavenumber; obtaining the retardation (φ₂(σ)) of the second retarder and the retardation variation (Δφ₂(σ)) of the second retarder from the reference value (φ₂ ^((i))(σ)) for calibration by use of the obtained spectral intensity; obtaining the retardation variation (Δφ₁(σ)) of the first retarder by use of the obtained retardation variation (Δφ₂(σ)) of the second retarder and data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder; obtaining the retardation (φ₁(σ)) of the first retarder from the reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and the obtained retardation variation (Δφ1(σ)) of the first retarder; and obtaining the spectroscopic quasi-Stokes parameter by use of each of the obtained spectral intensity components, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.
 19. A polarimetric spectroscope, comprising: a projection optical system, comprising a light source, a polarizer and a plurality of retarders, where the light source, the polarizer and the plurality of retarders are arranged such that light emitted from the light source is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order; an analyzer for allowing light to transmit therethrough, the light having been emitted from the projection optical system and reflected on or transmitted through the object to be measured; and a device for obtaining the spectral intensity of the light having transmitted through the analyzer.
 20. The polarimetric spectroscope according to claim 19, wherein the plurality of retarders that the projection optical system comprises are a first retarder and a second retarder, and the light source, the polarizer, the first retarder and the second retarder are arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder.
 21. The polarimetric spectroscope according to claim 20, wherein the polarizer and the second retarder are arranged such that an angle between the orientation of the transmission axis of the polarizer and the orientation of a fast axis of the second retarder is 45°.
 22. The polarimetric spectroscope according to claim 19, further comprising: an analyzer for calibration, detachably provided in a position in which light emitted from the projection optical system is received in a state where an object to be measured having an unknown spectropolarization characteristic does not exist in a light path between the projection optical system and the analyzer; and a device for obtaining the spectral intensity of the light having transmitted through the analyzer for calibration.
 23. The polarimetric spectroscope according to claim 19, further comprising an optical fiber for projecting light which guides the light emitted from the light source to the polarizer.
 24. The polarimetric spectroscope according to claim 23, wherein the device for obtaining a spectral intensity comprises a light-reception element or a spectrometer, and further comprises an optical fiber for receiving light which guides the light having transmitted through the analyzer to the light-reception element or the spectrometer.
 25. A spectroscopic polarimeter, comprising: a polarimetric spectroscope, which comprises a projection optical system, comprising a light source, a polarizer and a plurality of retarders, where the light source, the polarizer and the plurality of retarders are arranged such that light emitted from the light source is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order, an analyzer for allowing light to transmit therethrough, the light having been emitted from the projection optical system and reflected on or transmitted through the object to be measured, and a device for obtaining the spectral intensity of the light having transmitted through the analyzer; and an arithmetic unit for obtaining at least one of spectropolarization parameters of an object to be measured, by use of the spectral intensity.
 26. The spectroscopic polarimeter according to claim 25, wherein the plurality of retarders that the projection optical system comprises are a first retarder and a second retarder, the light source, the polarizer, the first retarder and the second retarder are arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder, the arithmetic unit is made capable of using data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder are obtained from the spectral intensity and the data showing the relation between the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder; and at least one of the spectropolarization parameters of the object to be measured is obtained by use of the spectral intensity, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.
 27. The spectroscopic polarimeter according to claim 25, wherein the plurality of retarders that the projection optical system comprises are a first retarder and a second retarder, the light source, the polarizer, the first retarder and the second retarder are arranged such that light emitted from the light source transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of the transmission axis of the polarizer disagrees with the orientation of the principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of the principal axis of the first retarder, the arithmetic unit is made capable of using data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder, a reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and a reference value (φ₂ ^((i))(σ)) for calibration of retardation of the second retarder, the retardation (φ₂(σ)) of the second retarder and the retardation variation (Δφ₂(σ)) of the second retarder from the reference value (φ₂ ^((i))(σ)) for calibration are obtained from the spectral intensity, the retardation variation (Δφ₁(σ)) of the first retarder is obtained by use of the obtained retardation variation (Δφ₂(σ)) of the second retarder and data showing the relation between the retardation variation (Δφ₁(σ)) of the first retarder and the retardation variation (Δφ₂(σ)) of the second retarder, the retardation (φ₁(σ)) of the first retarder is obtained from the reference value (φ₁ ^((i))(σ)) for calibration of retardation of the first retarder and the obtained retardation variation (Δφ₁(σ)) of the first retarder, and at least one of the spectropolarization parameters of the object to be measured is obtained by use of the spectral intensity, the retardation (φ₁(σ)) of the first retarder and the retardation (φ₂(σ)) of the second retarder.
 28. An optical device, comprising: a projection optical system, comprising a polarizer and a plurality of retarders, where the polarizer and the plurality of retarders are arranged such that light incident on the polarizer is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order; and an analyzer for allowing light to transmit therethrough, the light having been emitted from the projection optical system and reflected on or transmitted through the object to be measured.
 29. A light-projection device, comprising a polarizer and a plurality of retarders, wherein the polarizer and the plurality of retarders are arranged such that light incident on the polarizer is irradiated on the object to be measured after passing through the polarizer and the plurality of retarders in this order.
 30. The light-projection device according to claim 29, wherein the plurality of retarders are a first retarder and a second retarder, and the polarizer, the first retarder and the second retarder are arranged such that light incident on the polarizer transmits through the polarizer, the second retarder and the first retarder in this order, the orientation of a transmission axis of the polarizer disagrees with the orientation of a principal axis of the second retarder, and the orientation of the principal axis of the second retarder disagrees with the orientation of a principal axis of the first retarder.
 31. The spectroscopic polarimetry according to claim 1, wherein the polarimetric spectroscope prepared in the step of preparing a polarimetric spectroscope further comprises a device for changing the azimuth angle of the analyzer, the step of obtaining a spectral intensity is a step of obtaining a spectral intensity regarding the object to be measured in a plurality of states where azimuth angles of the analyzer are made different from one another, by use of the spectroscopic polarimetry, and the spectroscopic polarimetry further comprises a step of obtaining at least one of the spectropolarization parameters of the object to be measured by use of the spectral intensity obtained in the plurality of states.
 32. The spectroscopic polarimetry according to claim 1, wherein the step of preparing an object to be measured is a step of preparing an object to be measured which includes a sample and a polarization element on which light emitted from the sample is incident, the polarimetric spectroscope to be prepared in the step of preparing a polarimetric spectroscope further comprises a device for changing the characteristic of the polarization element, the step of obtaining a spectral intensity is a step of obtaining a spectral intensity regarding the object to be measured in a plurality of states where characteristics of the polarization element are made different from one another, by use of the spectroscopic polarimetry, and the spectroscopic polarimetry further comprises a step of obtaining at least one of spectropolarization parameters of the sample by use of the spectral intensity obtained in the plurality of states.
 33. The spectroscopic polarimetry according to claim 1, wherein the step of preparing an object to be measured is a step of preparing an object to be measured which includes a sample and a polarization element on which light emitted from the sample is incident, the polarimetric spectroscope to be prepared in the step of preparing a spectroscopic polarimetry further comprises a device for changing the characteristic of the polarization element and a device for changing the azimuth angle of the analyzer, the step of obtaining a spectral intensity is a step of obtaining a spectral intensity regarding the object to be measured in a plurality of states where characteristics of the polarization element, or azimuth angles of the analyzer, are made different from one another, by use of the spectroscopic polarimetry, and the spectroscopic polarimetry further comprises a step of obtaining at least one of spectropolarization parameters of the sample by use of the spectral intensity obtained in the plurality of states.
 34. The polarimetric spectroscope according to claim 19, further comprising a device for changing the azimuth angle of the analyzer.
 35. The spectroscopic polarimetry according to claim 19, further comprising a device for changing the characteristic of the polarization element in a case where the object to be measured includes a sample and a polarization element on which light emitted from the sample is incident.
 36. The spectroscopic polarimetry according to claim 19, further comprising a device for changing the characteristic of the polarization element, and a device for changing the azimuth angle of the analyzer, in a case where the object to be measured includes a sample and a polarization element on which light emitted from the sample is incident.
 37. The optical device according to claim 28, further comprising a device for changing the azimuth angle of the analyzer.
 38. The optical device according to clam 28, further comprising a device for changing the characteristic of the polarization element in a case where the object to be measured includes a sample and a polarization element on which light emitted from the sample is incident.
 39. The optical device according to claim 28, further comprising a device for changing the characteristic of the polarization element, and a device for changing the azimuth angle of the analyzer, in a case where the object to be measured includes a sample and a polarization element on which light emitted from the sample is incident. 